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Preface 


The  purpose  of  this  study  was  to  develop  a  microcomputer- 
based  engineering  workstation  for  the  analysis  and  design  of 
digital  satellite  links.  While  there  is  lots  of  software 
for  performing  link  budget  analysis  already  on  many 
mainframe  computers,  it  is  not  particularly  convenient  for 
engineers  doing  real-time  design  work.  There  is  also  some 
software  available  for  personal  computers,  but  it  exists 
largely  as  a  poorly-documented  hodgepodge  of  machine- 
dependent  code  developed  by  various  organizations  for  their 
own  specific  purposes.  I  perceived  a  need  for  a 
transportable  package  of  link  analysis  software  which  could 
run  on  a  wide  variety  of  personal  computers.  I  decided  to 
design  my  software  for  the  IBM  PC,  since  it  represents  by 
far  the  largest  class  of  compatible  machines. 

As  far  as  coding  goes,  I  needed  a  programming  language 
that  had  a  good  set  of  color  graphics  commands,  as  well  as 
one  that  could  be  obtained  at  low  cost.  The  obvious  choice 
was  Turbo  Pascal. 

I  would  like  to  thank  Major  Dale  Hibner  for  his  valuable 
advice  and  assistance  in  the  planning  and  development  of 
this  rather  massive  undertaking.  I  also  wish  to  thank  my 
wife  Pookie  for  graciously  allowing  me  to  put  the  color 
monitor  on  VISA.  Now,  about  that  Winchester  drive . 


William  C.  Jackson 
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ABSTRACT 

'  J_  'nZSt* 

The  product  of  this  -research  -gffort  is  a  comprehensive 

and  integrated  software  package  for  the  analysis  of  digital 
satellite  links.  The  system  is  designed  to  run  on  the  IBM 
PC,  and  should  also  run  on  most  compatibles. 

The  system  performs  three  basic  classes  of  functions: 
satellite  orbital  analysis,  antenna  gain  pattern  plotting, 

and  link  analysis.  The  first  class^f  functions  includes 

I 

the  computation  of  such  quantities  as  velocity,  orbital 
period,  and  coverage  area  for  satellites  in  circular  and 
elliptical  orbits.  The  second  class /of  functions  is 
concerned  with  plotting  the  gain  patterns  for  horns, 
helixes,  parabolic  reflectors,  and  phased  arrays  of 
dipoles.  The  last  class  of  function's"  represents  the  major 
thrust  of  the  system,  and  entails  computing  such  items  as 
the  G/T  figure  of  merit,  received  useful  power,  carrier-to- 
noise  ratio,  bit  error  rate,  maximum  data  rate,  and  power 
margin.  Inherent  within  this  class  are  mathematical  models 
for  computing  the  attenuation  due  to  rainfall  and 
atmospheric  absorption. 

The  link  budget  itself  appears  as  a  color-coded  display 
with  two  columns:  one  for  the  uplink  path,  and  one  for  the 
downlink  path.  The  user  also  has  the  capability  to  change 
certain  key  inputs,  and  then  have  the  system  automatically 
recompute  the  entire  link  budget  with  the  modified  data. 


DEVELOPMENT  OF  A  MICROCOMPUTER-BASED  WORKSTATION  FOR 
ANALYSIS  OF  DIGITAL  SATELLITE  LINKS 

I.  INTRODUCTION 

Orbiting  communications  satellites  play  a  vital  role  in 
information  transmission  well  beyond  line-of-sight .  One 
frequently  hears  the  satellite  transponding  process 
described  as  "bouncing”  a  signal  off  of  a  satellite.  This 
is  really  a  misnomer,  since  modern  communications  satellites 
are  a  far  cry  from  the  early  ECHO  series  passive  reflector 
satellites.  The  transponding  process  is  actually  an  active 
one,  and  consists  of  two  distinct  components:  the  uplink  and 
the  downlink.  The  uplink  signal  is  transmitted  from  a 
ground  station  to  the  satellite,  where  it  is  translated  in 
frequency  and  re-transmitted  as  a  downlink  signal  to  another 
ground  station.  The  link  itself  can  best  be  described  as 
the  total  transmission  path  from  the  first  ground  station 
transmitter,  through  the  channel,  to  the  receiver  at  the 
second  ground  station.  This  includes  all  the  modulators, 
transmitters,  receivers,  demodulators,  and  antennas  both  on 
the  ground  and  onboard  the  satellite.  The  channel  is  the 
total  electromagnetic  path  connecting  a  given  transmitter- 
receiver  pair,  and  includes  the  atmosphere  as  well  as  free 
space . 

In  a  typical  satellite  communications  link,  the  received 
signal  may  degrade  in  one  of  two  ways:  through  a  loss  in 
desired  waveform  power,  or  through  the  addition  of  noise  to 
the  signal.  Losses  occur  when  a  portion  of  the  signal  is 


diverted  or  reflected  from  its  intended  route,  or  is 
attenuated  along  its  intended  route.  Noise  arises  from 
unwanted  signal  energy  being  injected  into  the  link,  or  from 
thermal  noise  generated  within  the  link.  In  a  typical  link, 
free  space  propagation  loss  is  the  single  largest  loss  in 
signal  power.  Prominent  sources  of  noise  include  the 
internal  noise  of  receiver  amplifiers  and  feed  lines,  as 
well  as  the  noise  introduced  by  receiving  antennas. 

One  of  the  most  important  tools  available  to  the  engineer 
involved  in  the  design  and  analysis  of  communications 
satellite  links  is  link  budget  analysis.  A  system  link 
budget  is  essentially  a  balance  sheet  of  gains  and  losses. 

It  is  comprised  of  a  detailed  listing  of  transmission  and 
reception  sources,  noise  sources,  and  signal  attenuators  as 
measured  along  the  entirety  of  the  link.  A  good  part  of  the 
link  budget  is  derived  from  the  calculation  of  received 
useful  power.  Link  budget  analysis  is  basically  an 
estimation  technique  for  evaluating  communications  system 
performance.  Link  budgets  are  especially  useful  for  the 
following  (10:1-3): 

1.  Determining  hardware  constraints 

2.  Predicting  system  performance 

3.  Discovering  system  design  flaws 

4.  Experimenting  with  various  design  tradeoffs 

5.  Predicting  system  availability 

6.  Highlighting  system  nuances 

7.  Searching  for  an  optimal  design 


The  primary  purpose  of  this  research  effort  is  to  develop 
a  comprehensive  set  of  software  routines  for  performing  link 
budget  analysis,  and  to  integrate  them  into  a  transportable 
package  for  use  on  a  personal  computer.  Secondary  goals 
include  developing  software  to  compute  various  satellite 
orbital  parameters,  as  well  as  software  to  plot  antenna  gain 
patterns.  This  software  in  effect  turns  a  personal  computer 
into  a  dedicated  satellite  link  analysis  workstation. 

It  is  assumed  that  the  user  has  at  least  a  cursory 
knowledge  of  digital  satellite  communications,  and  is 
familiar  with  link  budgets.  While  this  project  might  well 
be  useful  as  a  learning  tool,  it  is  not  intended  to  be  a 
tutorial  for  the  novice.  No  particular  knowledge  of  the 
software  is  required  to  successfully  run  it,  since  it  is 
largely  menu  driven.  However,  a  user’s  manual  is  included 
as  Appendix  A  for  those  seeking  more  detailed  information. 

It  is  also  assumed  that  the  user's  machine  is  an  IBM  PC  or 
compatible  with  128K  of  memory,  at  least  one  floppy  disk 
drive,  and  a  color  monitor. 

Several  other  technical  assumptions  were  made  regarding 
the  links 

1.  All  noise  is  thermal  additive  white  gaussian  noise, 
with  a  flat  power  spectral  density  of  Nq  for  all 
frequencies  of  interest. 

2.  Signal  bandwidth  is  taken  to  be  equal  to  noise 
bandwidth . 

3.  The  mathematical  models  for  computing  atmospheric 
and  rainfall  attenuation  neglect  such  ubiquitous 


factors  as  polarization  losses,  ice  crystal  effects, 
multipath  propagation,  and  atmospheric  scintillation. 

4.  No  provision  is  made  for  error-correcting  codes. 

5.  Pulse  shaping  techniques  are  not  taken  into 
account,  so  information  transmission  is  subject  to  the 
constraints  dictated  by  the  Nyquist/ Shannon  criteria. 

In  a  software  project  of  this  magnitude,  it  is  necessary 
to  use  a  highly  structured  design  methodology.  The  first 
step  is  to  establish  an  initial  set  of  system  requirements, 
recognizing  that  these  requirements  are  fluid,  and  subject 
to  updating  as  development  proceeds.  It  is  important  to 
generate  specific  and  meaningful  requirements  before  any 
software  design  takes  place,  so  that  the  requirements  will 
drive  the  design,  and  not  vice-versa.  Once  this  is  done, 
the  overall  system  architecture  must  be  defined.  This  step 
establishes  the  primary  conceptual  issues  upon  which  the 
rest  of  the  system  is  based.  In  addition,  the  software 
modules  comprising  the  system  are  identified.  Next,  the 
detailed  design  is  accomplished,  where  the  software  modules 
themselves  are  designed.  It  is  at  this  level  that 
fundamental  algorithms  and  data  relationships  are  defined. 
Finally,  the  completed  software  must  be  validated  through  a 
comprehensive  testing  procedure. 

While  the  design  is  performed  in  a  top-down  manner,  the 
software  itself  is  coded  in  a  bottom-up  manner.  The  general 
methodology  is  that  the  lowest  level  modules  are  written 
first,  followed  by  any  higher-level  modules  (modules  which 


call  other  modules).  The  main  program  is  written  last,  and 
the  completed  modules  are  integrated  into  it.  Module 
testing  may  occur  before  or  after  the  integration  phase. 

The  final  system  check  is  performed  with  all  modules  tested 
and  in  place  in  the  main  program. 


II.  THEORY  OF  SATELLITE  COMMUNICATIONS 


This  chapter  presents  a  brief  overview  of  satellite 
communications  theory,  to  include  selected  aspects  of  space 
communications  as  well  as  orbital  mechanics.  This  chapter 
is  not  intended  as  a  tutorial  on  the  principles  of 
electronic  communications  per  se,  but  rather  as  an 
introduction  to  the  anatomy  of  a  satellite  communications 
system. 

Link  Description 

For  the  purposes  of  this  paper,  a  satellite 
coimnunication8  system  consists  of  two  earth  stations 
separated  by  some  distance,  and  an  orbiting  satellite.  The 

complications  of  satellite-to-satellite  communications  or 
multiple  users  are  not  considered.  Since  the  satellite 
transmits  the  downlink  by  responding  to  the  uplink,  it  is 
referred  to  as  a  transponder.  If  this  transponding  process 
were  accomplished  passively,  then  the  downlink  power  level 
would  be  extremely  low  due  primarily  to  the  large  path 
losses.  An  active  satellite  transponder  assists  this 
process  by  adding  power  amplification  to  the  signal  prior  to 
downlink  transmission.  However,  transmitting  the  downlink 
at  the  same  frequency  as  the  uplink  would  cause  unwanted 
feedback,  so  the  transponder  must  perform  some  sort  of 
frequency  translation  prior  to  amplification  and  re¬ 
transmission.  Separation  of  the  uplink  and  downlink 
frequency  bands  also  allows  the  same  antenna  to  be  used  for 
both  receiving  and  transmitting  (2t7-8). 

6 


System  Components.  The  behavior  of  a  satellite  link  is 
best  understood  by  analyzing  each  of  its  components,  the 
term  "component"  being  taken  to  mean  a  source  of  gain,  loss, 
or  noise.  Each  component  belongs  to  one  of  three  sections 
of  the  link:  the  transmitter,  channel,  or  receiver.  Figure 
2-1  illustrates  the  major  components  in  a  typical  satellite 
link.  This  figure  could  represent  either  an  uplink  or  a 
downlink.  The  transmitter  section  consists  of  a  high-power 
microwave  transmitter  and  an  antenna,  with  some  type  of 
cable  or  waveguide  between  the  two.  The  circuit  losses 
include  the  effects  due  to  both  the  lossy  line  and 
imperfectly  matched  couplings.  The  efficiency  of  the 
antenna  takes  into  account  losses  due  to  reradiation, 
scattering,  and  spillover. 

The  channel  section  consists  of  free-space  path  losses, 
atmospheric  losses,  and  rainfall  losses.  Path  loss  is  due 
to  the  decrease  in  electric  field  strength  with  distance, 
and  is  unavoidable.  Atmospheric  losses  are  due  to 
absorption  of  the  signal  by  atmospheric  oxygen  and  water 
vapor.  Rainfall  losses  are  due  to  the  attenuation  of  the 
signal  by  rain  at  the  uplink  or  downlink  site,  or  both.  The 
effects  of  atmospheric  and  rainfall  attenuation  can  be 
minimized  by  an  appropriate  choice  of  frequency. 

The  receiver  section  consists  of  the  receiving  antenna,  a 
low-noise  receiver,  and  some  type  of  coupling  between  the 
two.  This  side  of  the  link  is  similar  to  the  transmitter 
side,  with  some  additional  complications.  The  receiving 
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Figure  2-1.  Major  Components  of  a  Typical  Satellite  Link 
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antenna,  noise  temperature  now  contributes  to  system  noise. 
Sources  of  antenna  noise  may  include  galactic  noise,  man¬ 
made  noise,  as  well  as  atmospheric  absorption.  For  earth 
coverage  satellite  antennas,  290°K  is  generally  used  as  the 
antenna  temperature;  for  ground  station  antennas  it  may  be 
much  less,  since  the  antenna  is  looking  out  toward  a  cooler 
sky  temperature.  Receiver  noise  is  also  a  fundamental 
source  of  system  noise.  For  this  reason,  the  receiver 
amplifier (s)  are  sometimes  cryogenically  cooled  to  lower 
their  noise  temperature  (10sl-3). 

Design  Considerations.  Because  size  and  weight  are  not 
limiting  factors,  ground  station  transmitters  may  be 
arbitrarily  large.  Power  outputs  of  40-60  dBW  are  common. 

In  addition,  the  antenna  transmitting  the  uplink  may  also  be 
quite  large.  Large  earth  stations  may  employ  parabolic 
dishes  with  diameters  of  13-19  m.  Conversely,  satellite 
transmitters  and  antennas  are  subject  to  severe  size  and 
weight  constraints.  Spaceborne  transmitter  powers  may  be 
less  than  10  W;  spaceborne  parabolic  dishes  are  typically 
less  than  one  meter  in  diameter. 

Uplink  and  downlink  frequencies  are  not  chosen 
arbitrarily.  Above  10  GHz,  atmospheric  attenuation  is  a 
growing  concern.  Satellite  communications  bands  are 
specifically  allocated  to  avoid  the  resonance  absorption  due 
to  water  vapor  at  22  GHz,  and  oxygen  at  about  60  GHz.  In 
general,  it  is  easiest  and  most  economical  to  use  the  lowest 
available  frequencies.  However,  bandwidth  is  limited  in  the 


$ 
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lower  frequency  bands  due  to  the  large  number  of  other 
users,  and  interference  is  more  likely  (2:106-112;  9:104- 
105) . 


Signal  Analysis 

In  order  to  understand  the  concept  of  satellite  link 
analysis,  it  is  necessary  to  analyze  the  signals  at  various 
points  in  the  link.  The  behavior  of  signals  in  the 
transmitter,  antenna,  channel,  and  receiver  is  examined 
separately,  and  various  system  performance  criteria  are 
given . 

Transmitter  (2:83).  The  transmitter  field  is 
characterized  by  its  effective  isotropic  radiated  power 
(EIRP),  which  is  given  by 


EIRP 


P  +  TxGain  +  TxCircuitLoss 
out. 


(1) 


where  P  .-transmitter  power  (dB) 
out 


TxGain-transmitting  antenna  gain  (dB) 
TxCircuitLoss-transmitter  circuit  losses  (dB) 


Antenna  (2:87-90).  The  antenna  is  a  key  element  of  a 
satellite  communication  system  since  its  gain  has  a  direct 
effect  on  the  amount  of  received  power.  An  antenna  may  be 
specified  in  terms  of  its  gain  and  3-dB  beamwidth.  Gain  is 
proportional  to  the  square  of  the  frequency  and  to  the 
antenna  size,  while  beamwidth  varies  inversely  to  frequency 
and  size.  Formulas  for  the  gains  and  beamwidths  of  several 
different  antenna  types  are  given  below: 
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where  d«antenna  dimension 
L-wavelength 
c*helix  circumference 
l*helix  length 
N-nuraber  of  elements 
8*element  separation  (2*89) 


Channel  (2 t 15-16, 84-85) .  The  most  significant  source  of 


signal  loss  in  the  channel  is  free-space  propagation  or  path 


loss,  which  is  given  by 


Path  Loss  -  -36.6  -  20  logCSf] 


where  S«slant  range  (mi) 
f»frequency  (MHz) 


The  slant  range  itself  is  computed  from 


(  R  2sin2a  +  2R  h  +  h2  )1/2  -  R  sin  a 
e  e  e 


where  Re»radius  of  the  Earth 


a*terminal  elevation  angle 
h*satellite  altitude 


Atmospheric  attenuation  and  rainfall  also  contribute  to 


total  channel  losses,  and  are  discussed  in  detail  in 


Appendices  B  and  C. 


Receiver  (9i 113-118).  The  amount  of  receiver  noise 


present  is  given  by  the  system  equivalent  temperature.  This 


is  the  effective  equivalent  temperature  that  an  external 
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noise  source  would  need  to  have  in  order  to  produce  the  same 
amount  of  receiver  noise.  It  is  given  by 

Teq  *  RxAntTemp  +  [ (F/RxCircuitLoss)  -  1]*290  (8) 

where  RxAntTemp»receiving  antenna  temperature  (deg  K) 
F=noise  figure 

RxCircuitLosssreceiver  circuit  losses  (W) 

A  receiver  may  also  be  characterized  by  its  G/T  figure  of 
merit,  which  is  given  by 

G/T  *  RxGain/T  (9) 

eq 

where  RxGain»receiving  antenna  gain  (W) 

T  “system  equivalent  temperature  (deg  K) 
eq 

System  Performance  (2:104-105;  9:110-113).  There  are 
several  criteria  which  indicate  the  quality  of  the  link. 
These  include  received  useful  power  (P^.),  carrier  power-to- 
noise  power  density  ratio  (CNR),  bit  energy-to-noise  power 
density  ratio  (Eb/NQ),  bit  error  rate  (BER),  and  maximum 
data  rate . 

Received  useful  power  is  given  by 

Pr  *  p0ut  +  TxCircuitLoss  +  TxGain  +  TxEff  + 

PathLoss  +  AirAtten  +  RainAtten  +  RxEff  + 

RxGain  +  RxCircuitLoss  (10) 

where  TxEff»transmitting  antenna  efficiency 
RxEf f«receiving  antenna  efficiency 
PathLoss“loss  due  to  propagation  over  a  distance 
AirAtten-atmospheric  attenuation 
RainAtten*attenuation  due  to  rainfall 
(all  units  are  dB) 

The  carrier  power-to-noise  power  density  ratio  is  given  by 
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where  Pr=received  useful  power  (W) 
k*Boltzman 1  a  constant 

Te^*system  equivalent  temperature  (deg  K) 
Available  is  given  by 

VNo  *  Pr/RNo  <12> 


where  Pr=received  useful  power  (W) 

Rj=data  rate  (bits/sec) 

No=noise  power  density  (W/Hz) 

The  bit  error  rate  is  determined  from  the  modulation  type 

and  available  Ej/N  •  The  expressions  for  bit  error  rate  for 

eleven  different  modulation  formats  are  listed  below,  where 


Q(x)  is  Marcum*  s  Q-function  (2:53-68). 

BPSK,  QPSK,  OQPSK,MSK:  Q(  [2^^]  *  5 )  (13) 
MPSK:  Q([21og2M  Eb/NQ]*5sin  pi/M)/log2M  (14) 
Coherent  FSKs  Q([2E./N  ]*5)  (15) 
Non -Coherent  FSK:  . 5exp(-Eb/2No)  (16) 
MFSKi  M/4  exp(-log2M  Eb/NQ)  (17) 
ASK:  Q([Eb/N0]*5)  (18) 
MASK:  C(M-1)/M]  Q{[61og2M;j/[  (M2-l  JE^N^  }/log2M  (19) 
DPSK:  . 5exp(-Eb/NQ)  (20) 


By  algebraically  manipulating  the  above  expressions,  the 
Eh/N0  required  for  a  given  bit  error  rate  and  modulation 
type  may  be  found;  these  new  expressions  now  involve  inverse 
Q-functions.  The  difference  between  the  available  Ej^/^  and 
the  required  Ej^/^  is  called  the  power  margin.  Finally,  the 
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maximum  data  rate  is  given  by 


maxrate  *  Pr/(kTeq  E^l^)  (21) 

where  Pr*»received  power  (W) 

kTeg*noise  power  density  (W/Hz) 

Ej^/N^^-required  bit  energy-to-noise  power  density 
ratio 


Orbital  Aspects 


This  section  presents  some  of  the  salient  physical 
aspects  of  orbiting  satellites#  beginning  with  a  brief 
discussion  of  orbital  mechanics*  In  addition#  a  methodology 


for  determining  look  angles  is  also  presented. 

Orbital  Mechanics  (9 t 11-22).  In  order  to  be  placed  into 
orbit#  a  satellite  must  escape  the  influence  of  the  Earth's 
gravitational  field.  This  escape  velocity  is  a  function 


only  of  altitude,  and  is  given  by 


Vesc  *  C2GMe  /  (^  +  h)]*5 


(22) 


where  G=Earth  gravitational  constant 

h**altitude  of  launch  site 

M=mass  of  Earth 
e 

R  =radius  of  Earth 
e 

Once  in  space,  the  satellite  is  placed  into  either  a 
circular  or  an  elliptical  orbit.  For  the  case  of  a  circular 
orbit#  the  satellite  must  achieve  a  velocity  of 

v  -  [u  /  (Re  +  h)]*5 

where  u«GM  (gravitational  coefficient) 

®  14 


(23) 


The  orbital  period,  or  time  to  complete  one  revolution, 
is  given  by 

T  =  2*pi*[Re  +  h]1,5/u*5  (24) 

Another  important  consideration  is  coverage  area.  This 
refers  to  the  area  on  the  Earth  from  which  the  satellite  is 
visible  from  some  pre-determined  minimum  elevation  angle. 
Refer  to  the  set-up  depicted  in  Figure  2-2.  The  subtended 
angles  may  be  determined  from  simple  trigonometry  to  be 

s  =  cos”1 [R  cos  El  /  (R  +  h)]  -  El  (25) 

where  h=satellite  altitude 

El=minimum  elevation  angle 

The  coverage  area  is  therefore  the  area  of  the  spherical  cap 
subtended  by  the  angle  s,  and  is  given  by 

Acov  =  (2*pi*Re2)(l  -  cos  s)  (26) 

The  maximum  length  of  time  that  an  orbiting  satellite  is 
visible  to  an  observer  on  the  Earth's  surface  is  referred  to 
as  the  satellite's  duration  of  visibility.  Duration  of 
visibility  is  a  function  of  subtended  angle  and  orbital 
period.  It  is  expressed  as 

tmax  =  C2s/360]CT/(1  ±  T/TJ]  (27) 

where  Te*Earth's  orbital  period 

+  *  retrograde  orbit 
-  ”  prograde  orbit 

Note  that  for  the  special  case  of  a  prograde  geostationary 
orbit,  the  satellite  is  visible  continuously. 


For  the  case  of  an  elliptical  orbit,  the  satellite 
velocity  is  not  constant,  but  varies  with  position. 
Velocity  is  highest  at  the  orbital  perigee  (low  point  of 
orbit),  and  lowest  at  the  orbital  apogee  (high  point  of 
orbit).  The  major  and  minor  axes  a  and  b  of  the  orbital 
ellipse  are  given  respectively  by 


a  =  (R  R  ) 
p  a ' 

b  =  (R  +  R  )/2 
p  a " 


where  Rp=altitude  at  perigee 
R  =altitude  at  apogee 

cl 


(28a) 

(28b) 


From  these,  the  eccentricity  of  the  orbit  may  be  expressed 


E  =  (1  -  b2/a2 ) *  5 


The  polar  equation  for  the  orbit  is  then 


r  =  a(l  -  E  )/(l  +  E  cos  0) 


(29) 


(30) 


where  r=radial  distance  from  geometric  center  of  ellipse 
(Preference  angle  (zero  at  perigee) 

The  satellite's  velocity  at  any  given  radial  distance  r 


is  given  by 


v  =  [u( 2/r  -  1/a)] 


The  orbital  period  is  given  by 
T  -  2*pi*a1 *  5/u*  5 


(31) 


(32) 
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The  formula  for  coverage  area  is  similar  to  the  one  given 
for  a  circular  orbit,  with  the  exception  that  the  subtended 
angle  is  now  given  by 

s  =  cos_1[Re  cos  El  /  (Rfi  +  r)]  -  El  (33) 

Look  Angle  Determination  (9:22-30).  The  coordinates  to 
which  an  earth  station  antenna  must  point  in  order  to 
communicate  with  a  satellite  are  referred  to  as  the  look 
angles.  They  are  usually  specified  as  azimuth  and  elevation 
angle.  The  geometry  of  the  problem  is  given  in  Figure  2-3. 
Given  the  coordinates  of  the  earth  station  (lat, long)  and 
the  coordinates  of  the  sub-satellite  point  (lat  .long  ),  the 
great  circle  distance  between  the  two  is  given  by 

R  =  cos-1 [cos (lat  -  lat  )  cos (long  -  long  )]  (34) 

s  s 

The  azimuth  is  given  by 

Az  *  tan  ^[tandong  -  long  )  /  sin(lat  -  lat  )]  (35) 

The  tilt  angle  is  the  angle  between  the  satellite  vertical 
to  the  Earth  and  the  line-of-sight  to  the  ground  station. 
Normalizing  the  satellite's  altitude  h  with  respect  to  the 
Earth's  radius,  one  may  define  p*»l+h/Re.  The  tilt  angle  is 
then  expressed  as 

Tilt  *  tan_1[sin  R  /  (p  -  cos  R)]  (36) 

Finally,  the  elevation  angle  may  be  expressed  in  terms  of 
the  tilt  angle  by 


V’ 


Figure  2-3.  Geometry  For  Determining  Look  Angles 
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El  ■  cos-1[p  sin(Tilt)]  (37) 

This  chapter  has  presented  the  user  with  the  basic  theory 
needed  to  understand  satellite  link  analysis.  The  major 
components  of  a  typical  satellite  link  were  described 
qualitatively,  and  several  important  design  considerations 
were  discussed.  In  addition,  the  behavior  of  the  signal  in 
each  principal  section  of  the  link  was  described 
quantitatively,  and  several  system  performance  criteria  were 
stated.  Finally,  a  brief  synopsis  of  orbital  mechanics  was 
presented,  along  with  a  method  of  determining  look  angles. 


III.  SYSTEM  REQUIREMENTS 


This  section  addresses  the  specification  of  required 
functions/  interfaces/  and  performance  of  the  end  software 
product.  The  intent  is  to  present  a  complete  and 
unambiguous  set  of  functional  specifications  describing  what 
the  software  is  required  to  do,  rather  than  how  it  is  to  be 
implemented  at  the  coding  level.  This  phase  of  software 
development  is  crucial,  because  inadequate  specification  of 
requirements  often  leads  to  system  deficiencies  which  are 
difficult  to  correct  at  a  later  date.  In  addition,  system 
testing  tends  to  have  little  real  meaning  without  a  concrete 
set  of  requirements  as  a  basis. 

The  essential  elements  of  a  good  requirements  definition 
are  testability  and  flexibility.  This  means  that  it  must  be 
possible  to  verify  that  the  end  product  satisfies  the  given 
specifications.  Also,  the  specifications  must  not  be  so 
rigid  as  to  inhibit  program  development. 

Functional  Description 

The  following  is  a  global  description  of  the  system  in 
terms  of  the  functions  that  it  must  accomplish.  No  attempt 
will  be  made  here  to  derive  equations  or  develop  algorithms, 
as  this  is  merely  a  statement  of  intended  capability.  For 
the  purposes  of  this  chapter,  required  functions  are 
classified  as  primary,  secondary,  or  housekeeping 
functions.  Primary  functions  are  those  which  represent  the 
major  thrust  of  the  thesis,  and  include  antenna  analysis. 


<*  receiver/transmitter  analysis,  and  various  other  system 

design  calculations.  Secondary  functions  are  those  which 
supplement  the  overall  system,  but  are  not  required  for 
system  operation,  such  as  calculation  of  satellite  orbital 
parameters  and  plotting  antenna  gain  patterns.  Housekeeping 
functions  are  those  which  perform  no  calculations,  but  are 
responsible  for  systems-level  tasks  such  as  storing  and 
manipulating  data. 

Finally,  design  constraints  such  as  portability, 
expandability,  and  audit  capability  are  discussed.  Note 
that  design  constraints  specify  how  the  system  is  to  be 
constructed  and  implemented,  not  necessarily  which  functions 
^  are  to  be  included  in  the  system. 

Primary  Functions.  As  a  minimum,  this  set  performs  the 
functions  listed  below.  Note  that  the  first  item  listed 
applies  equally  to  horns,  helixes,  parabolic  dishes,  or 
phased  arrays  of  dipoles. 

1.  Compute  the  maximum  antenna  gain  and  half-power 

beamwidth . 

2.  Compute  the  effective  isotropic  radiated  power. 

3.  Compute  the  total  system  noise  temperature  and  G/T 

figure  of  merit. 

4.  Calculate  the  attenuation  due  to  rainfall  and 


atmospheric  absorption. 

5.  Calculate  received  useful  power  as  governed  by  the 

£  link  equation. 

6.  Compute  the  carrier  power- to-noise  power  density  ratio 
(CNR) . 
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7.  Determine  the  maximum  data  rate  as  a  function  of  CNR, 
desired  bit  error  rate,  and  type  of  modulation. 

8.  Compute  the  bit  error  rate  given  the  CNR,  data  rate, 
and  type  of  modulation. 

9.  Compute  the  link  margin. 

Secondary  Functions.  This  set  of  functions  should  be 
considered  optional,  and  has  absolutely  no  bearing  on  link 
budget  calculations.  These  functions  concern  satellite 
orbit  parameters  and  antenna  gain  pattern  plotting,  and 
perform  the  following  i 

1.  Calculate  escape  velocity  at  a  given  altitude  above 
the  Earth's  surface. 

2.  Compute  the  steady-state  temperature  of  a  satellite  in 
Earth  orbit. 

3.  Determine  velocity,  orbital  period,  coverage  area,  and 
duration  of  visibility  for  satellites  in  circular  orbit. 

4.  Determine  velocity  at  apogee  and  perigee,  average 
velocity,  orbit  eccentricity,  orbital  period,  and 
coverage  area  for  satellites  in  elliptical  orbit. 

5.  Compute  the  azimuth  and  elevation  angles  to  a 
satellite  from  a  given  ground  site. 

6.  Plot  the  antenna  gain  patterns  for  horns,  helixes, 
parabolic  dishes,  and  phased  arrays  of  dipoles. 
Housekeeping  Functions.  These  functions  are  primarily 

concerned  with  disk  access  and  database  management.  A 
capability  is  included  to  store  and  retrieve  files  of  data. 


as  well  as  to  delete  unwanted  files.  The  user  is  also 


provided  the  capability  to  change  the  default  drive  to  which 
data  is  written,  and  to  review  the  files  residing  on  a  given 
data  disk.  These  functions  are  designed  to  mimic  the 
corresponding  functions  in  DOS. 

Design  Constraints.  There  are  several  design  constraints 
which  drastically  affect  a  programming  effort  of  this  type. 
The  first  is  that  of  portability.  Because  the  resultant 
software  must  run  on  a  variety  of  machines,  it  is  necessary 
to  incorporate  a  certain  amount  of  flexibility  in  the  early 
stages  of  system  design  in  order  to  avoid  the  problem  of 
machine  dependence.  All  code  is  written  in  TURBO  Pascal, 
which  is  available  in  both  MS-DOS  and  CP/M  versions,  and 
therefore  can  run  on  a  wide  variety  of  microcomputers.  For 
this  reason,  one  copy  of  the  program  is  left  as  source  code, 
in  order  to  facilitate  transportability  from  an  MS-DOS 
machine  to  a  CP/M  machine  and  vice-versa.  The  only  possible 
problem  area  (with  regard  to  machine  dependence)  is  that  of 
display  graphics.  In  order  to  circumvent  this  problem,  all 
graphics  are  generated  by  a  graphics  driver  module.  While 
this  module  itself  is  highly  machine  dependent,  it  can  be 
easily  removed  intact  and  replaced  by  another  custom  driver 
module  for  a  different  machine. 

The  second  constraint  is  that  of  expandability.  Since 
system  requirements  tend  to  be  in  a  constant  state  of  flux 
in  a  project  of  this  nature,  it  is  important  not  to  make  the 
system  design  too  rigid.  In  other  words,  as  the  need  for 
new  functions  arises,  one  should  not  have  to  completely 


revamp  the  entire  architecture  in  order  to  accommodate  it* 
The  key  to  expandability  is  therefore  modular  design.  This 
allows  logical  units  to  be  inserted  and  deleted  at  will. 

The  last  constraint  to  be  addressed  is  that  of  audit 
capability.  This  refers  to  the  calculation  and  storage  of 
intermediate  results  during  data  entry  or  calculations. 
Incorporating  an  audit  capability  serves  several  purposes. 
First,  it  facilitates  system  testing  by  storing  an  audit 
trail  of  both  module  function  and  module  interaction. 
Secondly,  it  permits  the  user  to  go  back  through 
intermediate  calculations  when  seemingly  erroneous  results 
are  obtained  and  discover  where  the  error  occurred. 

Finally,  in  the  unlikely  event  of  a  program  crash, 
significant  data  entered  prior  to  that  point  will  be  stored 
in  memory  for  subsequent  recovery. 

Man-Machine  Interface 

The  success  of  any  interactive  software  design  effort  is 
measured  by  how  well  the  system  insulates  the  user  from 
technical  details  of  the  software.  While  the  programmer 
must  concern  himself  with  system  architecture,  data 
structure  and  flow,  and  control  structure  and  flow,  the  user 
is  only  concerned  with  these  issues  as  far  as  they  affect 
operations  at  the  terminal.  As  far  as  the  user  is 
concerned,  the  system  is  a  black  box.  This  section 
addresses  some  of  the  issues  encountered  in  designing  such  a 
system. 


User  Interaction.  A  variety  of  user-friendly  features 


are  incorporated  into  a  typical  interactive  software 
system.  The  most  obvious  feature  is  the  use  of  menus.  A 
menu  is  simply  a  display  of  a  set  of  program  options,  from 
which  the  user  is  asked  to  select  at  least  one.  This  means 
that  the  user  is  able  to  move  freely  through  the  hierarchy 
of  the  system  simply  by  entering  a  response  to  a  list  of 
options.  The  use  of  menus  eliminates  guesswork  on  the  part 
of  the  user,  and  consolidates  all  necessary  user  decisions 
into  convenient,  easy  to  use  subsets. 

Closely  associated  with  menus  is  the  use  of  prompting  and 
validity  checking.  Prompting  is  the  use  of  printed  test 
messages  asking  the  user  to  enter  some  specific  data,  or  to 
perform  some  other  particular  function.  This  eliminates 
confusion  as  to  exactly  what  data  is  required,  what  format 
the  data  should  be  in,  and  what  units  (if  any)  the  data 
should  be  in.  Validity  checking  refers  to  the  practice  of 
subsequently  analyzing  user-entered  data  and  ensuring  that 
it  conforms  to  some  set  standard  (e.g.  within  a  specified 
range  of  validity).  If  the  data  fails  the  validity  check, 
then  the  user  is  alerted  and  asked  to  re-enter  it.  This 
ensures  that  all  data  going  in  to  the  program  is  valid. 

At  this  point  something  should  be  said  regarding 
development  of  the  user's  guide.  During  the  initial  stages 
of  system  planning  and  design,  the  programmer  must 
constantly  make  tentative  decisions  and  assess  their  impact 
on  the  overall  system.  As  time  progresses  and  previously 
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unseen  constraints  and  limitations  appear,  it  may  become 
necessary  to  change  earlier  decisions.  The  programmer  is 
now  falling  into  the  trap  of  letting  the  code  dictate  system 
operation.  This  pitfall  can  be  avoided  by  developing  a 
preliminary  version  of  the  user's  guide  in  parallel  with  the 
design.  This  practice  is  a  form  of  self-imposed  programming 
discipline,  and  ensures  that  the  original  system 
requirements  will  in  fact  drive  the  implementation. 

Visual  Display  of  Data.  A  significant  part  of  the  man- 
machine  interface  is  the  manner  in  which  information  is 
presented  to  the  user.  One  aspect  of  this  is  the  use  of 
color.  Used  intelligently,  color  can  greatly  enhance  the 
readability  of  both  text  and  data.  For  example,  a  table 
consisting  of  several  different  sets  of  data  might  use  a 
different  color  for  each  set  in  order  to  more  easily 
differentiate  them.  With  a  standard  RGB  color  monitor,  the 
IBM-PC  places  certain  restrictions  on  the  use  of  color.  In 
the  text  mode,  any  one  of  sixteen  foreground  colors  and 
eight  background  colors  may  be  selected;  in  the  medium- 
resolution  graphics  mode,  the  user  may  choose  colors  from 
one  of  four  color  palettes,  each  containing  three  colors 
plus  the  background  color;  in  high-resolution  graphics  mode, 
the  background  color  is  always  black,  but  the  user  may 
choose  one  of  sixteen  foreground  colors. 

In  order  to  be  readily  interpreted  by  the  user,  numerical 
data  must  be  presented  in  a  clear  and  intelligible  format. 
Merely  filling  the  screen  with  numbers  is  not  necessarily 
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the  most  efficient  way  to  accomplish  this.  One  way  to 
improve  readability  is  to  group  related  data  together; 
different  groups  may  be  separated  by  collimating  the  data. 
Finally,  all  data  should  be  labeled  and  units  stated  where 
appropriate . 

Similar  requirements  exist  for  graphical  data,  but  now 
there  are  two  modes  to  worry  about.  When  multiple  curves 
are  drawn  on  the  same  graph,  the  medium- resolution  mode 
might  be  more  appropriate  since  it  allows  multiple  colors. 
However,  if  the  curves  are  highly  detailed,  the  high- 
resolution  mode  might  better  emphasize  these  details.  So 
there  is  a  trade-off  between  color  differentiation  among 
data  and  resolution.  Note  also  that  it  takes  much  longer  to 
obtain  hardcopy  of  high-resolution  displays  than  of  medium- 
resolution  displays. 

Software  Validation  (5:36-119) 

It  is  absolutely  essential  that  all  software  be  tested 
not  only  extensively,  but  intensively  as  well.  In  other 
words,  while  the  programmer  must  obviously  test  all  routines 
in  the  program  under  scrutiny,  he/she  should  also  "pick 
apart"  each  routine  to  examine  how  it  functions  under  all 
possible  circumstances.  Despite  all  efforts  at  objectivity, 
this  testing  philosophy  nevertheless  tends  to  be  counter¬ 
intuitive  to  the  programmer  testing  his  own  code.  For  this 
reason,  a  rigorous  testing  protocol  should  be  developed  and 
adhered  to.  This  section  will  address  the  basic  principles 
of  test  case  formulation  and  program  testing  procedures. 


Test  Case  Formulation.  There  are  five  general  test  case 
formulation  methodologies  of  interest  to  the  programmer. 

The  first  is  a  so-called  "white-box"  methodology;  the  rest 
are  "black-box"  methodologies.  The  term  "White-box"  implies 
that  the  tester  is  concerned  with  the  internal  structure  of 
the  program.  White-box  testing  is  concerned  with  the  degree 
to  which  the  established  test  cases  span  the  logic  of  the 
program.  Therefore/  the  test  data  itself  is  generated  by  a 
detailed  examination  of  program  logic.  The  ultimate  white- 
box  test  is  the  execution  of  every  control  path  in  the 
program.  However,  this  is  generally  not  feasible  for 
programs  with  many  levels  of  loops.  One  type  of  testing, 
logic-coverage,  represents  an  attempt  to  establish  alternate 
criteria  for  program  correctness. 

There  are  several  approaches  to  logic-coverage  testing. 
The  simplest  is  known  as  statement  coverage.  This  criterion 
requires  that  every  statement  in  the  program  be  executed  at 
least  once.  A  more  stringent  test  criterion  is  that  of 
decision  coverage,  which  requires  that  test  data  be 
generated  such  that  each  decision  has  at  least  one  true  and 
one  false  outcome.  However,  decision  coverage  is  inadequate 
if  multiple  conditions  are  involved  in  a  decision. 

Therefore  a  more  complete  criterion  is  that  of  condition 
coverage.  This  criterion  requires  that  each  condition  in  a 
decision  takes  on  all  possible  values  at  least  once.  The 
above  criteria  can  also  be  used  together  in  a  hybrid 
approach . 


vy 


The  remaining  discussion  of  test  case  formulation 
concerns  "black-box"  methodologies.  The  term  "black-box" 
implies  that  the  tester  is  not  concerned  with  the  internal 
structure  or  behavior  of  the  program,  but  rather  with  the 
conditions  under  which  the  code  fails  to  meet  its  initial 
requirements.  Therefore,  the  test  data  is  derived  solely 
from  performance  specifications.  The  ultimate  black-box 
test  is  to  perform  exhaustive  input  testing  over  all 
possible  inputs.  For  more  complicated  programs  this  number 
becomes  astronomical.  However,  this  seemingly  infinite  set 
of  test  cases  may  be  reduced  to  a  finite  and  quite 
manageable  set  if  certain  assumptions  are  made  regarding  the 
program' 8  response  to  carefully  chosen  data.  This  is  the 
goal  of  black-box  testing.  Black-box  methodologies  include 
equivalence  partitioning,  boundary- value  analysis,  and  error 
guessing . 

In  selecting  a  small  subset  of  inputs  for  testing,  it  is 
obviously  important  to  choose  one  which  has  the  highest 
probability  of  finding  the  most  errors.  A  well-selected 
test  case  should  cover  a  large  subset  of  other  possible  test 
cases,  and  thereby  reduce  the  number  of  other  test  cases 
needed.  This  may  be  accomplished  by  partitioning  the  input 
into  a  finite  number  of  equivalence  classes.  An  equivalence 
class  is  a  set  of  input  data  chosen  such  that  a  test 
performed  on  a  representative  value  of  each  class  is 
equivalent  to  a  test  performed  on  any  other  value  in  that 
class.  All  possible  input  conditions  should  be  partitioned 
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into  valid  and  invalid  equivalence  classes.  Test  cases 
should  be  written  covering  as  many  valid  equivalence  classes 
as  possible  for  the  sake  of  thoroughness.  Individual  test 
cases,  each  covering  a  single  invalid  equivalence  class, 
should  also  be  written.  This  is  necessary  because  some 
erroneous  input  checks  may  mask  others. 

In  equivalence  partitioning,  a  random  member  of  a  class 
is  the  subject  of  the  test.  However,  this  methodology 
ignores  the  case  of  a  value  directly  on  the  edge  of  a  valid 
equivalence  class.  The  methodology  of  specifically 
examining  the  edges  of  valid  equivalence  classes  is  known  as 
boundary-value  analysis.  Test  cases  are  derived  by 
considering  both  input  and  output  equivalence  classes.  One 
should  check  the  range  limits  of  test  variables  as  well  as 
just  beyond  these  limits. 

The  last  black-box  testing  methodology  to  be  discussed  is 
error  guessing.  Error  guessing  is  an  extremely  informal 
technique  based  on  programmer  intuition  and  experience.  One 
approach  is  to  simply  list  possible  error  situations,  and 
then  write  test  cases  based  upon  that  list.  Another 
approach  is  to  identify  test  cases  based  on  prograirmer 
assumptions  made  regarding  the  initial  specifications.  In 
any  case,  the  development  of  test  cases  generally  revolves 
around  the  programmer's  hunches  concerning  what  types  and 
combinations  of  input  will  generate  an  error. 

It  should  be  noted  that  none  of  the  above  methods  is 
complete  in  itself.  A  prudent  test  strategy  would  certainly 


invoke  several,  if  not  all,  of  these  methods.  A  logical 
order  of  testing  would  be  boundary-value  analysis  first, 
followed  by  equivalence  partitioning,  error  guessing,  and 
logic-coverage . 

Program  Testing.  A  large  program  is  generally  not  tested 
all  at  once,  but  rather  functional  modules  are  tested 
individually.  By  narrowing  the  scope  of  testing  to  a  more 
manageable  level,  it  becomes  possible  to  achieve  a  certain 
element  of  parallelism  by  testing  several  modules 
simultaneously.  Debugging  is  also  facilitated,  since  it  is 
easier  to  discover  and  trace  errors. 

There  are  two  considerations  in  module  testings  module 
function  and  module  interconnection.  Corresponding  to  this 
are  two  major  philosophies  of  program  testing:  non- 
incremental  testing  and  incremental  testing.  In  non- 
incremental  testing,  all  modules  are  first  tested 
individually  and  then  combined  to  form  the  program.  In 
order  to  do  this,  special  driver  modules  must  be  written  to 
transmit  test  cases  to  the  module  under  test  and  display 
results.  In  addition,  one  or  more  stub  modules  must  be 
written  to  receive  control  when  a  module  under  test  invokes 
another.  The  stubs  should  return  a  meaningful  result  to  the 
calling  module  as  well  as  returning  control. 

In  incremental  testing,  the  module  under  test  is  combined 
with  previously  checked  and  validated  modules.  As  modules 
are  progressively  linked  together,  the  individual  module 
functions  as  well  as  the  soundness  of  their  interconnections 


are  tested,  until  the  entire  program  has  been  pieced 
together.  Incremental  testing  has  several  advantages  over 
non-incremental  testing.  First,  there  are  fewer  drivers  and 
stubs  to  write.  Also,  errors  related  to  interface  problems 
will  be  detected  earlier,  since  combinations  of  modules  are 
tested  together  early.  Finally,  incremental  testing  results 
in  a  more  thorough  test,  since  executing  a  module  as  a 
result  of  an  invocation  by  another  module  might  trigger  a 
previously  undetected  error  condition. 

There  are  two  basic  strategies  for  performing  incremental 
testing:  top-down  testing  and  bottom-up  testing.  They  each 
have  their  relative  advantages  and  disadvantages,  and  the 
choice  of  which  method  to  use  is  highly  dependent  upon  the 
characteristics  of  the  program  being  tested. 

In  top-down  testing,  one  starts  at  the  top  of  the  program 
and  tests  downward.  No  drivers  are  needed  since  we  are 
beginning  execution  at  the  top  of  the  program,  but  stubs 
must  be  prepared  for  each  non-terminal  module  at  the  current 
hierarchical  level  of  test.  A  terminal  module  is  simply  one 
which  does  not  call  any  other  modules.  The  test  proceeds 
downwards  until  the  program  has  been  formed.  The  only 
criterion  for  determining  the  next  module  to  test  is  that  at 
least  one  of  that  module's  calling  modules  must  have  been 
previously  tested.  Top-down  testing  provides  an  initial 
program  skeleton  that  allows  early  program  demonstrations 
and  proves  the  overall  design. 

In  bottom-up  testing,  no  stubs  are  needed  since  testing 
begins  at  the  terminal  modules,  but  drivers  must  be  written 


for  the  nodules  under  test.  The  only  criterion  for  choosing 
the  next  module  to  test  is  that  all  of  its  subordinate 
modules  must  have  been  previously  tested.  Testing  proceeds 
upwards  until  the  program  has  been  formed.  Test  conditions 
are  easier  to  create  in  bottom-up  testing,  and  the 
observation  of  test  results  is  easier.  However,  the  program 
does  not  exist  as  an  entity  or  even  as  a  functional  subset 
until  the  last  module  has  been  added. 

After  module  testing  has  been  completed,  some  type  of 
higher-order  testing  is  appropriate.  Software  development 
is  largely  a  communication  process,  where  information 
regarding  the  eventual  program  is  translated  into  various 
forms.  For  example,  system  requirements  are  translated  into 
program  objectives,  and  these  are  translated  into  external 
specifications,  and  so  forth.  Breakdowns  in  this 
communication  process  may  readily  lead  to  errors. 

The  first  higher-order  test  usually  performed  is  function 
testing.  This  activity  is  designed  to  find  discrepancies 
between  a  program  and  its  external  specification.  Test 
cases  are  derived  from  the  external  specifications,  and 
testing  is  accomplished  using  the  previously  mentioned 
techniques . 

Another  class  of  higher-order  tests  is  system  testing. 

Its  purpose  is  to  compare  the  program  to  its  original 
objectives.  Here,  specifications  cannot  be  used  to  derive 
test  cases,  since  this  would  be  defeating  the  purpose  of  the 
test.  Instead,  test  cases  are  designed  by  analyzing  general 


program  objectives.  Some  categories  of  system  testing 
include  facility  testing  (to  see  if  every  facility  mentioned 
in  the  objectives  is  indeed  implemented),  usability  testing 
(to  determine  how  well  the  user  interacts  with  the  system), 
and  compatibility  testing  (to  see  how  easily  the  software 
can  be  installed  on  other  systems). 

This  section  has  addressed  the  issue  of  system 
requirements.  The  problem  of  stating  complete  requirements 
is  threefold,  since  one  must  address  functional  descriptions 
of  the  system,  the  user  interface  with  the  system,  as  well 
as  a  methodology  for  validating  the  end  product.  The 
requirements  stated  in  this  chapter  are  considered  to  be 
only  a  minimum  set  of  criteria. 
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IV.  SYSTEM  DESIGN 


The  purpose  of  this  section  is  to  present  the  overall 
system  design  from  a  global  perspective.  System  design  may 
be  thought  of  in  terms  of  two  distinct  categories :  system 
architecture  and  design  structure.  System  architecture 
refers  to  the  breakdown  of  the  system  into  its  major 
subsystems,  illustrating  the  scope  and  relative  hierarchy  of 
each  subsystem,  as  well  as  the  relationships  between 
subsystems.  System  design  concerns  the  functional  mapping 
into  modules,  and  the  interaction  between  modules. 

Representation  Techniques  (8:44-62) 

In  this  section,  certain  graphical  techniques  are  used  to 
represent  system  architecture  and  design  structure.  The 
following  is  a  brief  explanation  of  those  techniques . 

System  Architecture.  System  architecture  may  be 
represented  by  a  modified  Leighton  diagram.  This  type  of 
diagram  presents  the  issues  of  scope,  hierarchy,  and 
external  interfaces  without  the  complexity  associated  with 
excessive  detail.  Its  notation  consists  simply  of  labels 
and  line  segments.  The  labels  represent  logical  modules, 
and  the  line  segments  represent  logical  connection.  The 
horizontal  dimension  portrays  overall  hierarchy,  while  the 
vertical  dimension  portrays  scope  of  control. 

Design  Structure.  Design  structure  is  represented  via 
structure  charts.  The  structure  chart  illustrates  the 
overall  module  hierarchy,  as  well  as  the  interaction  between 


modules.  Its  notation  consists  of  rectangles  representing 
the  modules,  vectors  representing  control  relationships,  and 
arrows  with  circular  tails  representing  transfer  of 
data/control  (solid  circles  represent  boolean  flags). 

Synopsis  of  Overall  System  Architecture 

The  overall  architecture  of  the  system  is  concisely 
illustrated  in  Figure  4-1.  The  global  module  MAIN 
corresponds  to  the  user's  main  menu.  Via  the  appropriate 
menu  input,  control  then  passes  to  one  of  five  subsystems: 
the  satellite  parameters  group,  the  antenna  gain  patterns 
group,  the  link  analysis  group,  the  database  support  group, 
and  the  exit  routine.  Except  for  the  exit  routine,  each 
subsystem  consists  of  a  number  of  functional  modules,  and 
hence  are  referred  to  as  "groups."  Each  subsystem  will  now 
be  discussed  individually  in  detail. 

Satellite  Parameters  Group.  The  function  of  this  group 
is  to  perform  several  commonly  used  calculations  concerning 
a  satellite's  orbit  and  position.  Five  modules  within  this 
group  perform  the  following  tasks: 

1.  Compute  Earth  escape  velocity  at  an  altitude. 

2.  Compute  the  steady-state  temperature  of  a  satellite  in 
Earth  orbit. 

3.  Compute  various  parameters  associated  with  a  specified 
circular  orbit. 

4.  Compute  various  parameters  associated  with  a  specified 
elliptical  orbit. 
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Figure  4-1.  Overall  System  Architecture 
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5.  Compute  azimuth  and  elevation  angles  to  a  satellite 

from  two  different  ground  sites. 

( 

Only  the  last  module  has  any  external  connection  with  the  1 

rest  of  the  system.  The  computed  elevation  angles  are 
stored  in  memory,  and  may  be  used  as  input  during  link 

t 

analysis  calculations.  The  remainder  of  the  satellite  j 

parameter  data  is  discarded  upon  return  to  the  main  menu. 

Antenna  Gain  Patterns  Group.  The  function  of  this  group 
is  to  plot  the  gain  patterns  of  various  types  of  antennas . 

Four  modules  handle  the  pattern  plotting  for  horns,  helixes, 
parabolic  reflectors,  and  phased  arrays  of  dipoles.  The 
first  three  modules  also  display  pertinent  design 
information  along  with  the  gain  pattern  for  more  meaningful 
hardcopy.  The  last  module  does  not  display  this 
information,  but  does  display  each  of  the  six  standard  polar 
antenna  plots  separately. 

Link  Analysis  Group.  This  functional  group  constitutes 
the  heart  of  the  system.  Although  there  are  a  large  number 
of  individual  modules  within  this  group,  the  user  perceives 
only  three  distinct  processes:  link  analysis,  display  link 
budget,  and  recompute  link  budget.  The  entire  system  is,  of 
course,  organized  around  the  link  analysis  process.  In  this 
process,  the  user  is  guided  step  by  step  through  a  complex 
series  of  link  calculations,  with  intermediate  results  being 
displayed  along  the  way.  The  engineer  is  provided  with  all 
of  the  data  necessary  to  design  and  analyze  a  satellite 


link 


The  display  link  budget  process  displays  all  the 
pertinent  data  computed  above,  and  organizes  it  in  an  easily 
readable  and  understandable  format.  The  display  has 
separate  columns  for  uplink  and  downlink  data,  and  is  color- 
coded  according  to  function  and  location  in  the  link. 

The  recompute  link  budget  process  is  included  to  allow 
the  design  engineer  the  opportunity  to  alter  various  key 
parameters,  in  order  to  determine  the  overall  effect  on  the 
link.  Given  a  new  value  of  a  parameter,  the  entire  set  of 
link  calculations  is  recomputed  with  the  new  parameter 
without  the  user  having  to  manually  re-enter  all  of  the 
previous  input  data.  This  feature  becomes  especially 
valuable  when  one  is  engaged  in  the  process  of  making  design 
tradeoffs,  for  example,  allowing  a  smaller  receiving  antenna 
while  increasing  the  receiver  noise  figure. 

Database  Support  Group.  The  purpose  of  this  group  is  to 
perform  all  of  the  housekeeping  and  disk  management  tasks 
associated  with  the  link  database.  Thus,  the  user  is  spared 
the  tedium  of  keeping  track  of  what  data  is  in  what  file, 
what  files  are  open,  and  on  which  disk;  the  system  also 
prevents  the  user  from  inadvertently  erasing  or  not  saving 
an  important  file.  The  user  interfaces  with  this  group  via 
five  options  in  the  main  menu.  These  options  are  concerned 
with  storing  newly  computed  link  files,  recalling  old  link 
files,  deleting  unwanted  link  files,  listing  the  existing 
files  in  the  link  file  directory,  and  changing  the  default 
data  drive.  These  options  are  so  complete  in  their  function 
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and  behavior  that  they  strongly  resemble  the  corresponding 
DOS  commands . 

Exit  Routine.  This  routine  is  not  associated  with  any 
distinct  module,  and  is  actually  embedded  within  the  master 
program.  All  it  does  is  allow  the  user  a  clean  way  to 
return  to  DOS.  If  a  recently  computed/recomputed  link  file 
has  not  yet  been  saved,  this  routine  gives  the  user  the 
opportunity  to  do  so  before  exiting. 

Synopsis  of  Overall  Design  Structure 

At  this  point  the  system  design  structure  is  presented  in 
detail.  The  design  structure  may  be  considered  one  level  of 
detail  down  from  the  architecture.  The  emphasis  here  is  on 
module  functionality  and  module  interaction. 

Explanation  of  Module  Functionality.  For  future 
reference,  the  functional  modules  on  the  system  disk  are 
listed  alphabetically  below  and  described  according  to  basic 
function . 

AIR  -  Computes  attenuation  due  to  atmospheric  absorption. 
ALOG  -  Converts  dB  to  numbers . 

ANTENNA 1  -  Computes  gains  and  beamwidths  for  earth 
station  antennas. 

ANTENNA2  -  Same  as  above,  only  for  satellite  antennas. 
ARRAY  -  Plots  gain  pattern  for  phased  array  of  dipoles. 
AZEL  -  Computes  azimuth  and  elevation  angle  to  a 
satellite  from  two  different  ground  stations. 

BESSEL  -  Computes  first-order  Bessel  functions  of  the 
first  kind  (  J^(x)  ). 


BITRATE  -  High-level  module  which  computes  available 
E^/N^  required  Ej^/N^  and  max  data  rate. 

BLACKBOD  -  Computes  steady-state  temperature  of  a 
satellite  in  Earth  orbit. 

CAPACITY  -  Computes  Shannon  rate  and  Nyquist  rate,  and 
determines  whether  the  link  is  channel-limited. 

CIRCULAR  -  Computes  velocity,  orbital  period,  coverage 
area,  and  duration  of  visibility  for  a  satellite  in 
circular  Earth  orbit. 

CUSTOM  -  Contains  custom  plotting  and  drawing  routines. 
DISH  -  Plots  gain  pattern  for  parabolic  dish  antennas. 
ELLIPSE  -  Computes  velocity  at  apogee  and  perigee, 
average  velocity,  eccentricity,  orbital  period,  and 
coverage  area  for  a  satellite  in  elliptical  Earth  orbit. 
ESCAPE  -  Computes  escape  velocity  at  altitude. 

EXIST  -  Function  that  determines  whether  a  given  file 
exists  on  disk. 

FUNCTION  -  Low-level  module  used  by  many  other  modules  to 
perform  actual  calculations  (e.g.  gain,  beamwidth,  EIRP, 
system  temperature,  figure  of  merit,  received  useful 
power,  CNR,  available  Ej^/N^,  required  E^/N^,  BER,  and  max 
data  rate) . 

GETALT  -  Obtains  satellite  altitude  from  the  user. 

HELIX  -  Plots  gain  pattern  for  helical  antennas. 

HORN  -  Plots  gain  pattern  for  horn  antennas . 

INTEGER  -  Polls  the  user  for  an  integer  input. 

LOG10  -  Computes  log(x). 


LOG 2  -  Computes  log2(x). 

MASTER  -  Main  program. 

NOISE  -  Makes  beeping  sound  (used  to  signal  error 
conditions) . 

POWER  -  High-level  module  which  computes  received  useful 
power  and  CNR. 

PWR  -  Function  to  perform  exponentiation. 

Q  -  Computes  Marcum's  Q-function. 

QINVERSE  -  Computes  inverse  Q-function. 

RAIN  -  Computes  attenuation  due  to  rainfall. 

REAL  -  Polls  the  user  for  a  real  input. 

REC  -  High-level  module  which  computes  system  temperature 
and  figure  of  merit. 

SAVEIT  -  Routine  to  save  a  link  file  on  disk. 

TRIG  -  Package  of  functions  to  compute  arccos(x), 
arcsin(x),  and  tan(x). 

WAIT  -  Suspends  execution  until  a  key  is  pressed. 

XMIT  -  High-level  module  which  computes  EIRP. 

Module  Interactions.  The  interactions  among  the  above 
modules  can  best  be  expressed  by  structure  charts.  For  the 
sake  of  brevity,  only  those  structure  charts  illustrating 
significant  module  inter-relationships  are  discussed.  For 
example,  the  case  of  one  module  calling  another  module  to 
compute  a  simple  trigonometric  function  is  not  considered 
significant.  All  graphics  routines  call  the  module  CUSTOM; 
by  replacing  this  module  with  an  alternate  graphics  driver, 
the  system  should  be  capable  of  plotting  on  any  MS-DOS 


machine.  In  addition,  all  modules  requiring  numeric  inputs 
employ  INTEGER  and/or  REAL  to  screen  user  input.  Figure  4-2 
presents  the  most  important  structure  charts. 

From  this  figure  it  may  be  seen  that  the  modules 
CIRCULAR,  BLACKBOD,  and  AZEL,  as  well  as  the  main  program, 
all  employ  the  module  GETALT  to  obtain  a  value  for  satellite 
altitude  from  the  user.  The  main  program  interfaces  with 
AZEL  for  the  purpose  of  obtaining  values  for  elevation 
angle.  If  the  flag  HAVEALPHA  is  true,  then  the  user  may 
choose  to  continue  using  the  elevation  angles  previously 
computed  in  AZEL.  The  modules  HORN,  HELIX,  and  DISH  all  use 
the  module  FUNCTION  to  obtain  values  for  maximum  gain.  DISH 
also  uses  the  module  BESSEL  to  actually  calculate  values  for 
J^x).  The  modules  ANTENNA  1  and  ANTENNA2  both  employ  the 
module  FUNCTION  to  obtain  values  for  gain  and  beamwidth. 

The  module  REC  uses  FUNCTION  to  compute  system  equivalent 
temperature  and  figure  of  merit;  XMIT  uses  it  to  compute 
EIRP.  The  module  POWER  calls  AIR  and  RAIN  to  obtain  values 
for  attenuation  due  to  the  atmosphere  and  to  rainfall, 
respectively,  and  FUNCTION  to  compute  received  power.  The 
module  BITRATE  uses  FUNCTION  to  compute  available  E^/Nq, 
required  E^/N^  and  max  data  rate.  It  also  calls  CAPACITY 
to  compute  the  Shannon  rate  and  the  Nyquist  rate;  if  one  of 
these  rates  is  lower  than  the  computed  max  data  rate,  then 
an  appropriate  flag  is  set  and  BITRATE  makes  the 
corresponding  substitution  for  the  max  data  rate. 

This  section  has  given  the  user  an  overview  of  the 
overall  system  design.  There  are  two  fundamental  facets  to 


system  design t  system  architecture  and  design  structure. 
System  architecture  provides  the  user  with  insight  regarding 
the  decomposition  of  the  overall  system  into  its  constituent 
subsystems;  design  structure  examines  the  functions  and 
inter-relationships  of  the  modules  which  serve  as  building 
blocks  for  the  system.  A  working  knowledge  of  the  overall 
system  design  gives  the  user  a  feel  for  how  the  system 
operates,  without  the  frustration  of  having  to  dredge 
through  complicated  equations  or  piles  of  source  code. 


V.  DETAILED  DESIGN 


In  the  previous  section  the  system's  functional  modules 
were  considered  black  boxes,  and  no  attempt  was  made  to 
explain  their  inner  workings.  It  is  the  purpose  of  this 
section  to  examine  in  detail  the  structure  and  behavior  of 
these  modules.  The  issues  of  control  flow  and  database 
management  are  also  covered. 

Representation  Techniques 

In  this  section,  certain  graphical  techniques  are  used  to 
represent  control  flow  and  database  structure.  The 
following  is  a  brief  explanation  of  those  techniques. 

Control  Graphs  (8*96-98).  Control  graphs  are  used  to 
show  the  mechanisms  by  which  a  program  may  transition  from 

one  logical  state  to  another.  It  is  implicit  that  a  program 
can  only  exist  in  one  state  at  any  given  time.  The  notation 
consists  of  node  symbols  representing  a  given  state,  solid 
lines  representing  a  control  path,  and  dashed  lines 
representing  input  from  some  external  condition.  Control 
graphs  are  borrowed  from  automata  theory,  and  are  used  to 
establish  that  all  nodes  in  a  system  can  be  reached,  and 
that  once  reached  they  can  be  exited.  The  nodes  themselves 
may  be  only  conceptual  entities,  and  the  resulting  control 
graph  becomes  an  abstract  program  model  rather  than  a 
physical  model. 

Database  Structure  (81 77-81).  The  purpose  of  database 
structure  representation  is  twofold*  to  show  flow  of  data 
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within  the  system,  and  to  illustrate  the  basic  file 
structure.  There  are  many  types  of  data  structure 
representation  techniques,  each  with  its  own  particular 
merits.  Because  of  the  relatively  simple  database 
constructs  in  this  system,  only  a  rudimentary  technique  is 
needed.  The  data  structure  diagram  used  to  describe  this 
system  consists  simply  of  boxes  representing  disk  files,  and 
vectors  representing  directional  access.  If  a  key  is 
required  to  access  a  particular  file,  then  it  is  shown  in  a 
box  adjoining  the  corresponding  file. 


Key  Detailed  Design  Issues 

The  previous  chapter  on  system  design  examined  the  system 
from  a  logical  point  of  view,  and  the  items  discussed 
therein  do  not  necessarily  correspond  to  physical, 
implement able  sections  of  the  program.  In  addition,  all 
functional  modules  were  treated  as  black  boxes.  It  is  the 
purpose  of  this  section  to  examine  the  system  from  a 
physical  perspective,  and  to  describe  the  system  in  terms  of 
software  behavior.  Three  major  topics  are  discussed:  flow 
of  control,  database  management,  and  software  design. 

Flow  of  Control .  Figure  5-1  illustrates  the  system  flow 
of  control.  There  are  13  states  in  the  system,  which 
correspond  to  the  main  menu,  the  11  main  menu  selections, 
and  DOS.  Although  DOS  is  technically  not  a  program  state, 
it  is  considered  a  system  state .  It  may  be  seen  from  this 
figure  that  the  initial  program  state  is  MAIN,  which 
corresponds  to  the  main  menu  from  the  user's  perspective. 
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The  states  Sat,  Gain,  Disp,  and  Drive  correspond 
respectively  to  four  selections  from  the  main  menus 
satellite  parameters,  antenna  gain  patterns,  display  link 
budget,  and  change  default  drive.  These  states  are  all 
characterized  by  the  fact  that  they  are  entered  only  from 
MAIN,  and  once  exited,  control  again  returns  to  MAIN. 

All  of  the  remaining  states  also  possess  this  control 
path,  but  are  also  affected  by  other  states  or  external 
conditions  (such  as  program  flags).  For  example,  when  the 
the  state  LINK  is  entered,  it  provides  an  external  input  to 
RECOMP  indicating  that  link  analysis  has  just  been 
performed.  Likewise,  when  the  state  RCL  is  entered,  it 
provides  an  external  input  to  RECQMP  indicating  that  a  link 
file  has  been  recalled  from  disk.  The  purpose  of  this 
arrangement  is  to  disallow  recomputing  a  link  budget  unless 
a  link  budget  has  either  been  just  computed  or  just  recalled 
from  disk.  This  prevents  the  fatal  errors  which  would  occur 
if  the  recompute  module  were  directed  to  recompute  a  link 
file  which  consisted  only  of  garbage  data.  Multiple 
recomputing  of  the  same  file  is  allowed.  Note  that  the 
external  condition  EXIST  (not  a  state)  will  not  allow 
control  to  pass  to  RCL  unless  the  file  designated  for  recall 
does  in  fact  exist  on  the  default  drive.  Similarly,  control 
cannot  pass  to  DEL  unless  the  file  designated  for  deletion 
exists  on  the  default  drive. 

The  state  LIST  cannot  be  entered  unless  the  link  file 
directory  is  non-empty.  This  prevents  the  display  of  a 
blank  screen  in  the  case  of  an  empty  directory. 

50 


The  remaining  states  are  closely  intertwined.  The  state 
STO  has  the  usual  control  path  to  and  from  MAIN.  The  state 
EXIT  can  always  be  entered  from  MAIN;  however,  control  does 
not  return  to  DOS  unless  either  there  is  no  unsaved  link 
data,  or  there  is  unsaved  link  data  but  the  user  has  given 
permission  to  discard  it.  If  neither  of  these  two 
conditions  exists,  then  control  passes  to  STO  so  that  the 
data  can  be  stored  on  disk,  and  then  to  DOS.  This 
arrangement  eliminates  the  possibility  of  exiting  the  system 
without  saving  important  data. 

Database  Management.  The  system  database  is  relatively 
simple,  and  consists  of  three  basic  filetypest  link 
directory,  buffer,  and  link  file.  In  the  current  system 
there  is  one  link  directory,  one  buffer,  and  numerous  link 
files.  The  link  directory  is  designated  LINK. DIR,  and 
contains  the  names  of  the  link  files  currently  stored  on 
disk.  The  buffer  file  TEMP. STO  is  created  only  during  the 
file  deletion  process,  and  is  used  to  temporarily  hold 
directory  entries.  Each  link  file  contains  all  the  data  for 
one  link  budget,  and  may  be  given  any  name  from  one  to  eight 
characters  long.  All  of  the  above  files  reside  on  the 
default  drive  under  the  directory  path  drive xdata\ 
[filename].  Figure  5-2  graphically  shows  the  system 
database  structure. 

When  a  link  file  is  to  be  stored  under  a  given  filename, 
the  system  first  checks  to  see  if  there  is  already  a  link 
file  on  disk  with  the  same  filename.  If  so,  a  warning 


Figure  5-2.  System  Database  Structure 
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notice  is  issued  and  the  operation  is  disallowed.  Otherwise 
the  file  is  written  to  disk  with  the  extension  . LNK,  and  an 
appropriate  entry  is  made  in  the  directory. 

When  a  link  file  is  to  be  recalled  from  disk,  the  system 
first  checks  to  see  if  the  requested  file  exists  on  disk. 

If  not,  a  warning  notice  is  issued  and  the  operation  is 
again  disallowed  in  order  to  prevent  a  disk  I/O  error.  If 
it  does  exist,  then  it  is  opened  and  read.  The  user  is 
encouraged  to  check  the  directory  for  the  existence  of  a 
file  before  trying  to  recall  it. 

When  an  unwanted  link  file  is  to  be  deleted  from  disk, 
the  same  check  for  existence  is  performed  as  above.  If  the 
file  does  indeed  exist,  the  user  is  asked  to  confirm  the 
decision  to  delete.  If  the  user  responds  in  the 
affirmative,  then  the  designated  link  file  is  erased  from 
disk,  and  the  corresponding  entry  is  removed  from  the 
directory.  However,  since  the  actual  location  of  the 
deleted  filename  in  the  directory  is  unknown,  the  directory 
LINK. DIR  is  first  copied  into  the  buffer  TEMP.STO.  Next, 
LINK. DIR  is  purged  of  all  records.  Finally,  all  records  in 
TEMP.STO  with  the  exception  of  the  one  being  deleted  are 
written  back  into  LINK. DIR. 

When  the  directory  is  to  be  listed,  a  check  is  made  to 
determine  whether  there  is  anything  in  LINK. DIR  to  display. 
If  not,  an  appropriate  message  is  printed  and  no  action  is 
taken.  If  LINK. DIR  is  non-empty,  then  it  is  simply  listed 
record  by  record. 


Changing  the  default  drive  merely  changes  the  path  to  the 


link  data  from  olddrive :data\[ filename]  to  newdrive :data\ 
[filename].  This  permits  multiple  data  disks  to  be  used  at 
the  same  time. 

Software  Design 

The  purpose  of  this  section  is  not  to  derive  or  verify 
equations  used  in  the  calculations,  but  rather  to  explain  in 
a  somewhat  general  sense  the  algorithms,  engineering 
principles,  and  assumptions  involved  in  the  design  of  each 
functional  module.  Flowcharting  every  module  and  discussing 
every  line  of  source  code  does  little  to  further  the  user's 
understanding  of  the  system.  Hence,  the  software  is 
discussed  more  qualitatively  than  analytically.  For  the 
purposes  of  this  discussion,  all  software  is  grouped  into 
one  of  six  packages:  utility,  math,  graphics,  low-level 
functional,  high-level  functional,  and  main  program.  For 
further  detail,  the  reader  should  refer  to  the  software 
listings  in  Appendix  D. 

Utility  Package.  The  utility  package  is  a  set  of  small 
modules  which  perform  miscellaneous  program  tasks.  These 
modules  are  described  below: 

CUSTOM  -  Contains  the  procedures  CUSTOMPLOT  and 
CUSTOMDRAW,  which  in  turn  use  the  standard  Pascal 
procedures  PLOT  and  DRAW  to  plot  points  and  draw  lines  in 
graphics  mode.  If  the  system  is  transported  to  another 
machine  which  has  different  graphics  protocols,  CUSTOM 
can  be  rewritten  to  acconmodate  the  new  machine.  All 


graphics  routines  use  CUSTOM,  so  modifying  it  will 
automatically  modify  all  the  graphics  routines. 

EXIST  -  Function  which  checks  to  see  whether  a  specified 
file  exists  on  disk.  It  works  by  turning  on  the  $1- 
compiler  directive  (which  activates  I/O  error  checking) 
and  then  attempting  to  reset  the  file.  If  the  file  does 
not  exist,  this  operation  will  result  in  an  I/O  error. 

The  result  of  this  test  is  returned  to  the  calling 
routine  by  a  boolean  variable. 

GETALT  -  This  module  simply  prompts  the  user  for  a  value 
of  satellite  altitude.  The  default  value  "G"  is 
geosynchronous  altitude.  This  particular  routine  is  made 
into  a  module  because  it  is  called  by  many  other  modules 
(it  is  essentially  a  subroutine). 

INTEGER  -  Polls  the  user  for  integer  numeric  data.  If 
the  user  does  not  enter  a  legal  integer,  then  an  error 
flag  is  set  which  can  be  detected  by  the  calling  module. 
NOISE  -  Uses  the  Pascal  procedure  SOUND  to  make  the 
computer 's  speaker  emanate  a  high-pitched  tone.  It  is 
used  where  an  audible  warning  signal  is  desired. 

REAL  -  Similar  to  INTEGER,  only  this  module  polls  the 
user  for  real  numeric  data. 

SAVEIT  -  Saves  a  specified  file  on  disk.  In  order  to 
prevent  overwriting  an  existing  file,  the  directory  is 
first  checked  to  see  if  the  specified  filename  already 
exists.  If  so,  a  warning  is  issued,  and  control  returns 
to  the  main  menu.  Otherwise,  the  current  link  file  is 
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written  to  the  default  drive  under  the  specified 
filename,  with  the  extension  M.LNK.M  This  module  is 
accessed  in  two  places:  once  in  the  main  "Store  Link 
File"  routine,  and  once  at  the  end  of  the  program  in  the 
case  that  the  user  has  mistakenly  chosen  to  exit  without 
storing  a  needed  link  file. 

WAIT  -  Suspends  execution  until  a  key  is  pressed. 

Control  enters  an  infinite  loop  until  the  Pascal  variable 
KEYPRESSED  is  true.  This  routine  is  used  to  temporarily 
halt  program  execution  so  that  intermediate  results  may 
be  displayed. 

Math  Package.  The  math  package  is  a  set  of  modules  which 
perform  various  basic  mathematical  functions.  For  more 
information,  refer  to  any  text  of  standard  math  tables. 

These  modules  are  described  below: 

ALOG  -  Converts  dB  to  numbers  by  raising  10  to  the  power 
of  dB/ 10 . 

BESSEL  -  Computes  first  order  Bessel  functions  of  the 
first  kind,  This  is  accomplished  by  evaluating 

the  first  200  terms  of  the  series  expansion  for  J1(x)  (as 
given  in  any  handbook  of  standard  math  tables,  such  as 
the  CRC) .  Because  of  the  extreme  magnitude  of  some  of 
the  terms  in  the  series,  an  overflow  condition  will 
result  if  they  are  computed  directly.  To  circumvent 
this,  each  term  is  broken  up  into  smaller  pieces  which  by 
themselves  cannot  exceed  the  capacity  of  the  machine; 
these  pieces  are  then  combined  to  complete  the  series. 


If  more  than  two  consecutive  terms  do  not  contribute 
anything  to  the  series  (i.e.  they  go  to  zero  in  the 
limit)  then  the  series  is  truncated  before  the  full  200 
terms. 

LOG10  -  Function  to  compute  base  ten  logarithms  via 
natural  logarithms. 

L0G2  -  Function  to  compute  base  two  logarithms  via 
natural  logarithms. 

PWR  -  Performs  real  exponentiation  by  using  laws  of 
logarithms.  The  cases  of  zero  and  negative  arguments  are 
taken  into  account . 

Q  -  Computes  Marcum's  Q-function  of  the  argument.  If  the 
argument  is  greater  than  13  then  the  value  zero  is 
returned  to  avoid  an  overflow  condition.  If  the  argument 
is  between  4.24  and  13,  then  the  exponential 
approximation  of  the  Q-function  is  used.  Otherwise,  the 
Q-function  is  computed  by  first  calculating  the  error 
function  of  the  argument,  and  then  deriving  the  Q- 
function  from  that.  This  is  done  to  facilitate  the 
Simpson ' 8  rule  numerical  integration  which  computes  the 
error  function  (16:656). 

QINVERSE  -  Computes  the  inverse  Q-function  of  the 
argument.  The  equation  for  the  inverse  Q-function  is 
derived  from  the  same  exponential  Q-function 
approximation  used  above.  However,  since  this  is  a 
nonlinear  equation,  and  it  must  be  solved  numerically. 
Initially,  the  method  of  halving  the  interval  is  used  to 


isolate  the  general  region  of  the  root .  This  method  is 
predicated  upon  the  fact  that  if  the  sign  of  a  continuous 
function  changes  between  any  two  points,  then  a  root  must 
exist  within  that  interval;  the  region  of  sign  change  may 
then  be  isolated  as  closely  as  possible.  Once  this 
initial  estimate  of  the  root  is  made,  the  result  is 
passed  to  a  routine  which  uses  Newton ' s  method  to  more 
exactly  determine  the  root.  Newton's  method  extrapolates 
along  a  tangent  from  a  point  near  the  root  to  its 
intersection  with  the  x-axis.  This  process  is  iterated 
until  the  tangent  lines  are  sufficiently  close  together; 
the  final  point  of  x-axis  intersection  is  the  root. 

TRIG  -  Computes  some  simple  trigonometric  functions  not 
included  in  the  standard  Pascal  library.  The  tangent 
function  is  computed  by  dividing  the  sine  function  by  the 
cosine  function.  The  inverse  sine  and  inverse  cosine 
functions  are  computed  by  evaluating  simple  formulas 
involving  the  inverse  tangent  function. 

Graphics  Package.  The  graphics  package  is  a  set  of  four 
modules  which  plot  antenna  gain  patterns.  This  section  does 
not  delve  deeply  into  electromagnetic  theory,  but  rather 
explains  the  fundamentals  of  how  these  plotting  routines 
work.  All  of  the  plotting  routines  employ  the  utility 
routine  CUSTOM  to  perform  actual  plotting.  These  modules 
are  described  below: 

HORN  -  Plots  the  gain  pattern  for  a  square  horn  antenna. 
The  gain  function  for  a  horn  antenna  is  given  by 
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G=[(sin  x)/x]  ,  where  x=pi*sin(theta) .  This  module 
allows  theta  to  sweep  from  -pi/2  to  pi/2  in  increments  of 
.05,  and  calculates  corresponding  values  for  G.  Since 
these  are  in  terms  of  polar  coordinates,  they  are 
converted  to  rectangular  coordinates  for  the  purposes  of 
plotting  on  the  screen.  Appropriate  offsets  and  scaling 
coefficients  are  also  used  to  insure  that  the  resultant 
plot  fits  on  the  screen  as  desired. 

HELIX  -  Plots  the  gain  pattern  for  helical  antennas.  The 
gain  expression  for  these  antennas  is  not  quite  as 
simple,  and  is  derived  via  antenna  theory.  To  do  this, 
the  helix  is  modeled  mathematically  as  an  array  of 
circular  current  loops.  Two  quantities  are  derived,  a 
pattern  factor  due  to  the  circular  loops,  and  an  array 
factor  due  to  the  iteration  of  the  loops.  Then,  using 
pattern  multiplication  techniques,  the  overall  gain 
pattern  expression  is  formed.  This  expression  is  plotted 
in  a  manner  similar  to  the  above. 

DISH  -  Plot 8  the  gain  pattern  for  parabolic  reflector 
antennas.  The  gain  function  for  parabolic  reflectors  is 
given  by  G«(2J^(x)/x)  ,  where  x=pi*d*sin(theta)/L,  L 

is  wavelength,  and  d  is  dish  diameter.  Note  that  the 
gain  function  is  similar  to  a  sine-over-argument  squared 
function,  but  instead  uses  a  Bessel  function. 
Consequently,  DISH  calls  the  module  BESSEL  to  actually 
compute  these  values.  The  resultant  gain  function  values 
are  then  plotted  as  before. 


ARRAY  -  Plots  the  gain  pattern  for  a  planar  phased  array 
of  dipoles.  The  derivation  of  this  depends  heavily  upon 
array  theory.  This  module  first  defines  a  number  of 
functions  to  be  used  later  on  in  the  computations;  these 
being  the  three  rectangular  direction  cosines,  the  theta 
and  phi  components  of  the  element  factors  corresponding 
to  the  three  possible  vector  current  directions,  and  the 
pattern  factors  corresponding  to  the  seven  possible 
current  distributions.  The  actual  values  of  these 
quantities  depend  on  the  user  input.  Next,  the  main 
pattern  computation  routine  computes  the  array  factors 
for  the  x  and  y  directions,  and  multiplies  them  together 
along  with  the  pattern  factor  and  appropriate  element 
factor  to  obtain  an  overall  expression  for  the  gain 
pattern.  Finally,  a  plotting  routine  is  used  to  actually 
display  the  computed  gain  patterns.  Six  polar  plots  are 
displayed:  the  theta-pi/2  plane,  the  phi«0  plane,  and  the 
phi»pi/2  plane  (both  theta  and  phi  components  for  each). 
The  user  also  has  the  capability  to  magnify  or  shrink  the 
resultant  plots  to  any  size  desired. 

Low-Level  Functional  Package.  The  modules  in  this  set 
constitute  the  basis  of  the  system,  since  they  perform  the 
bulk  of  the  actual  link  analysis.  These  modules  are 
described  below: 

AIR  -  Computes  the  total  path  attenuation  due  to 
atmospheric  absorption.  This  attenuation  is  a  function 
of  two  variables:  elevation  angle  and  frequency.  The 


computations  in  this  module  are  based  upon  a  mathematical 
model  of  atmospheric  attenuation.  For  detailed 
information  on  this  model,  the  user  is  referred  to 
Appendix  B. 

AZEL  -  Computes  the  azimuth  and  elevation  angles  from  two 
specified  ground  sites  to  a  satellite  at  a  given 
position.  The  user  is  required  to  enter  the  latitude  and 
longitude  of  the  uplink  site,  downlink  site,  and  sub¬ 
satellite  point;  as  well  as  the  satellite's  altitude. 

The  azimuth  is  computed  in  a  straightforward  manner  using 
basic  spherical  trigonometry.  The  elevation  angle  is 
derived  via  simple  plane  trigonometry  from  the  tilt 
angle,  which  is  the  angle  between  a  line  from  the 
satellite  to  the  ground  site,  and  a  normal  from  the 
satellite  to  the  ground.  The  altitude  is  normalized  with 
respect  to  the  Earth's  radius,  in  effect  making  the  Earth 
a  unit  circle.  The  tilt  angle  follows  from  basic  plane 
and  spherical  trigonometry  (refer  to  Chapter  II). 

BLACKBOD  -  Computes  the  steady-state  temperature  of  a 
satellite  given  its  altitude,  total  area,  and  fraction  of 
total  area  exposed  to  the  sun.  The  value  returned 
only  a  crude  estimate,  and  neglects  such  factors  as 
rotation  of  the  satellite,  internal  heat  sources,  and 
surface  material  considerations.  The  module  models  the 
sun  as  a  perfect  blackbody,  and  calculates  the  total  flux 
at  the  satellite's  altitude.  The  steady-state 
temperature  is  found  by  setting  the  flux  gained  equal  to 
the  flux  lost. 


CAPACITY  -  Computes  the  theoretical  Shannon  rate  and 
Nyquist  rate  for  the  specified  channel,  and  determines 
whether  one  of  these  rates  is  limiting  system  performance 
with  the  parameters  specified.  The  Shannon  and  Nyquist 
rates  are  given  by 

^shannon  “  <38> 

Rnyqui.t  *  <39> 

where  R  is  channel  capacity  (bits/sec),  B  is  bandwidth 
(Hz),  P  is  received  power  (W),  Nq  is  the  single-sided 
noise  power  spectral  density  (W/Hz),  and  M  is  the  number 
of  signalling  levels  (16x7).  If  one  of  these  quantities 
is  less  than  the  previously  computed  bit  rate,  then  an 
appropriate  flag  is  set. 

CIRCULAR  -  Computes  the  velocity,  orbital  period, 
coverage  area,  and  duration  of  visibility  for  a  satellite 
in  circular  Earth  orbit.  Formulas  for  the  first  two 
quantities  may  be  found  in  any  orbital  mechanics  text. 

The  last  two  quantities  are  derived  via  simple 
trigonometry.  For  more  detail,  refer  to  Chapter  II. 
ELLIPSE  -  Computes  the  velocities  at  apogee  and  perigee, 
average  velocity,  velocity  at  a  specified  orbital  point, 
eccentricity,  orbital  period,  and  coverage  area  for  a 
satellite  in  elliptical  Earth  orbit.  Formulas  for  these 
quantities  follow  from  basic  trigonometry  and  analytic 
geometry,  and  may  be  found  in  any  orbital  mechanics  text. 
For  more  detail,  refer  to  Chapter  II. 


ESCAPE  -  Computes  escape  velocity  at  a  given  altitude. 

The  formula  for  escape  velocity  may  be  found  in  any 
physics  or  orbital  mechanics  text.  For  more  detail, 
refer  to  Chapter  II. 

FUNCTION  -  Contains  a  number  of  link  analysis  functions 
which  are  called  by  several  other  modules.  FUNCTION 
contains  functions  to  compute  gains  and  3-dB  beamwidths 
for  four  different  antenna  types,  effective  isotropic 
radiated  power  (EIRP),  system  equivalent  temperature,  G/T 
figure  of  merit,  received  useful  power,  carrier  power-to- 
noise  power  density  ratio,  available  and  required  bit 
energy-to-noise  power  density  ratios,  bit  error  rate,  and 
maximum  data  rate.  For  more  information  on  the 
computation  of  these  quantities,  refer  to  Chapter  II. 

RAIN  -  Computes  the  attenuation  due  to  rainfall  in  the 
channel.  Like  AIR,  this  module  relies  upon  a 
mathematical  model  of  rainfall  attenuation  as  a  function 
of  rain  rate,  frequency,  and  elevation  angle.  For 
detailed  information,  the  reader  is  referred  to  Appendix 
C. 
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High-Level  Functional  Package.  The  modules  in  this 
package  do  not  perform  link  calculations  per  se,  but  rather 
organize  user  input  and  call  the  appropriate  low-level 
modules.  The  high-level  system  modules  are  discussed  below: 
ANTENNA 1  -  Allows  the  user  to  select  the  antenna  for 
uplink  transmitting  or  downlink  receiving  at  the 
corresponding  ground  site.  The  user  may  choose  a  horn, 
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helix,  parabolic  reflector,  or  phased  array  of  dipoles. 
The  module  then  computes  the  gain  and  beamwidth  for  the 
particular  antenna  chosen.  These  quantities  are  computed 
by  calling  the  appropriate  routine  in  the  module 
FUNCTION. 

ANTENNA2  -  Allows  the  user  to  select  the  antenna  for 
uplink  receiving  or  downlink  transmitting  on  the 
satellite.  These  antennas  are  currently  restricted  to  be 
parabolic  reflectors.  However,  the  user  may  specify 
Earth  coverage  or  spot  coverage.  The  gain  and  beamwidth 
are  computed  as  before,  along  with  the  gain  and  beamwidth 
required  to  achieve  the  specified  coverage  format.  The 
user  is  alerted  if  either  the  actual  gain  or  actual 
beamwidth  is  insufficient  to  achieve  25-dB  sidelobes. 

REC  -  Prompts  the  user  for  noise  figure,  receiver  circuit 
losses,  and  receiving  antenna  temperature;  and  then  calls 
appropriate  routines  in  FUNCTION  to  compute  the  system 
equivalent  temperature  and  G/T  figure  of  merit. 

XMIT  -  Prompts  the  user  for  output  power  and  transmitter 
circuit  losses,  and  then  calls  an  appropriate  routine  in 
FUNCTION  to  compute  the  effective  isotropic  radiated 
power . 

POWER  -  Organizes  all  the  input  necessary  to  compute  the 
received  useful  power  and  CNR.  First,  the  user  enters 
the  rain  rate,  and  RAINATTEN  is  invoked  to  compute  the 
attenuation  due  to  rainfall.  A I RATTEN  is  also  invoked  to 
compute  the  atmospheric  attenuation.  Next,  the  slant 


range  from  ground  site  to  satellite  is  computed ,  and  the 
total  path  loss  is  determined.  FUNCTION  is  then  invoked 
to  compute  Pf  and  CNR. 

BITRATE  -  Organizes  all  the  input  necessary  to  compute 
available  Ej^/N^  required  E^/N^  margin,  and  maximum  data 
rate.  The  user  enters  the  desired  data  rate  (which  may 
or  may  not  be  achievable),  the  type  of  modulation  to  be 
used,  and  the  desired  bit  error  rate.  FUNCTION  is  then 
invoked  to  compute  available  Ejj/N^  required  E^/N^  and 
max  data  rate;  the  margin  is  computed  by  subtracting  the 
required  Eb/NQ  from  the  available  E^/Nq.  The  module 
CAPACITY  is  then  called  to  determine  whether  the  Nyquist 
or  Shannon  criteria  will  preclude  transmission  at  the 
desired  rate.  If  so,  the  system  will  not  allow  the  user 
to  use  a  bit  rate  which  exceeds  the  capacity  of  the 
channel,  and  either  the  Nyquist  rate  or  the  Shannon  rate 
is  used  instead. 

Main  Program  The  main  program  is  highly  modularized  to 
facilitate  integration  of  the  functional  modules.  Of 
particular  interest  in  the  type  definition  section  is  the 
variable  type  'link',  which  is  the  type  assigned  to  all 
variables  in  the  link  files.  It  consists  of  a  two-part 
record i  one  for  uplink  and  one  for  downlink.  In  addition, 
all  link  file  variables  are  assigned  to  a  single  record,  so 
that  an  entire  link  file  may  be  accessed  by  a  single 
variable.  All  of  the  functional  modules  are  included  for 
compilation  imnediately  before  the  program  body  by  using  the 
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$1  compiler  directive.  Next,  there  is  an  initialisation 
section  where  various  flags  and  strings  are  initialised  and 
files  are  assigned.  The  main  menu  is  then  displayed,  and 
the  user  is  prompted  for  a  menu  item.  The  remainder  of  the 
program  consists  of  11  segments  which  correspond  to  the  11 
main  menu  selections.  Following  completion  of  a  segment, 
control  always  returns  to  the  main  menu  unless  the  EXIT 
option  is  selected. 

It  is  neither  practical  nor  desirable  to  analyse  each  and 
every  line  of  code  in  a  document  of  this  nature.  Instead, 
the  underlying  basis  of  operation  of  each  module  was 
presented  qualitatively.  This  approach  was  chosen  to  give 
the  user  a  better  appreciation  for  the  lower- level  design  of 
the  system. 
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VI.  TESTING  AND  EVALUATION 


This  chapter  addresses  the  general  procedures  followed 
during  system  testing.  It  is  not  the  purpose  of  this 
chapter  to  document  the  exact  testing  procedure  for  every 
minor  program  error  detected,  but  rather  to  give  the  reader 
an  understanding  of  the  overall  testing  process,  and  to 
explain  how  major  program  flaws  were  discovered.  The  end 
product  is  also  evaluated  qualitatively  according  to  the 
criteria  specified  in  the  system  requirements. 

Code  Testing 

This  section  addresses  the  process  of  software  testing, 
starting  with  an  explanation  of  the  overall  testing 
strategy.  Module  testing  and  higher-order  testing  are  then 
discussed  separately. 

Testing  Strategy.  The  system  requirements  chapter 
discussed  in  detail  the  various  types  of  test  case 
formulations,  as  well  as  the  specific  activities  performed 
during  module  and  higher-order  testing.  This  section  does 
not  reiterate  this  material,  but  instead  discusses  the 
overall  testing  strategy  used  to  validate  the  code. 

In  all  cases,  black-box  testing  was  performed  prior  to 
white-box  testing.  In  addition,  by  carefully  choosing  both 
the  type  of  black-box  tests  and  the  order  in  which  modules 
were  tested,  it  was  possible  to  greatly  simplify  the  testing 
process.  In  most  cases,  it  was  unnecessary  to  perform  white 
box  testing  after  black-box  testing,  because  the  process  of 
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black-box  testing  itself  satisfied  all  of  the  conditions  of 
logic  coverage.  In  particular*  exercising  the  input 
screening  routines  (a  type  of  black-box  test)  generally 
guaranteed  statement*  decision*  and  condition  coverage.  The 
main  program  (treated  as  a  separate  module  during  testing) 
was  an  exception  to  this  rule. 

The  functional  packages  were  tested  in  the  following 
order*  utility,  math*  graphics*  low-level*  high-level*  and 
main  program.  However,  the  modules  within  each  package  were 
tested  non-incrementally .  This,  in  effect*  is  a  two-tiered 
testing  mechanism.  This  procedure  allowed  the  modules 
within  each  package  to  be  tested  independently,  with  the  aid 
of  small  driver  modules.  The  only  exception  to  this  was  in 
the  math  package,  where  PWR  had  to  be  tested  before  ALOG 
because  the  latter  calls  the  former. 

The  main  program  itself  was  tested  using  both  black-  and 
white-box  methods.  In  addition,  as  new  modules  were 
completed,  it  was  tested  incrementally  in  a  top-down 
fashion.  This  was  an  ongoing  process  throughout  program 
development . 

Lastly,  function  testing  and  system  testing  were 
performed  to  insure  that  the  final  product  satisfied  all  of 
its  intended  objectives. 

Module  Testing .  Module  testing  was  performed  in  six 
stages*  corresponding  to  the  five  functional  packages  plus 
the  main  program.  The  utility  package  was  tested  first. 
Because  these  modules  are  so  simple*  testing  formalism 
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was  dispensed  with,  and  they  were  tested  simply  by  running 
them  to  see  if  they  worked.  No  errors  were  found  here. 

In  the  math  package,  PWR  and  ALOG  were  tested  in  order 
for  the  reasons  mentioned  above.  All  modules  were  tested 
using  boundary-value  analysis  and  equivalence  partitioning. 
Overflow  errors  occurred  in  BESSEL  due  to  large  floating 
point  values  in  the  series  representation.  This  problem 
could  have  been  alleviated  by  drastically  reducing  the 
number  of  terms  in  the  series,  but  this  would  have  reduced 
accuracy.  The  solution  finally  decided  upon  was  to  break 
both  the  numerator  and  denominator  of  each  term  into  smaller 
pieces,  and  then  perform  arithmetic  on  these  pieces  in  such 
a  way  as  to  avoid  an  overflow  condition.  The  modules  LOG10, 
L0G2,  Q,  and  QINVERSE  all  fail  when  given  a  negative 
argument;  however,  input  screening  prevents  this  situation 
from  occurring.  Inaccuracies  in  the  resulting  output  by  Q 
led  to  the  discovery  of  two  errors  in  the  numerical 
integration  routine.  First,  a  term  was  missing  from  the 
Simpson's  rule  approximation.  Second,  the  curve  being 
integrated  was  not  necessarily  partitioned  into  an  even 
number  of  panels,  as  required  by  Simpson's  rule.  Finally, 
boundary-value  analysis  discovered  that  the  inverse  cosine 
function  in  TRIG  produced  a  division  by  zero  for  an  argument 
of  zero,  as  did  the  inverse  sine  function  for  an  argument  of 
one .  All  these  errors  were  easy  to  correct. 

The  graphics  package  was  tested  next,  since  the  modules 
within  it  call  a  number  of  the  modules  in  the  two  previously 


tested  packages.  The  plots  made  by  HORN  and  DISH  do  not 
vary  with  the  input  due  to  the  ambiguities  that  would  result 
if  the  plots  were  scaled  to  fit  within  the  screen.  In 
HELIX,  stray  lines  unrelated  to  the  gain  pattern  would 
always  appear  before  plotting  began.  This  module  was 
modified  to  suppress  these  lines;  however,  for  some 
combinations  of  input  (usually  near  the  range  limits),  a  few 
points  near  the  x-axis  may  be  omitted.  Lastly,  when  the 
object  code  exceeded  the  64K  limit  set  by  Pascal,  a  memory 
overflow  condition  occurred.  This  problem  was  remedied  by 
converting  the  module  ARRAY  into  an  independent  program  to 
be  compiled  separately  from  the  rest  of  the  system.  It  is 
invoked  in  the  main  program  with  the  'chain'  command. 

Next,  the  low-level  functional  package  was  tested. 

Because  the  modules  in  this  package  are  so  simple  in  design, 
very  few  errors  were  found.  There  were  a  few  problems 
encountered  with  CIRCULAR  and  ELLIPSE  in  particular,  due  to 
the  fact  that  the  sine  and  cosine  functions  expect  the 
argument  to  be  in  radians,  and  degrees  were  sometimes 
mistakenly  used.  Slight  inaccuracies  in  AIR  and  RAIN  were 
traced  to  some  interpolating  polynomials  not  being 
normalized  properly.  Finally,  FUNCTION  contained  a  litany 
of  errors  due  to  inconsistent  units  (e.g.  bandwidth  is 
specified  in  MHz,  while  frequency  is  in  GHz,  etc.). 

Next,  the  six  high-level  functional  modules  were  tested. 
No  problems  were  found  in  REC,  XMIT,  or  BITRATE.  ANTENNA 1 
and  ANTENNA2  presented  some  difficulties.  Both  modules 


initially  had  problems  correctly  storing  the  input 
parameters  associated  with  a  given  antenna.  Because  there 
are  four  different  antenna  types,  each  with  a  different  set 
of  input  parameters,  the  modules  had  to  first  correctly 
identify  the  antenna  type,  and  then  store  the  input 
parameters  in  the  correct  order.  ANTENNA2  also  contained  a 
ubiquitous  logical  error.  While  it  was  able  to  correctly 
compute  the  gain  and  beamwidth  required  for  global  or  spot 
coverage,  there  was  no  reason  to  assume  that  the  satellite 
antenna  was  capable  of  achieving  the  desired  performance. 
Lastly,  it  was  found  that  POWER  failed  when  the  rain  rate 
was  identically  zero.  The  above  anomalies  were  easily 
remedied. 

The  last  module  to  be  tested  was  the  main  program 
itself.  Both  black-  and  white-box  testing  methods  were 
used.  As  new  modules  were  completed  and  tested,  they  were 
incrementally  added  to  the  main  program,  which  was  then 
tested  in  a  top-down  fashion.  Of  course  numerous  minor 
errors  were  found  and  corrected,  but  by  far  the  most 
significant  problem  was  the  implementation  of  the  link 
budget  recomputation  function.  The  difficulty  arose  from 
the  fact  that  the  modules  ANTENNA1  and  ANTENNA2  prompt  the 
user  for  antenna  type,  then  for  various  antenna  parameters, 
and  then  compute  the  gain  and  beamwidth.  However,  when  a 
recomputation  is  performed,  the  antenna  type  has  already 
been  chosen,  and  the  user  is  just  changing  its  original 
parameters.  This  means  that  the  recomputation  software  has 


to  know  the  original  antenna  type,  where  the  original 
parameters  are  stored,  and  whether  it  is  to  work  on  the 
uplink  or  on  the  downlink.  It  took  a  great  deal  of 
restructuring  before  a  suitable  scheme  for  accomplishing 
this  was  devised.  Finally,  it  was  discovered  that  the  order 
in  which  the  inclusion  files  are  listed  at  the  beginning  of 
the  program  is  crucial,  since  this  governs  the  order  of 
compilation.  Consequently,  modules  must  always  be  listed 
before  any  other  modules  which  call  it. 

Higher-Order  Testing .  Two  types  of  higher-order  testing 
were  performed:  function  testing  and  system  testing.  This 
type  of  testing  is  much  less  rigorous  than  module  testing, 
and  tends  to  be  somewhat  subjective.  The  purpose  of 
function  testing  is  simply  to  determine  how  well  the  end 
product  satisfies  the  original  program  requirements.  The 
program  does,  in  fact,  satisfy  the  primary,  secondary,  and 
housekeeping  functional  requirements  as  specified  in  Chapter 
III. 

System  testing,  on  the  other  hand,  is  designed  to 
qualitatively  compare  the  product  to  its  original 
objectives.  First,  every  facility  mentioned  as  a  program 
objective  was  indeed  implemented.  Second,  the  use  of  menus, 
display  of  intermediate  results,  and  color  displays  greatly 
enhanced  the  usability  of  the  system.  Lastly,  the  high 
degree  of  modularity  (especially  the  graphics  driver) 
assures  compatibility  on  other  MS-DOS  machines.  Thus, 
adherence  to  the  original  objectives  is  very  close. 
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Code  Evaluation 


This  section  comments  on  the  merits  of  the  end  product 
from  the  user's  perspective,  as  opposed  to  the  programmer's 
perspective.  Since  the  concept  of  usability  has  already 
been  addressed,  two  additional  aspects  of  the  code  are 
discussed:  accuracy  and  maintainability. 

Accuracy.  One  can  group  all  of  the  modules  which  return 
numerical  results  into  one  of  three  categories:  those  which 
rely  upon  a  formula  for  generating  the  result,  those  which 
rely  upon  numerical  methods,  and  those  which  rely  upon  a 
mathematical  model.  Most  modules  fall  into  the  first 
category,  and  if  the  formula  is  correct  (as  verified  by 
testing),  there  is  really  very  little  room  for  error. 

The  modules  BESSEL,  Q,  and  QINVERSE  comprise  the  second 
category.  BESSEL  evaluates  an  infinite  series,  so  its 
accuracy  depends  upon  the  number  of  terms  included.  Q 
performs  a  Simpson's  rule  numerical  integration,  and  its 
accuracy  depends  upon  the  number  of  panels  that  the 
integrand  function  is  partitioned  into.  QINVERSE  solves  a 
non-linear  equation,  and  its  accuracy  depends  upon  the  value 
of  the  error  tolerance  specified  in  the  code.  It  should  be 
noted  that  this  non-linear  equation  is  based  on  the 
exponential  approximation  to  the  Q-function,  and  thus  has 
inherent  inaccuracies. 

The  modules  AIR  and  RAIN  fall  into  the  last  category. 

Both  of  these  modules  suffer  from  the  inevitable  errors  that 
result  from  ill-fitting  interpolating  polynomials.  While 


the  vast  majority  of  data  points  could  be  connected  by  a 
polynomial,  in  several  cases  the  concavity/convexity  of  the 
curve  was  improper.  In  addition,  both  models  make 
simplifying  assumptions  (see  Appendices  B  and  C)  which  limit 
their  accuracy  somewhat. 

Maintainability.  A  significant  percentage  of  total 
software  cost  is  due  not  to  development,  but  to 
maintenance.  For  this  reason,  software  must  be 
intentionally  designed  to  be  maintainable.  The  high  degree 
of  modularity  of  the  system  greatly  facilitates  expansion 
and  maintenance.  If  a  new  function  is  to  be  added,  then  the 
modules  needed  to  implement  it  are  written  and  tested 


separately,  and  then  included  for  compilation  in  the  main 
program  (taking  care  to  place  them  in  the  proper  order).  A 
stub  is  then  inserted  in  the  proper  place  in  the  main 
program.  Testing  with  the  stub  in  place  assures  that  the 
overall  structure  is  still  sound.  Finally,  the  stub  is 
replaced  by  the  code  required  to  call  the  previously  written 
inclusion  modules.  In  this  manner,  the  system  can  be 
expanded  continuously  without  having  to  alter  the  system 
architecture.  The  system  modularity  also  makes  it  much 
eas'er  to  trace  program  bugs. 

This  chapter  has  addressed  the  issues  of  system  testing 
and  evaluation.  Since  the  program  itself  is  highly  modular, 
an  efficient  testing  process  must  also  be  modular, 
validating  a  piece  of  software  of  this  size  requires  a 
rigorous  and  thorough  testing  plan  in  order  to  assure 
prograsi  reliability. 
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VII.  C0MCLU8 10B8  AMD  RKONMUDATIOM 


The  resultant  and  product  of  this  research  effort  is  a 
comprehensive,  integrated,  and  transportable  software 
package  for  the  analysis  of  digital  satellite  links.  While 
the  prisiary  objective  of  this  software  package  is  link 
budget  analysis,  the  package  also  includes  routines  for 
coeluting  various  satellite  orbital  parameters,  as  well  as 
for  plotting  antenna  gain  patterna.  All  of  the  noeiinal 
functions  specified  in  the  system  requi resient a  have  been 
implemented . 

System  portability  to  other  MS-DOS  machines  is 
facilitated  by  a  rearcvable  graphics  driver.  Expandability 
is  straightforward!  a  new  function  is  written  as  a  separate 
module  and  tested  independently:  it  is  then  included  for 
compilation  in  the  s^in  program  via  the  51  cospller 

directive.  In  addition,  the  problem  of  tracing  perceived 
errors  is  mitigated  by  the  display  of  intermediate  results 
following  each  major  calculation. 

The  man-machine  interface  is  designed  to  maximise  the 
user's  efficiency  as  well  as  to  minimise  confusion.  Menus 
are  used  extensively  to  guide  the  user  through  the  program, 
and  the  use  of  prompting  and  input  validity  checking  insure 
the  accuracy  of  all  data  entries.  In  addition,  color 
displays  are  used  to  present  the  data  in  an  organised,  easy- 
to-read  format . 

The  system  has  successfully  withstood  the  rigors  of  both 
module  testing  and  higher-order  testing.  It  is  highly 


improbable  that  tha  usar  can  cauaa  tha  system  to  craah  by 
entering  inappropriate  data. 

Thia  softwara  package  la  far  froai  coaiplata.  Thara  ara  a 
n unbar  of  recoamendations  for  futura  Modification  and 
expanaion  t 

1.  Davalop  a  mathematical  nodal  for  calculating 
antanna  noiaa  tanparatura.  Tha  graataat  challanga  to 
thia  will  ba  computing  tha  aky  noiaa  tanparatura 
coaiponant ,  which  la  a  function  of  both  fraquancy  and 
alavation  angla. 

2.  Examine  tha  af facta  of  dapolarisat ion ,  ica 
crystala,  multipath  propagation,  and  atmoapharic 
acinti 1 lation ;  and  dataraina  if  thay  can  ba  quantised 
and  incorporatad  into  tha  rain  and  air  attanuation 
modal  a . 

3.  Apply  aoaia  of  tha  raaulta  of  Navaa  (6),  ao  that 
link  availability  may  ba  aatimatad  baaed  on  local  rain 
rate  atatiatica. 

4.  Give  the  uaar  the  option  to  employ  pulaa  ahaping 
and  error  correcting  codea . 

5.  Allow  multiple  accaaaea  to  the  link,  and  introduce 
a  noiaa  term  due  to  intarfarance  from  other  uaara. 

Tha  incorporation  of  theae  recommendations  should  make 
tha  packaga  much  more  flexible,  and  therefore  more 
applicable  to  real-world  design  problems. 
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User ' a  Manual  ] 

i 

i 

While  the  program  is  highly  modular  and  highly  menu- 
driven,  there  are  a  few  points  regarding  its  operation  which 

i 

require  amplification.  It  is  not  the  purpose  of  this  manual  < 

I 

to  teach  link  analysis,  but  rather  to  show  the  user  how  to  j 

i 

use  the  program  correctly  and  efficiently.  j 

The  host  computer  must  be  an  IBM  PC  or  equivalent  ! 

i 

compatible  with  at  least  one  floppy  disk  drive.  A  dot 
matrix  printer  is  also  suggested  for  obtaining  hardcopy  of 
plots  or  other  displays.  The  system  disk  contains  36  source 

l 

files,  which  occupy  92K.  When  compiled,  the  object  code 

i 

occupies  63K.  Since  DOS  requires  an  additional  27K,  it  is  | 

recommended  that  the  host  computer  be  equipped  with  at  least 

I 

128K  of  memory.  Also,  the  displays  will  only  be  readable  on  | 

I 

a  color  monitor.  | 

I 

If  only  one  disk  drive  is  available,  then  the  system  disk 
and  the  data  disk  will  be  one  in  the  same.  Otherwise,  it  is  j 

I 

recommended  that  the  data  disk  reside  on  a  separate  drive.  j 

I 

This  is  especially  useful  on  machines  with  Winchester-type 
hard  drives,  which  can  hold  very  large  link  analysis 

databases.  The  program  initially  assumes  drive  A  to  be  the  I 

default  data  drive,  until  told  otherwise  by  the  user.  At 

any  rate,  it  is  imperative  that  the  path  drive »data\  exists 

on  the  data  disk,  since  this  constitutes  the  directory  under 

which  all  link  files  are  stored.  Failure  to  include  this 

path  on  the  data  disk  will  result  in  a  fatal  I/O  error. 
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To  start-up  the  system,  boot  the  operating  system,  and 
then  place  the  system  disk  in  drive  A  and  type  'run.'  The 
main  menu  will  then  appear  (see  Figure  A-l).  The  default 
drive  should  be  switched  now  if  desired.  A  listing  of  all 
link  files  on  the  default  data  disk  can  be  obtained  by 
choosing  menu  option  9,  'List  Link  File  Directory.' 

Unwanted  files  may  be  deleted  by  choosing  menu  option  8, 
'Delete  Link  File,'  and  specifying  the  filename  to  be 
deleted  without  the  .Ink  extension. 

Menu  item  1,  'Satellite  Orbital  Parameters,'  contains  a 
submenu  which  permits  the  user  to  choose  from  computing 
escape  velocity,  satellite  steady-state  temperature, 
circular  orbit  parameters,  elliptical  orbit  parameters,  or 
azimuth  and  elevation  angles.  User  input  in  this  section  is 
very  straightforward.  The  circular  orbit  parameters  include 
satellite  velocity,  orbital  period,  coverage  area,  and 
duration  of  visibility.  The  elliptical  orbit  parameters 
include  satellite  velocity  at  apogee  and  perigee,  average 
velocity,  velocity  at  a  given  reference  angle  phi  (where 
phi-0  at  perigee  and  phi-pi  at  apogee),  orbital  period,  and 
coverage  area.  The  elevation  angles  computed  in  this 
section  may  be  used  later  on  during  link  analysis.  This  is 
the  only  connection  to  the  rest  of  the  program. 

Main  menu  item  2,  'Antenna  Gain  Patterns,'  also  has  no 
connection  to  the  rest  of  the  program.  It  contains  a 
submenu  from  which  the  user  may  choose  to  plot  the  gain 
patterns  of  horn,  helical,  parabolic  reflector,  or  phased 


MAIN  MENU 


1.  Satellite  Orbital  Parameters 

2.  Antenna  Gain  Patterns 

3.  Link  Analysis 

4.  Display  Link  Budget 

5.  Recompute  Link  Budget 

6.  Store  Link  File 

7.  Recall  Link  File 

8.  Delete  Link  File 

9.  List  Link  File  Directory 

10.  Change  Default  Drive 

11.  EXIT 


Figure  A-l.  Main  Menu  Display 


IV  J 

L*  V 


% 


y 


•V .  j 


p'. 


array  antannas.  Tha  firat  thraa  aalactiona  raquira  tha  uaar 
to  antar  tha  frequency  and  varioua  itmmm  relating  to  tha 
geometry  of  tha  particular  antenna,  and  thair  uaa  ia  alao 
at raight forward .  Tha  phaaad  array  aalaction,  however,  ia 
conaidarably  nora  coaiplicatad . 

Thia  aalaction  plota  tha  gain  pattarn  for  a  planar  phaaad 
array  of  idantical  dipolaa.  Firat,  tha  uaar  aiuat  antar  tha 
langth  of  tha  individual  dipole  a lament a  normal lead  to  tha 
wavelength.  Next,  tha  uaar  muat  chooae  from  ona  of  aavan 
currant  diatributiona  on  tha  dipolaa.  Tha  daairad  array 
muat  than  ba  dacompoaad  into  homoganaoua  aubarraya  in  both 
tha  x-  and  y-directiona .  Tha  program  prompta  tha  uaar  for 
tha  numbar  of  x-  and  y-aubarraya,  tha  numbar  of  alamanta  in 
aach  aubarray,  tha  apacing  between  thoaa  alamanta  normalisad 
to  wavelength,  and  tha  phaaa  tapar  batwaan  conaacutiva 
alamanta.  For  many  caaea,  tha  aaparationa  and  phaaa  tapara 
batwaan  any  two  adjacant  alamanta  may  wall  ba  idantical; 
conaequantly,  thara  will  only  ba  ona  aubarray  in  aach 
diraction.  Givan  thia  data,  tha  program  will  ganarata  aix 
polar  plota,  thoaa  baing  tha  thata  and  phi  componanta  of  tha 
planaa  theta«pi/2  (xy-plana),  phi-0  (xz-plana),  and  phi«pi/2 
(ys-plana).  After  aach  plot,  tha  uaar  ia  prompted  for  a 
magnification  factor.  Entering  a  value  greater  than  unity 
will  cauaa  tha  aama  gain  pattarn  to  ba  plotted,  only  larger 
than  before.  Tha  converae  holda  for  valuaa  laaa  than 
unity.  Entering  a  value  of  zero  will  cauaa  tha  acreen  to  be 
cleared,  and  tha  naxt  polar  plot  to  ba  drawn. 


Main  M(iu  it  mu  3,  ‘Link  Analysis,  *  is  ussd  to  psrfonn 
link  budget  analysis.  Ths  ussr  is  guided  stsp  by  stsp 
through  svsry  calculation,  and  intansadiats  rasults  ara 
displayad  along  tha  way.  Tha  final  link  budgat  can  ba 
displayad  by  aalacting  menu  itM  4,  'Display  Link  Budgat.' 
Manu  itasi  5,  '  Recosq>ute  Link  Budgat,'  allows  tha  usar  to 
racoagwta  all  tha  link  calculations  with  naw  valuas  of 
fraquancy,  antanna  dasign,  transaiittar  powar,  racaivar  noisa 
figura,  bit  arror  rata,  or  data  rata.  It  also  pannits  tha 
usar  to  parforai  tha  ravaraa  procass  of  coaq>uting  bit  arror 
rata  givan  tha  data  rata.  One a  a  das i rad  link  budgat  is 
obtainad,  tha  usar  may  stora  it  on  tha  data  disk  by 
salacting  manu  ltam  6,  'Stora  Link  Fila.'  Praviously  storad 
link  filas  may  ba  racallad  for  raviaw  or  racomputation  by 
salacting  manu  itam  7,  'Racall  Link  Fila.’ 

Tha  following  is  an  ordarad  listing  of  instructions  for 
dasigning  and  analysing  a  typical  link,  starting  from  tha 
main  manu . 

1.  Entar  tha  fraquancias  for  tha  uplink  and  downlink. 

2.  Entar  tha  alavation  anglas  and  latitudes  for  the 
transmitting  and  receiving  ground  sites.  Elevation 
anglas  praviously  computed  in  the  satellite  parameters 
section  may  ba  used  again  here. 

3.  Enter  the  mean  sea  level  altitude  of  the  transmitting 
and  receiving  ground  sites. 

4.  Enter  the  satellite  altitude  (enter  'G'  for 
geosynchronous  altitude). 


5.  Enter  the  system  bandwidth  (The  signal  bandwidth  is 
taken  to  be  equal  to  the  noise  bandwidth). 

6.  Specify  the  type  of  antenna  for  the  uplink 
transmitter. 

7.  Enter  the  appropriate  antenna  parameters  asked  for. 

8.  Specify  the  type  of  coverage  for  the  uplink  receiver, 
along  with  the  antenna  parameters  asked  for. 

9.  Repeat  steps  6-8  for  the  downlink  transmitter  and 
receiver.  The  maximum  gain  and  3-dB  beamwidth  for  each 
antenna  is  displayed  as  an  intermediate  result. 

10.  Enter  the  output  amplifier  power  level  and 
transmitter  circuit  losses  for  the  transmitters  at  the 
first  ground  site  and  onboard  the  satellite.  The 
effective  isotropic  radiated  power  is  then  calculated  and 
displayed . 

11.  Enter  the  amplifier  noise  figure,  circuit  losses,  and 
antenna  noise  temperature  for  the  receivers  at  the  second 
ground  site  and  onboard  the  satellite.  The  effective 
system  temperature  and  G/T  figure  of  merit  are  then 
displayed. 

12.  Enter  the  rain  rate  at  the  uplink  and  downlink 
sites.  If  it  is  not  known  exactly,  choose  one  of  the 
qualitative  descriptors  that  best  matches  the  estimated 
rain  rate. 


13.  Specify  the  desired  data  rate,  modulation  type,  and 
bit  error  rate  for  the  uplink  and  downlink  signals. 


After  a  slight  delay,  the  available  Eb/No,  required 


Eb/N0,  and  maximum  data  rate  will  be  displayed. 
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14.  This  is  the  end  of  user  data  entry.  The  link  budget 
may  now  be  displayed  by  selecting  menu  item  4,  'Display 
Link  Budget.'  The  display  consists  of  two  columns:  one 
for  the  uplink  and  one  for  the  downlink.  The  display  is 
also  organized  into  color-coded  groups  for  convenience. 
Cyan  corresponds  to  general  information,  blue  to  the 
transmitter,  yellow  to  the  channel,  green  to  the 
receiver,  magenta  to  link  performance  criteria,  and  red 
to  critical  link  constraints.  If  the  computed  data  rate 
exceeds  the  maximum  data  rate,  then  it  will  blink.  The 
same  holds  for  negative  power  margins.  This  is  to  call 
attention  to  the  fact  that  the  link  cannot  exist  as 
specified . 

15.  The  user  may  now  select  menu  item  5,  'Recompute  Link 
Budget'  in  order  to  tailor  the  system  to  his  needs. 

16.  The  final  link  budget  may  be  stored  on  disk  for  by 
selecting  menu  item  6,  'Store  Link  File.'  Previously 
computed  link  files  may  be  retrieved  for  review  or 
additional  computation  by  selecting  menu  item  7,  'Recall 
Link  File. ' 

This  manual  should  enable  the  user  to  efficiently  use  the 
link  analysis  program.  Because  all  user  input  is  screened, 
it  is  doubtful  that  erroneous  input  can  cause  the  system  to 
crash.  In  the  unlikely  event  that  this  does  occur,  the  user 
should  restart  the  program.  The  only  data  potentially  lost 
will  be  that  data  entered  subsequent  to  the  last  'Store  Link 
File1  command. 
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APPENDIX  B 


Mathematical  Model  for  Determining  Atmospheric  Attenuation 

At  sufficiently  high  frequencies,  atmospheric  attenuation 
may  significantly  degrade  a  satellite  link.  This 
attenuation  is  due  to  absorption  of  electromagnetic  waves  by 
various  constituent  gases  in  the  atmosphere.  These 
interactions  occur  at  the  gas  resonance  frequencies.  In  the 
1-100  GHz  frequency  band,  the  only  significant  contributors 
to  atmospheric  attenuation  are  water  vapor  and  oxygen.  The 
water  vapor  absorption  band  is  centered  at  22.235  GHz,  while 
the  oxygen  absorption  lines  extend  from  53.5  to  65.2  GHz. 
Satellite  communications  bands  are  generally  chosen  to  avoid 
these  frequencies,  so  resulting  attenuation  values  are 
typically  less  than  1  dB  (2x94-95;  9x323). 

The  total  atmospheric  attenuation  depends  upon  a  number 
of  variables;  namely  frequency,  elevation  angle,  ground 
station  altitude,  relative  humidity,  and  surface 
temperature.  Of  these,  only  frequency  and  elevation  angle 
play  a  significant  role.  Since  the  uncertainty  of  the  link 
may  be  on  the  order  of  several  dB  anyway,  it  is  perfectly 
acceptable  to  neglect  the  other  variables.  Hence,  the  model 
described  below  determines  atmospheric  attenuation  as  a 
function  of  frequency  and  elevation  angle. 

Figure  B-l  is  a  logarithmic  plot  of  atstospheric 
attenuation  versus  frequency  for  a  vertical  path  over  the  1- 
1N  GHb  range.  Tike  mathemat  i ca  1  model  for  atsiospheric 
attenuation  was  developed  by  fitting  approximating 


polynomials  to  various  sections  of  this  curve.  The 
equations  of  the  polynomials  were  obtained  via  a 
commercially  available  curve-fitting  program.  Thus,  a 
quantitative  relationship  between  attenuation  and  frequency 
was  established.  Because  of  the  complexity  of  this  curve, 
it  was  partitioned  into  four  regions:  1-22  GHz,  22-40  GHz, 
40-60  GHz,  and  60-100  GHz. 

Figure  B-2  shows  the  first  region.  Note  that  both  axes 
have  logarithmic  scales,  even  though  attenuation  is  in  units 
of  dB.  While  taking  the  logarithm  of  a  quantity  already 
expressed  in  dB  may  not  make  sense  intuitively,  it  was 
necessary  in  order  to  generate  the  curves  properly.  A 
single  approximating  polynomial  joining  all  the  data  points 
in  this  region  could  not  be  found,  so  the  region  was 
decomposed  into  three  subregions.  The  first  two  subregions 
correspond  to  areas  of  the  curve  where  the  relationship 
between  attenuation  and  frequency  is  approximately  linear, 
so  two  different  linear  regression  lines  are  used  to  model 
this  relationship.  The  last  subregion  corresponds  to  an 
area  of  approximately  cubic  data  dependence,  so  a  third- 
degree  polynomial  is  used. 

Figure  B-3  shows  the  second  region.  The  data  points  in 
this  region  display  a  good  quadratic  relationship,  so  a 
second-degree  polynomial  provides  a  good  approximation. 

Figure  B-4  shows  the  third  region.  Several  higher-degree 
potynoeiials  were  fit  to  this  data,  but  only  the  second- 
degree  polynomial  provided  an  acceptable  fit.  Note  that 
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there  is  some  measurable  error  along  the  center  of  the 
curve. 

Figure  B-5  shows  the  fourth  region.  This  region  also  had 
to  be  decomposed  into  three  subregions,  and  proved  to  be  the 
most  difficult  to  model.  Only  two  data  points  occupy  the 
first  subregion,  since  no  way  could  be  found  to  connect 
these  points  with  the  rest  of  the  curve.  These  points  are 
joined  by  a  regression  line.  Similarly,  the  second 
subregion  is  also  occupied  by  only  two  points.  However,  a 
quadratic  polynomial  is  used  to  connect  them  in  order  to 
more  closely  approximate  the  true  shape  of  the  curve.  The 
last  subregion  is  modeled  quite  nicely  by  a  fourth-degree 
polynomial . 

The  above  discussion  has  explained  in  some  detail  how  the 
mathematical  model  for  vertical  path  atmospheric  attenuation 
was  developed.  To  include  the  elevation  angle  dependence, 
it  is  necessary  only  to  multiply  the  vertical  path 
attenuation  by  the  cosecant  of  the  elevation  angle. 


Figure  B-l  .  Atmoepher  i  ,  a*  • 
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APPENDIX  C 


Mathematical  Model  for  Determining  Rain  Attenuation 

Rainfall  is  the  most  serious  degrading  factor  in  a 
satellite  link.  Rain  droplets  scatter  and  absorb  impinging 
electromagnetic  radiation,  causing  potentially  severe  signal 
attenuation,  particularly  at  wavelengths  approaching  the 
size  of  the  droplets.  These  interactions  depend  upon  the 
number  of  droplets  encountered  (and  thus  upon  the  rain  rate) 
as  well  as  the  geometry  of  the  droplets. 

There  is  a  wide  variety  of  rain  attenuation  models  from 
which  to  choose.  However,  many  of  them  rely  heavily  upon 
the  knowledge  of  extensive  local  rain  rate  statistics  or 
rain  droplet  distributions,  which  are  unavailable  to  many 

users.  The  model  presented  here  estimates  rain  attenuation 
as  a  function  of  frequency,  elevation  angle,  rain  rate, 
ground  station  altitude,  and  ground  station  latitude. 

The  attenuation  A  due  to  rainfall  has  been  commonly  given 
in  the  literature  by 

A  ■  aRbL  dB  (c-1 ) 

where  a  and  b  are  frequency-dependent  coefficients,  R  is 
rain  rate  in  mm/hr,  and  L  is  effective  path  length.  The 
quantity  aRb  by  itself  is  called  the  specific  attenuation, 
and  has  units  of  dB/km.  Note  that  the  attenuation  depends 
on  the  rate  of  rainfall,  rather  than  on  rainfall 
accumulation.  The  process  of  measuring  rain  rate  is  inexact 
at  best,  so  the  model  allows  the  user  to  qualitatively 


describe  the  rain  intensity  (light,  medium,  etc.)*  and  then 
inserts  a  ballpark  value.  Hence,  calculating  attenuation 
involves  two  basic  steps  *  determining  the  coefficients  a  and 
b,  and  then  determining  the  effective  path  length  (2:96-99; 
9:327-335). 

Approximations  for  the  coefficients  a  and  b  are  given  in 
[7]  as  functions  of  frequency  only.  However,  these 
coefficients  are  also  dependent  upon  the  microstructure  of 
the  rain,  such  as  the  rain  temperature,  droplet  shape,  and 
droplet  size  distribution.  With  a  little  extra  work,  it  is 
possible  to  obtain  a  very  accurate  approximation.  In  [7], 
there  is  extensive  tabular  data  for  the  coefficients  a  and  b 
up  to  1000  GHz,  for  rain  temperatures  of  20°C  and  0°C,  for 
several  different  droplet  distributions.  These  droplet 
distributions  include  those  of  Laws  and  Parson,  Marshall  and 
Palmer,  and  Joss  et  al.  While  it  is  beyond  the  scope  and 
intent  of  this  report  to  present  these  distributions  in 
detail,  it  should  be  mentioned  that  the  Laws  and  Parsons 
distribution  was  chosen  because  it  is  a  little  more  accurate 
than  the  others  above  30  GHz.  This  distribution  also  yields 
different  a  and  b  values  for  low  and  high  rain  rates.  As 
for  rain  temperature,  while  the  assumption  of  20°C  gives 
suitable  results  for  terrestrial  links,  0°C  is  more 
appropriate  for  satellite  links.  All  of  the  figures  in  this 
appendix  were  obtained  by  plotting  the  tabular  data  for  the 
coefficients  a  and  b  (12:1466-1467). 

Figure  C-l  is  a  plot  of  a(f)  versus  frequency  for  low 
rain  rates  (i.e.  less  than  25  nn/hr),  and  is  approximated  by 


,<^r 


a  sixth-degree  polynomial.  Figure  C-2  is  simply  a 
magnification  of  the  previous  figure  for  frequencies  below 
10  GHz,  and  is  also  modeled  by  a  sixth-degree  polynomial. 
Figure  C-3  is  a  plot  of  a(f)  versus  frequency  for  high  rain 
rates  (i.e.  greater  than  25  nun/hr),  and  is  approximated  by  a 
tenth-degree  polynomial.  Figure  C-4  is  also  a  magnification 
of  the  previous  figure  for  frequencies  below  10  GHz,  and  is 
given  by  a  cubic. 

Figure  C-5  is  a  plot  of  b(f)  versus  frequency  for  low 
rain  rates.  Due  to  the  nature  of  this  particular  curve,  a 
single  approximating  polynomial  could  not  be  found.  Hence, 
the  curve  was  divided  into  two  regions:  1-8  GHz,  and  8-100 
GHz.  Figure  C-6  shows  the  section  of  the  curve  from  1-8 
GHz,  which  is  approximated  by  a  fifth-degree  polynomial. 
Figure  C-7  shows  the  section  of  the  curve  from  8-100  GHz, 
which  is  approximated  by  an  eighth-degree  polynomial. 

Figure  C-8  is  a  plot  of  b(f)  versus  frequency  for  high  rain 
rates.  Due  to  the  considerable  complexity  of  this  curve,  it 
was  broken-up  into  three  regions:  1-6  GHz,  6-30  GHz,  and  30- 
100  GHz.  Figure  C-9  shows  the  section  of  the  curve  from  1-6 
GHz,  which  is  approximated  by  a  sixth-degree  polynomial. 
Figure  C-10  shows  the  section  of  the  curve  from  6-30  GHz, 
which  is  approximated  by  a  fifth-degree  polynomial.  Lastly, 
Figure  C-ll  shows  the  section  of  the  curve  from  30-100  GHz, 
which  is  approximated  by  fourth-degree  polynomial. 

Once  the  coefficients  a  and  b  have  been  determined,  the 
effective  path  length  through  the  rain  must  be  found.  For  a 
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given  elevation  angle  El,  the  path  length  L  through  the  rain 
is  given  by 

L  *  (H  -  HQ)/sin(El)  km  (C-2) 

where  H  is  the  effective  storm  height  in  km,  and  H  is  the 
e  o 

ground  station  elevation  in  km. 

The  temperature  gradient  where  T=0°C  is  called  the  zero- 
degree  isotherm,  and  is  generally  the  point  above  which  no 
liquid  water  exists.  The  height  of  the  zero-degree  isotherm 
H^,  (also  known  as  melting  layer  height),  is  related  to  the 
latitude  by 

Hi  =  4.8  km  abs(lat)<=30°  (C-3) 

*  7.8  -  0.1*abs(lat )  abs(lat)>30° 

Seasonal  variations  of  the  melting  layer  height  have  been 
neglected.  The  effective  storm  height  is  related  to  the 
melting  layer  height  by 

He  *  R<=10  nro/hr  (C-4) 

=  +  log(R/10)  R>10  nm/hr 

Note  that  C-4  accounts  for  the  fact  that  higher  rain  rates 
are  associated  with  more  active  convective  rain  cells,  which 
contain  strong  updrafts  that  can  carry  water  well  above  the 
zero-degree  isotherm. 

The  above  relations  are  generally  accurate  for  relatively 
low  rain  rates  (less  them  10  mm/hr).  However,  high  rain 
rates  create  another  wrinkle  in  the  problem.  The  same 


convective  rain  cells  which  carry  water  above  the  zero- 
degree  isotherm  also  render  the  distribution  of  rain  non- 
uniform  in  the  horizontal  direction.  Much  data  on  this 
problem  has  been  collected  using  rain  gauge  networks,  but 
these  experiments  are  beyond  the  scope  of  this  report.  The 
end  result  of  this  research  was  the  derivation  of  a  path 
profile  model  which  takes  the  presence  of  convective  rain 
cells  into  account.  Putting  all  the  pieces  together,  the 
total  path  attenuation  A  due  to  rain  is  given  by 

A  =  aRbL  R<=10  mm/hr  (C-5) 

A  =  aRb[l  -  exp{-PL}]/P  R>10  iran/hr 

where  P  =  gb  ln(R/10)cos  El 
and  g  =  1/22 

By  specifying  the  rain  rate  as  a  value  which  occurs  a  given 
percentage  of  time,  the  resultant  attenuation  will  also  be 
expressed  as  a  percent  time  occurrence.  Likewise,  if  a  link 
is  designed  to  operate  with  a  fixed  amount  of  rain 
attenuation,  then  the  link  availability  can  be  determined 
from  the  percentage  of  time  that  rain  rate  is  exceeded 
(12il466-1471) . 

The  above  discussion  has  presented  the  development  of  the 
rain  attenuation  model.  It  should  be  noted  that  this  is 
only  a  crude  model,  and  does  not  take  into  account  such 
factors  as  ice  crystals,  fog,  suspended  particulates,  and 
depolarization . 
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Frtqutncy  (GHz) 


Figure  C-l.  A(f)  vs.  Frequency  for  Low  Rain  Rates  (7:323) 
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Fr*qu«ncy  (GHz) 


Figure  C-2.  Magnification  of  Frequencies  Below  10  GHz 
(7:323) 


Figure  C-5.  B(f)  vs.  Frequency  for  Low  Rain  Rates  (7j323) 
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Frequency  (8Hx) 


Figure  C-10.  6-30  GHz  Partition  (7*323) 

K' 
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APPENDIX  D 


Source  Code  LI s 1 1 ng 


(*  SATELLITE  LIMC  ANALYSIS  WORKSTATION  FOR  IBM  PC  (Ver  1.1)  *) 

(*  Copyright  1986  by  Bill  Jhckscn  *) 

program  master(  input  .output); 

const  mu=3.964el4; 
fe=6.378e6; 
dr“1.7463e-2; 
rd=57.2958; 
k=1.38e-23; 

<H>  ********  «  0  *  *********  *********** ) 


( 

type  legal =set  of  char; 
f  i  lename=stri  ng[23] ; 
1 1nk=record 


(*  gravitational  constant  times  Earth  mass  *) 
(*  Earth  radius  *) 

(*  deg  rees-to- radians  conversion  *) 

(*  radians-to- degrees  conversion  *) 

(*  Boltanan's  constant  *) 


upl  :real ; 
dnl  :real ; 

end; 

1  inkrecord=record 

altitude  (*  satellite  altitude  *):real; 
lat  (*  ground  station  latitude  *), 
rt>  (*  ground  station  MSL  elevation  *) , 
alpha  (*  elevation  angle  *), 
freq  (*  frequency  *), 

TxEfficlency  {*  transmitting  antenna  efficiency  *), 
fa  Efficiency  {*  receiving  antenna  efficiency  *), 
TxG&in  (*  transmitting  antenna  gain  *), 
faGain  (*  receiving  antenna  gain  *), 
beamwidth  (*  3-dB  beamwidth  *), 
tout  (*  transmitter  output  power  *) , 

TxClixultLoss  (*  transmitter  circuit  losses  *), 
faCIrcuitloss  (*  receiver  circuit  losses  *), 

EIRP  (*  effective  Isotropic  radiated  power  *), 
NolseFIgure  (*  noise  figure  *), 
faArtfemp  (*  receiving  antenna  noise  temperature  *), 
Teq  (*  systan  equivalent  tarperature  *), 

GT  (*  (VT  figure  of  merit  *), 

Slanttenge  (*  slant  range  to  satellite  *), 
ralnrate  (*  rein  rate  *), 

RainAtten  (*  rainfall  attenuation  *), 

AirAtten  (*  atmospheric  attenuation  *), 

Pathless  (*  free  space  propagation  loss  *) , 

Pr  (*  received  useful  power  *), 

CJR  (*  carrier-to-noise  power  density  ratio  *), 
bandwidth  (*  system  bandwidth  *), 
datarate  (*  transmitted  data  rate  *) , 
modtype  (*  modulation  type  *), 
m  (*  renter  of  signalling  levels  *), 
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BER  (*  bit  error  rate  *), 

Avail  B>fto  (*  available  Q)/Nd  *), 

ReqEbto  (*  required  Eb/No  *), 
tex  Data  rate  (*  maximum  data  rate  *) , 
margin  (*  link  power  margin  *) , 

Antenna  Type  (*  antenna  type  *) , 
beamtype  (*  beam  type  *) , 

angl,arg2,arg3,arg4,arg5  (*  dunny  variables  ♦) :1  ink; 

end; 


var  InputError  (*  illegal  naneric  input  *):boolean; 
integerdata  (*  integer  input  data  *)  integer; 
real data  (*  real  input  data  *)  :real ; 
legal  integer  .legal  real  (*  valid  input  sets  *):legal; 
duiriun  (*  dunny  file  *):file; 

data  (*  link  file  record  containing  all  pertinent  variables  *):linkrecord; 

linkfile  (*  link  file  *):file  of  linkrecord; 

buffer  (*  storage  buffer  for  use  in  deleting  link  files  *), 

library  (*  directory  of  current  link  files  *):file  of  string[23]; 

1  (*  loop  variable  *) , 
option  (*  menu  item  option  *), 
choice  (*  submenu  item  choice  *) , 
site  (*  antenna  site  *):integer; 

Ebto  (*  temporary  storage  for  Eb/No  *), 

rate  (*  tanporary  storage  for  datarate  *) , 

err  (*  temporary  storage  for  bit  error  rate  *} , 

limit  (*  limiting  data  transnission  rate  *):real; 

junk  (*  durmy  variable  *)  :real ; 

path  (*  link  file  directory  path  *):string[ll]; 

oldiame,name,dujmyname  (*  link  filename  temporary  storage  *)  :string[23]; 

letter  (*  dunny  variable  *):char; 

inp  (*  dunrny  variable  *)  :string[8]; 

killnane  (*  name  of  file  to  be  deleted  *}:string[8]; 

havealpha  (*  elevation  angles  previously  computed  *), 

skip  (*  skip  reentry  of  elevation  angles  in  link  analysis  *), 

recompute  (*  it  is  possible  to  recompute  the  existing  link  budget  *), 

Already  There  (*  filerame  already  exists  *), 

uplink  (*  uplink  portion  is  being  confuted  *), 

newinfo  (*  new  unsaved  data  has  neen  generated  *) , 

standard  (*  FALSE  when  computing  BER  given  datarate  *), 

shannon  flag  (*  indicates  Sermon  rate  Is  limiting  factor  *), 

nyquistflag  (*  indicates  fyqulst  rate  Is  limiting  factor  *):boolean; 

title  (*  display  titles  *)  array  [1..10]  of  string[20]; 


§$I  A:1nteger.past 
§$1  A:real.past 
§$I  A:ex1st.past 
§$I  A:tr1g.past 
§$I  Arpwr.past 
§$I  A:bessel .past 
§$I  A:logl0.past 


§$I  Aalog.past 
§$I  A:log2.past 
§$I  Aralse.past 
§$I  A:custan.past 
§$I  A  .-wait,  past 
§$I  A:getalt.past 
§$I  A:escape.past 
§$I  Arblackbod.past 
§$I  Axircular.past 
§$I  A:ell ipse. past 
§$I  A:azel  .past 
§$I  A:q.past 
§$I  A:qinverse.past 
§$I  A:ilinction.past 
§$I  A:hom.past 
§$I  Arhelix.past 
§$I  Aidish.past 
§$I  A:antennal.past 
§$I  A  antenna  2.  past 
§$I  A:xmit.past 
§$I  A:rec.past 
§$I  A:rain.past 
§$I  A:a1r.past 
§$I  A:poMer.past 
§$I  A:save1t.past 
§$I  Arcapaclty.past 
^$1  Aibitrate.past 

(*  INITIALIZATION  AN)  SET-UP 

begin 

title[l]:*'UPLINK  TRANSMITTER'; 
title[2]:*'UPLINK  RECEIVER' ; 
title[3]:='D0MN.INK  TRANSMITTER'; 
title[4]:=' DOWNLINK  RECEIVER'; 
title[5]:*' GROUND  TERMINAL'; 
title[6]:='SATELLITE' ; 
title[ 7] UPLINK  SITE'; 
t itle[8] :* ' D0WN.INK  SITE'; 
tit1e[9]:*'U  P  L  I  N  K' ; 
title[10]:*'D  0  W  N  L  I  N  K' ; 
path:*'a:data0'; 
legalinteger :*[' O' . . ' 9' ] ; 
legalreal  :«C'0'..'9' ,' .' .'+'  ,'e’  ,'E']; 
havealpha:*false; 
skip:»false; 
recompute:  “false; 

A1  ready  The  re  :*fa1  se ; 
newinfb:*false; 
standard  :«true; 
ass1gn(dunduB,'aamay.chn' ) ; 
assign(  library  .path*’ link  .dir'); 
assign(buffer  jath+'tanp.sto* ) ; 
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vj.vv  *  -■ 


if  not  exist(path+*  link  .dir')  then 
begin 

rewrite( library); 
close(library); 

end; 


repeat 

textmode; 

textbackground(green) ; 
clrscr; 
gotoXY(30,5); 
textbackground(yel  low) ; 
textcolor(black) ; 

writeln('  MAIN  MENU':19); 
textbackground(green) ; 
textcolor(viiite); 
gotoXY(l.lO); 

writeln('  1.  Satellite  Orbital  ferameters'); 
writelnC  2.  Antenna  Gain  Patterns'); 

writelnC  3.  Link  Analysis'); 

writeln('  4.  Display  Link  Budget'); 

writeln('  5.  Recompute  Link  Budget'); 

writeln('  6.  Store  Link  File'); 

writelnC  7.  Recall  Link  File'); 

writeln('  8.  Delete  Link  File'); 

writelnC  9.  List  Link  File  Directory'); 

writelnC  10.  Change  Default  Drive'); 

writelnC  11.  EXIT'); 

repeat 

writeln  jwriteln ;write( '  ENTER  OPTION  :'); 

Get  Integer; 

option  :=integerdata ; 

if  not  (option  in  [1..11])  then  lnputError:=true; 

if  Input Error=true  then  noise; 
until  InputError=false; 
case  option  of 

l:begin  (*  Satellite  Orbital  farameters  *) 
textbackgnound( red) ; 
cl  rscr  * 
gotoXY(l.lO); 

writelnC  1.  Escape  Velocity'); 

writelnC  2.  Steady-State  Tsnperature'); 

writeln('  3.  Circular  Orbit  ft rameters' ); 

writelnC  4.  Elliptical  Orbit  fferaneters' ); 

wltelnC  5.  Azimuth  and  Elevation  Angle'); 

repeat 

writeln;wr1teln;wr1teC  ENTER  QCICE  :'); 

Get  Integer; 

choice  ;*1ntegerdata; 

If  not  (choice  In  [1..5])  then  lnputError:*true; 

if  In  put  Error*  true  then  noise; 
until  InputError*false; 
case  choice  of 


l:Escape  Velocity; 

2:blackbo4y; 

3:circular; 

4:ell iptical ; 

5azel; 

end; 

end; 

2:begin  (*  Aatenna  fein  fettems  *) 
textbackground(  red) ; 

rl rcrr  * 

gotoXY(l.lO); 

writelnC  1.  Ham'); 

writelnC  2.  Helix'); 

writelnC  3.  ffcrabolic  Dish'); 

writelnC  4.  Phased  Arr^y  of  Dipoles'); 

repeat 

writelnjwritelnjwriteC  ENTER  CHOICE  :'); 

Get  Integer; 
choice  :=integerdata ; 

if  not  (choice  in  [1..4])  then  InputError:=true; 
if  Input Error=true  then  noise; 
until  InputError=false; 
case  choice  of 
l:hom; 

2:helix; 

3:dish; 

4:chain(ckmdun) ; 

end; 

end; 

3:begin  (*  Link  Aialysis  *) 
recanpute:=true; 
newinfo  :*true ; 
textbackgnound(red) ; 
clrscr; 
goto  XY( 1,5); 
repeat 

writeC  ENTER  UPLINK  FREQUENCY  (1  -  100  Gtt)  :'); 

Get  Real; 

data.freq.upl  :*realdata; 

if  not((data.freq.upl>=1.0)  and  (data.freq.upl<=100.0)) 
then  Input  Error  :=true; 
if  Input  Error*  true  then  noise; 
until  InputError*false; 
repeat 

write('  ENTER  DOWN. INK  FREQUENCY  (1  -  100  Gfe)  :'); 

Get  Real; 

data.freq.dnl  :=realdata ; 

if  not{(data.fneq.dnl>*1.0)  and  (data.freq.dnl<=100.0)) 
then  Input  Error:  *true; 
if  In  put  Error*  true  then  noise; 
until  InputError=false; 
wrlteln; 

if  havealpha*true  then 


begin 

vriteC  USE  ELEVATION  ANGLES  TOM  MENU  ITEM  #1  ?  ’); 
read!  n(  letter); 

if  ((letter=*Y*)  or  (letter='y‘))  then  skip:=true  else 
skip:=false; 

end; 

if  skip=false  then 
begin 

repeat 

writeC  ENTER  TRANSMITTING  TERMINAL  ELEVATION  ANGLE' + 

'  (de9) 

GetReal; 

data.alpha.upl  :=rea1data; 

if  not((data.alpha.upl>=0)  and  (data. alpha. up! <=360)) 
then  Input  Error  :=true; 
if  InputError=true  then  noise; 
until  InputError=fa1se; 
repeat 

write('  ENTER  TRANSMITTING  TERMINAL  LATITUDE  (deg)  :'); 
Get  Real ; 

data.lat.upl  :=nealdata; 

if  (abs(data.lat.upl)>90)  then  InputError:=true; 
if  Input Error=true  then  noise; 
until  InputError=false; 
writeln; 
repeat 

write('  ENTER  RECEIVING  TERMINAL  ELEVATION  ANGLE  (deg)'  + 
1  :'); 

Getfeal; 

data.alpha.dnl  :=realdata; 

if  not((data.alpha.dnl>=0)  and  (data.alpha.dnl<=360)) 
then  Input  Error  :=t  rue; 
if  InputError=true  then  noise; 
until  InputError=fa1se; 
repeat 

writeC  ENTER  RECEIVING  TERMINAL  LATITUDE  (deg)  :’); 

Get  Real ; 

data.lat.dnl  :=realdata; 
if  (abs(data.lat.<*il)>90)  then  noise; 
until  InputError=false; 

end; 

writeln; 

repeat 

writeC  ENTER  TRANSMITTING  SITE  ALTITUDE  (m)  :'); 

Get  Real; 

data.HD.upl  :=realdata; 

if  (abs(data.tt).upl)>Re)  then  InputError:=true; 
if  Input Error^t rue  then  noise; 
until  InputError*false; 
repeat 

writeC  ENTER  RECEIVING  SITE  ALTITUDE  (m)  :'); 

Get  Real; 

data.HD.dnl  :=realdata; 


if  (abs(data.rt).dil)>te)  then  IhputError:=true; 
until  JhputError=false; 
vriteln; 
getaltitude; 
data.altitude  :=realdata ; 
writeln; 
repeat 

vrite('  ENTER  SYSTEM  BAMUIDTH  (f*iz) 

Get  Real; 

data.bandwidth.up1  :=nealdata; 
data  .bandwidth  .dnl  :=realdata; 

if  not((data.bandwidth.upl>0)  and  (data.bandwidth.upl<1000.0)) 
then  Input  Error  :=true; 
if  InputError=tnie  then  noise; 
until  Input Error=false; 
for  site:=l  to  4  do 
begin 

if  ( (site=l)  or  (site=2))  then  uplink:=true  else 
uplink:=false; 
case  site  of 
1,4  antennal; 

2, 3  antenna  2; 

end; 

wait; 

end; 

upl1nk:=true; 

transmitter; 

uplink:=false; 

transmitter; 

upl  ink:=true; 

receiver; 

uplink  :=false; 

receiver; 

tec  Rower; 

uplink:=true; 

bitrate; 

uplink:«false; 

bitrate; 

end; 

4:begin  (*  Display  Link  Budget  *) 
textbackground(black); 
clrscr; 

textcolor( white) ; 


wr1teln( 'UPLINK' :36,'D0WM.INK':25); 

with  data  do 

begin 

textcolor(lightcyan) ; 

wr1teln( 'Frequency  (GHz)  ’:21,freq.upl  :15:2,freq.dr»l  :25:2); 
wr1teln( 'Slant  tenge  (tan)  ':21,Slantrange.upl :  15:2, 
Slanttenge.dnl  :25:2); 

wr1teln( ' 3-dB  Beamwidth  (deg)  ':21,beamwidth.upl  :15:3, 
beamwidth.dnl  :25:3); 
textcolor(lightblue) ; 
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writeln('Tx  tower  (W)  ':21,tout.up1:15:2.tout.dnl:2S:2); 

write1n('Tx  Circuit  Losses  ':21,TxCircuitloss.upl :15:2, 

TxCircuitLoss.dnl  :25:2); 

wr1teln('Tx  Antema  Gain  '  :21,TxGaln.upl  :15:2, 

TxGain.dnl  :25:2); 

wr1teln('EIRP  '  :21,EIRP.up1 :15:2,EIRP.dm  :25:2); 

textcolor(yellow); 

writeln('Path  Loss  ' :21,tothLoss.upl :15:2, 

Pathloss.dnl  :2S:2); 

write1n('toin  Attenuation  ':21,to1nAtten.upl:15:2, 

RalnAtten.dnl  :25:2) ; 

vr 1  tel n(' Atmos  Attenuation  ' :21,A1rAtten.upl  :15:2, 

AirAtten.dnl  :25:2) ; 
textcolor(llghtgreen) ; 

writeln( 'Received  Useful  tower' :21,fV.upl  :15:2,f¥.dnl  :25:2); 
writeln('to  Antenna  Gain  '  :21,toGain.upl :  15:2, 

lfcGain.dnl :25:2); 

writeln('ffc  Circuit  Losses  ' :21,fcCircuitloss.upl :15:2, 

toCiraiftLoss.dnl  :25:2) ; 

writeln('f&  tolse  Figure  '  :21,toiseFigure.upl  :15:2, 

toiseFIgure.dnl  :25:2); 

writeln('Systan  Eq  Temp  (K)  ':21,Teq.upl  :15:2,Teq.dnl  :25:2); 
textcolor(l  Ightmagenta) ; 

wr1teln('G/T  Fig  of  terit  '  :21,GT.upl  :15:2,GT.dnl  :25:2); 

writeln('C  N  R  '  :21,CNR.upl  :15:2,Cf«.cAnl  :25:2); 

wr1te('Data  tote  (M)/sec)  *:21); 

If  datarate.upbttaxDatarate.upl  then 
textcolor(l  Ightmagentatol  Ink) ; 
write(datarate.upl :15:2) ; 
textcolor(l ightmagenta) ; 
if  datarate.dnl>ffexDatarate.<inl  then 
textcolor(l  Ightmagentatol  Ink) ; 
writeln(datarate.dnl  :25:2) ; 
textcolor(l  Ightmagenta); 

writeln('BER  ':21,'  ',BER.upl:8, 

'  '  BER  dnl  *8) ; 

writeln( 'Available  to/to  ' :21,/Va11toto.upl  :15:2, 

Avail  to  to. dnl  :25:2); 

wr1teln( 'Required  to/to  '  :21,toqtoto.upl  :15:2, 

Reqtoto.dnl  :25:2); 
textcol or (light red); 
wr1te('torg1n  ' :21) ; 

If  margln.upkO  then  textcolor(11ghtre<M>11nk) ; 
write(marg1n.ipl :  15:2) ; 
textcolor(llghtred); 

If  marg1n.dnl<0  then  textcolor(llghtredtollnk); 
writeln(marg1n.dnl  :25:2) ; 
textcolor(l  Ightred) ; 

wr1te(’tox  totarate  (Mo/sec)' :21,Max0etarate.upl  :15:2, 
tox  Data  rate,  dnl  :25:2) ; 

end; 

end; 

5:1f  reca^ute-true  then  (*  Recompute  Link  Budget  *) 


115 


new1nfo:=true; 
textbackground( red) ; 
clrscr; 
gotoXY(l,5); 

wrltelnC  1.  Frequency'); 

writeln('  2.  Antenna  Design'); 

wrltelnC  3.  Transmitter  tower'); 

write1n('  4.  Receiver  Noise  Figure'); 

writeln('  5.  Bit  Error  fete'); 

wrltelnC  6.  Data  Rate'); 

writeln('  7.  Compute  BER  given  data  rate' ) ; 

repeat 

wrlteln  jwriteln  ;write('  ENTER  CK)ICE  :'); 

Get  Integer; 
choice  :=integerdata ; 

if  not  (choice  in  [1..7])  then  InputError:=true; 
if  Input Error=true  then  noise; 
until  InputErrorcfalse; 
clrscr; 

if  cho1ce<>2  then 

begin 

repeat 

vrltelngwrlteC  (U)pllnk  or  (D)ownlink  ?'); 
readln( letter); 

If  not  ((letter  in  ['a'..'z'])  or  (letter  In  ['A'..'Z'])) 
^0rj  noise# 

until  (letter  In  ['a'..'z'])  or  (letter  In  ['A'..'Z']); 

If  ((letter='u')  or  (letter='U')j  then  uplink:=true 

else  uplink:=false; 

end; 

case  choice  of 

l:beg1n  (*  enter  new  frequency  *) 
repeat 

gotoXY(2,5); 

wr1te('  ENTER  NEW  FREQUENCY  (GHz)  :'); 

Get  Real; 
junk:*realdata; 

If  not((junk>*1.0)  and  (junk<= 100.0))  then 
Input  Error  :=t  rue; 

If  InputError*true  then  noise; 
until  InputError*false; 

If  upl1nk»true  then  data.freq.upl  :»junk  else 
data.freq.tfril  :*junk; 

end; 

2:beg1n  (*  enter  new  antenna  design  *) 
repeat 

gotoXY(l,5); 

»r1teln('  1.  Upl  Ink  Trananitter' ) ; 

wrltelnC  2.  Uplink  Receiver'); 

wrltelnC  3.  fowl  Ink  Transmitter'); 

wrltelnC  4.  fowl  Ink  Receiver'); 

wrlteln  jwrlteln prite( '  ENTER  ANTEW1A  SITE  :'); 


Get  Integer; 
site:=integendata; 

If  not(s1te  In  [1..4])  then  InputError  :*true; 

If  InputError=true  then  noise; 
until  IhputError=fa1se; 

If  site<=2  then  upl1nk:*true  else  upl1nk:=false; 

If  ((s1te*l)  or  (s1te=4))  then  antennal  else  antenna2; 
end; 

3:beg1n  (*  enter  new  transmitter  power  *) 
repeat 

gotoXY(2,5); 

write('  ENTER  OUTPUT  POWER  LEVEL  (W)  :'); 

Get  Real; 
junk:=realdata; 

if  junk<=0  then  InputError  :=true; 

If  Input Error=true  then  noise; 
until  InputError=false; 
if  upl  1nk=true  then  data. tout .upl  :=junk  else 
data.tout.dnl  :-junk; 

end; 

4:beg1n  (*  enter  new  receiver  noise  figure  *) 
repeat 

gotoXY(2,5); 

write{'  ENTER  NEW  NOISE  FIGURE  :'); 

Get  Real; 
junk  :*real  data; 

if  not((junk>0)  and  (junk<*1.0))  then 
Input  Error  :=true; 
if  InputError=true  then  noise; 
until  InputErrop=false; 

if  uplink=true  then  data,  ftoise  Figure  .upl  :=junk  else 
data.NolseFigure.fiM  :=junk ; 

end; 

5:begin  (*  enter  new  bit  error  rate  *) 
standard  :=true; 
repeat 

gotoXY(2,5); 

writeC  ENTER  fCW  DESIRED  BIT  ERROR  RATE  :'); 

Get  Real; 
junk  :=real  data; 

if  not((junk>0)  and  (junk<1.0))  then 
InputError  :=true; 
until  Input  Error*  false; 
if  upl1nk*true  then  data.BER.upl  :*junk  else 
data.BER.dnl  :=junk; 

end; 

6:begin  (*  enter  new  data  rate  *) 
standard  :*true; 
repeat 

gotoXY(2,5); 

writeC  ENTER  DATA  RATE  DESIRED  (to/sec)  :'); 

Get  Real; 

junk  :*real  data; 
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If  junk<le-3  then  InputError:=true; 
if  InputError=true  then  noise; 
until  Input Error*false; 

if  uplink=true  then  data.datarate.upl  :=junk  else 
data.datarate.dnl  :=jtnk ; 
end; 

7:begin  (*  compute  BER  given  data  rate  *) 
standard  :=false; 
repeat 

gotoXY(2,5); 

write('  ENTER  MAX  DATA  RATE  DESIRED  (H>/sec)  ;') 

Get  Real; 

rate  :=real  data; 

if  rate<le-3  then  InputError  :=true; 
if  Input  Error=true  then  noise; 
until  IhputError=false; 
end; 

end; 

with  data  do 

begin  (*  perform  actual  recomputations  *) 
case  trunc(AntenraType.upl )  of 
1 '.begin 

TxG&in.upl  :=10*1ogl0(gainl(argl.upl  .freq.upl)); 
beanwidth.upl  :=bwl(argl.upl  .freq.upl ) ; 
end; 

2:begin 

Tx&in.upl  ;»10*logl0(gain2(argl.upl  ,arg2.upl , 
freq.upl)); 

beanwidth.upl  :=bw2(argl.upl  ,arg2.upl  .freq.upl ) ; 

end; 

3:beg1n 

Tx&1n.upl  :=10*logl0(gain3(argl.upl  .freq.upl)); 
beanwidth.upl  :=bw3(argl.upl  .freq.upl ) ; 
end; 

4:begin 

TxG&in.upl  :=10*1ogl0(gain4(arg3.upl  .argl.upl , 
freq.upl)); 

beanwidth.upl  :«bw4(arg3.upl  .argl.upl  .freq.upl); 
end; 

end; 

case  trunc(  Antenna  Type  .Chi)  of 

l:fc&1n.chl  :«10*logl0(gainl(argl.dnl  .freq.dnl )) ; 
2:faGBln.dnl  :»10*logl0(gain2(argl.<lnl  ,arg2.dnl , 
freq.dnl)); 

3:RxGfcin.dnl  :«10*loglQ(gain3(argl.dnl  .freq.dnl)) ; 
4:ftGa1n.dnl  :*10Mogl0(ga1n4{arg3.dnl  .argl.dnl , 
freq.dnl)); 

end; 

feGfein.upl  :-10Mogl0(ga1n3(arg5.upl .freq.upl )) ; 
lx€toln.(*nl  :«IO*1oglO(ga1n3{arg5.<iil  .freq.dnl )) ; 

EIRP.upl  :-i$otropfc(Rout.upl  .TxGaln  jjpl  .TxCircuitloss.upl ) ; 
EIRP.Chl  :«1sotrop1c(R)ut.dnl  .TxGain.dll  .TxCIrcuitloss.chl ); 
leq.upl  :«eqta*>(faAntTemp.upl  .ffciseFigure.upl , 


faCIrcuitLoss.upl); 

Teq.dnl  :«eqtanp(RxMTerp.dnl  .ftlseFIgure.dnl , 
ftcCIrcurtLoss.dil); 

GT.upl  :=figner1t(lfcG&1n.upl  .Teq.upl ) ; 

GT.dil  :*f1gner1t(lfcGa1n.dnl  .Teq.dnl ); 

ftthLoss.upl  :-"36.6-20*logl0($1antrange.upl*le3*6.214e-4* 

freq.upl*le3); 

f^thloss  .dnl  :«-36. 6-20*logl0(sl  ant  range  .dnl*le3*6. 214e-4* 
freq.dnl*le3); 

AlrAtteruupl  :»a1r(freq.upl  .alpha.upl ) ; 

AlrAtten.dnl  :=air( freq.ckil  .alpha  .dnl ); 

RainAtten.upl  :=rain(ra1nnate.upl  .freq.upl  .alpha .upl  .Hb.upl , 
lat.upl); 

teinAtten.dnl  :=rain(rainrate.dnl  .fneq.dnl  .alpha.dnl  .ft  .dnl , 
lat.dnl); 

FV-.upl  :=power(Ftut.upl  .TxCircuitLoss.upl  .Txfein.upl , 
TxEfficiency.upl  .PathLoss.upl  .AlrAtten.upl  .RainAtten.upl , 
FfcEfficiency.upl  ,lfc&1n.upl  .fccCIrcuitLoss.upl ) ; 

Pr.dnl  :=power(  tout  .dnl  .TxCIrcuft Loss  .dnl  .TxGaln.dnl . 
TxEfficiency.dnl  .fteth  Loss  .dnl  ,A1rAtten.dnl  .telnAtten.dnl . 
^Efficiency  .dnl  .ffcGaln.dnl  .faC1rcuitass.dll); 

CffUupl  :=cnratio(Pr.upl  .Teq.upl ) ; 

CNR.dnl  :=>cnrat1o{ft\dn1  .Teq.dnl); 

If  standard=true  then 
begin 

/Wail  afb.upl  :*aft/Wail  (Pr  jipl  .Teq.upl  .datarate.upl ) ; 
/Wall  aft  .dnl  : -aft/Wall  (fV.dnl  .Teq.dnl  .datarate.dnl ) ; 
fteqEbto.upl  :=atofeq(trunc(modtype.up1 )  /n.upl  ,BER.upl ) 
Reqaft.dnl  :=Ebffcteq(trunc(modtype.dnl)  ^n.dnl  .BER.dnl ) 
ftexDBtaFfete.upl  :4fexfete(Pr.upl  .Teq.upl  ,Reqaft.upl ) ; 
capacity(bandw1dth.upl  .Pr.upl  .Teq.upl  .ffex  Datarate.upl . 
M.upl ) ; 

if  ((shannon  flag-true)  or  (nyqulstflag-true))  then 
ffexDatarate.upl  :*1  Imlt ; 

textfetatete.dnl :=flax  fete  (Pr.dnl  .Teq.dnl  .Reqaft.dnl ) ; 
capacity(bandw1dth.dnl  .Pr.dnl  .Teq.dnl  J%x  Datarate.dnl , 
H-dnl); 

If  ( (shannon flag-true)  or  (nyquistflag-true))  then 
ffexDatarate.dnl  :«1  Imlt ; 
margin.upl  : -Avail  aft  .upl-Reqaft.upl ; 
margin.dnl  : -/Wail  aft  .dnl-feqafb.dnl; 
end  else 
begin 

If  upl ink-true  then 
begin 

ait  :»CWT.upl-10*logl0(  rate*  le6) ; 
err  :-b1terror(trunc(modtype.upl )  .m.upl , 

end  else 
begin 

aft  :<f*.dnl-10*logl0(rate*le6) ; 
err:-b1terror(trunc(inodtype.dil )  .m.dnl , 

aft); 


end; 

clrscr; 

gotoXY(l,5); 

if  uplink=true  then  *r1teln(title[9]:46)  else 

writeln(title[10]:46) ; 

gotoXY(l.lO); 

textcol  or(  1  i  ghbnagenta ) ; 

writelnC  Desired  Data  tote  =  ',rate:0:2,'  M>its/sec'); 
writelnC  /V ailable  Eb/No  *  \B>lto:0:2,'  dB'); 
vriteln('  Bit  Error  fete  *  *,err:8); 

end; 

end; 
end  else 
begin 

noise; 

textcolor(>ellow); 

writelnjwritelnC  ***  UNABLE  TO  RECCMRJTE  ***'); 
textcolor( white) ; 

6nd  * 

6:begin  (*  Store  Link  File  *) 
new1nfo:*false; 
textbac kg round ( red) ; 
clrscr; 
gotoXY(l,5); 
saveit; 

6nd  * 

7:begin  (*  tecall  Link  File  *) 
recompute  :=true; 
newinfo:*false; 
textbackground(  red) ; 
clrscr; 
gotoXY(l,5); 
repeat 

writeC  ENTER  FILENAME  (Vo  extension)  :'); 
neadln(name); 

if  not((length(name)>*l)  and  (length(name)<*8))  then  noise; 
until  (1ength(name)>*l)  and  (length(name)<»8); 
name  :»path+name*'. Ink' ; 
if  exist(name)*true  then 
begin 

assign(11nkfile,name) ; 
reset(linkflle); 

with  data  do  nead(11nkfile,data); 
close(linkflle); 
end  else 
begin 

noise; 

textcolor(yellow4bl  ink) ; 
writeln  ;writeln('  ***  FILE  NdCXISTENT  ***'); 
textcolor(trfiite) ; 
end; 

end  * 

8:beg1n  (*  Delete  Link  File  *) 


// 


textbackground(red) ; 
clrscr; 
gotoXY(l.lO); 
repeat 

wr1te('  DELETE  WHICH  FILE  ?'); 
readln( kill none); 

if  not({1ength(killname)>=l)  and  (length(killname)<=8)) 
then  noise; 

until  (length(kil1name)>=l)  and  (length(killname)<=8); 
writeln;write('  ARE  YOU  SURE  ?  (YorN)'); 
readl  n(  letter) ; 

if  ((letter*'Y')  or  (lette r='y' ))  then 
begin 

name  :*patMk  ill  name*'.  Ink' ; 
if  exist(name)=true  then 
begin 

clrscr; 
gotoXY(l,  12); 
textcolor(yel  IowHjI  ink) ; 
vriteln('D  E  L  E  T  I  N  G  FILE* :52); 
assign(linkfile.rame); 
erase(linkfile);  (*  erase  disk  file  *) 
reset(library); 
revrite(  buffer) ; 
while  not  EDF( library)  do 
begin  (*  copy  all  other  filenames  to  buffer  *) 
read(  library  .name); 

If  pos(  kill  name  .name)  *0  then  write(buffer,name) ; 

end; 

reset(  buffer) ; 
revrite(llbrary); 

while  not  EDF(buffer)  do  (*  copy  buffer  into  library  *) 
begin 

read(buffer,name); 

write(library.name); 

end; 

close(buffer) ; 
close(library); 
clrscr; 
end  else 
begin 

textcolor(yello#+bl  ink) ; 
noise; 

writeln^ritelnC  ***  FILE  NONEXISTENT  *•*•); 
textcolor(»rfi1te) ; 

end; 

end; 


9:begin  (*  List  Link  File  Directory  *) 
reset(l ibrary) ; 
textbackground(red) ; 
clrscr; 

if  f11es1ze(library)«0  then 


»v«S 


v.v; 
v  •>' 


W 


begin 

clrscr; 

noise; 

textcolor(jellovHbl  Ink) ; 
witeln;»r1teln('  ***  LIBRARY  9*>TY  ***'); 

end  * 

vhile  not  B0F( library)  do 
begin 

read(  library, name); 
writelnC  '.nane); 

end; 

close( library); 

end; 

10:begin  (*  Change  Default  Drive  *) 
textbackground(ned); 
clrscr; 
gotoXY(l,5); 

vr1te1n('N0TE:  the  path  disk^data  mist  exist  on  the  DATA  disk  !'); 
repeat 

writeln prite( '  ENTER  KW  DEFAULT  DRIVE  :'); 
readln(letter) ; 

if  not((letter  in  ['a' ..'e'])  or  (letter  in  ['A'..'E'])) 
then  noise; 

until  (letter  in  ['a'./e'])  or  (letter  in  ['A'..'E']); 
path:=letterf'  :data»' ; 
ass1gn(l 1brary,path+' 1 Ink.dlr' ) ; 
assign(buffer  ,path+'  temp.sto' ) ; 
if  not  ex1st(path+'11nk.dir')  then 
begin 

revr1te(  library); 
close(llbrary); 

end; 

end; 

end;  (*  end  of  case  *) 
if  optionOIl  then  wait; 
until  optional;  (*  EXIT  routine  *) 
if  newinfO“true  then  (*  last  chance  to  save  data  *) 
begin 

writeln g»r1te( '  WANING:  DATA  NOT  SAVED  ...  save  ?  (Y,N)  ’); 
read (letter)  | 

If  ((letter*'Y')  or  (letter*'y'))  then  savelt; 

end; 

textcolor(vhite) ; 
textbackground(black) ; 
clrscr; 
end. 


(*  custom  graphics  driver  *) 
procedure  custcmp1ot(x  »y  jn  :integer) ; 

(*  plot  a  point  *) 
begin 

plot(x„y,n) ; 
end; 

(**) 

procedure  custandraw(xl,yl»x2»y2,n  :integer); 

(*  draw  a  line  between  two  points  *) 
begin 

draw(xl*yl,x2,y2,n); 

end; 

(**) 

(*  detennine  if  a  disk  file  exists  *) 
function  exist(name:filenare):boolean; 
var  fil(*  duimy  filename  *):file; 
begin 

assign(fil  .name) ; 

§$I-t 

reset(fil); 

§$I+t 

exist  :=(I0result=0); 
end; 

(**) 

(*  poll  user  for  satellite  altitude  *) 

procedure  getaltitude; 

begin 

repeat 

wr1te('  ENTER  SATELLITE  ALTITUDE  (m)  (  or  6  for  geosynchronous  )  :'); 
Get Real ; 

if  realdata<1.5e5  then  InputEiTor:=true; 
if  Input Error=t rue  then  noise; 
until  InputError=false; 
end; 

(**) 

(*  store  an  integer  input  *) 
procedure  Get  Integer; 
var  i  (*  loop  variable  *), 

code  (*  error  code  *) integer; 
ctaracter  (*  individual  input  character  *):char; 
datain  (*  input  data  string  *):string[5]; 
begin 

InputError:*false; 

readln(datain); 

for  1  :*1  to  length(datain)  do 

begin 

character  :=coRy(datain  ,i  ,1) ; 

If  not  (character  in  legal  Integer)  then  Input  Error  :=true  else 
val (datain .Integerdata  .code) ; 

end; 

end; 

(**) 

(*  make  beeping  sound  *) 


procedure  noise; 
begin 

sound(lOOO); 
delay(500) ; 
nosound; 
end; 

(**) 

(*  stone  a  real  input  *) 

procedure  Get Real ; 

var  i  (*  loop  variable  *), 

code  (*  error  code  *) integer; 
character  (*  individual  input  ctaracter  *):char; 
datain  (*  input  data  string  *)  :string[15]; 
begin 

lnputError:=false; 

readln(datain); 

if  ((copy(datain,l,l)='g')  or  (copy(datain,l,l)='G'))  then  realdata:=3.5784e7 

else  for  i :=1  to  length( datain)  do 

begin 

character  :=ooRy(datain  ,i  ,1) ; 

if  not(character  in  legalreal)  then  Input  Error  :=true  else 
val (datain ,nea1data .code) ; 

end; 

end; 

(**) 

(*  save  a  file  on  disk  *) 
procedure  saveit; 
begin 
repeat 

writeln *rite( 1  ENTER  FILENAME  (  1  -  8  chars  )  :'); 
neadln(rane) ; 

if  not((length(name)>=l)  and  (length(name)<=8))  then  noise; 
until  (length(name)>=l)  and  (length(nane)<=8); 
name :=path+rBme+'. Ink' ; 
dwmyname  :=name ; 

del  ete(  duimyname ,  1 ,1  ength  ( path )) ; 
neset(l ibrary) ; 

(*  see  if  file  already  in  library  *) 

while  not  EDF( library)  do 

begin 

read(  1  ibrary  ,ol  dname) ; 

if  pos(oldiame,duirni>name)00  then  Alreadylhere:=true; 

end; 

(*  if  not  already  in  library,  write  to  disk  and  update  library  *) 

if  Already  The  re=  false  then 

begin 

asslgn(llnkflle^ame) ; 

newrite(llnkflle); 

with  data  do  write(11nkfile ,/data); 

close(linkflle); 

reset  (library); 

If  filesize( library)  >  Othen  seek(l  1brary.fi lesize( library)); 
wr  1  te(  1 1  brary  ^lumnyrane ) ; 


c1ose(l ibrary) ; 
end  else 
begin 

A1  ready  The  re  :=fia1se; 
noise; 

textco1or(yellow+bl  ink) ; 

writelnjwrite1n('  ***  FILE  ALREADY  EXISTS  ***'); 
textcoloK  white) ; 

end; 

end; 

(**) 

(*  suspend  execution  pending  keyboard  entry  *) 
procedure  wait; 
label  here; 

var  wide(*  width  of  message  *)  integer; 
begin 

gotoXY(l,25); 
textcolor(yel  low) ; 

if  ((option=2)  and  ((choice=l)  or  (choice=3)))  then  wide:=34  else  wide:=54; 
write( '[press  any  key  to  continue]' .-wide); 
textcolor(white) ; 

hererif  keypressed=false  then  goto  here; 
end; 

(**) 

(*  convert  dB  to  nunbers  *) 
function  alog(d) -.real) -.real; 
begin 

alog:=pwr(  10.db/10) ; 
end; 

(**) 

(*  first-order  bessel  functions  of  the  first  kind  *) 
function  bessel (x: real): real; 

var  ctr  (*  counter  to  truncate  insignificant  series  terms  *), 

1  im  (*  max  #  terms  *) , 
sign  ,i  ,1  ,z:integer; 

a, subterm  .term  .old  .-real ;  (*  constants  and  loop  variables  *) 
begin 
ctr:=0; 

1  im  :=200; 
old:=0; 
a:=x/2; 
term:=l; 
sign :=1 ; 
i  :=0; 

(*  evaluate  series  until  200  terms  added  or  3  consecutive  terms  not 

contribute  to  sun  *) 

while  ( ( i<l im)  and  (ctr<=2))  do 

begin 

i  :=i+2; 
sign:=-s1gn; 
subterm  :=1; 

for  1  :=1  to  1  dlv  2  do  subterm  :=$ubternv'((  1+1  )*1); 
for  z:*l  to  1  do  subterm ^subternfx ; 


term  :=term*{  subtend  pwr(2,i )  )*sign ; 
if  a*teraw>ld  then  ctr:=ctr+l; 
o1d:=a*term; 

end; 

bessel  :=old; 
end; 

(**) 

(*  natural  logarithms  *) 
function  log2(x:real):real; 
begin 

log2:=ln(x)/ln(2. 0) ; 
end; 

(**) 

(*  comnon  logarithms  *) 
function  loglO(x:real )  :real ; 
begin 

loglO:=ln(x)/ln( 10) ; 
end; 

(**) 

(*  raise  x  to  the  power  of  n  *) 
function  pwr(x  />: real ): real ; 
begin 

if  x=0  then  pwr:=0  else  if  n=0  then  pwr:=l  else 
if  x<0  then 

if  (trunc(n)  mod  2)=0  then  p*r  :=exp(n*ln(-x))  else  p*r:=-exp(n*ln(-x)) 
else  pwr:=exp(n*ln(x}); 

end; 

(**) 

(*  ftarcun's  Q-function  *) 

function  q(u : real ): real ; 

var  1  (*  lower  limit  of  integration  *), 

N  (*  #  partitions  *), 

inc  (*  partition  width  *), 

sun, loop  (*  loop  variables  *), 

area  (*  intermediate  result  *), 

erf  (*  value  of  corresponding  error  function  *):real; 

(**) 

(*  approximation  of  Q-function  *) 
function  qapprox(x  :real )  :real ; 
begin 

if  x>13  then  qapprox  :=0  else 

gapprox  :=exp(-(sqr(x)/2))/(x*sqrt{2*pi )) ; 

end; 

(**) 
begin 
clrscr ; 

textcolor()el  low+bl  ink) ; 
gotoXY(l,12); 

writeln('S  T  A  N  D  B  Y* :46) ; 

(*  evaluate  error  function  using  Simpson's  rule  *) 
if  iKA.24  then 


u:=u/sqrt{2); 

N:=1000.0; 

inc:=(u-l)/N; 

(*  get  first  and  last  points  *) 
sum:=exp(-sqr(1  ))+exp(-sqr(u) ) ; 
loop:=l+inc; 

(*  get  remaining  points  *) 
while  loop<=1+(N-l)*inc  do 
begin 

sun  :=sun+4*exp(-sqr(  loop) ) ; 
loop:=loop+2*inc; 

end; 

loop:=l+2*inc; 

while  loop<=l+(N-2)*inc  do 

begin 

sun :  =suiH-2*exp( -sqr  ( 1  oop) } ; 
loop:=loop+2*inc; 

end; 

area  :=sun*inc/3.0; 
erf  :=2/sqrt(pi  )*area ; 

(*  extract  Q-ftinction  from  corresponding  error  function  *) 
q:=(l-erf)/2.0; 

end  else  q:=qapprox(u) ;  (*  if  u>=4.24  then  use  approximation  *) 

textcolor(v^iite); 

clrscr; 

end; 

(**) 

(*  inverse  Q-function  *) 

function  qinverse(q:real ) :real ; 

var  xl,x2,x3  (*  duimy  variables  *):real; 

(**) 

(*  equation  for  argunent  given  the  corresponding  Q-function  *) 

function  f(x  :real ) :real ; 

begin 

f : =sqr(  x ) +2*1  n  ( x ) +2*1  n  ( cf  sqrt  ( 2*pi ) ) ; 

end; 

(**) 

(*  derivative  of  above  equation  *) 

function  g(x:real):real; 

begin 

g  :=2*x+2/x; 

end; 

(**) 

begin 

clrscr; 

textcolor(>el  low+bl  ink) ; 
goto  XY(  1,12); 

writeln('S  T  A  N  0  B  Y':46); 

(*  establish  search  range  [xl*x2]  *) 

xl  :=0»  1 ; 

x2:=13; 

(*  use  method  of  halving  the  interval  to  localize  root  *) 
while  Q.5*abs{xl-x2)>=le-6  do 


begin 

x3:=(xl*x2)/2.0; 

if  abs(f(x3))  +  abs(f(xl))  O  abs(f(x3)+f(xl))  then  x2:=x3  else  xl:=x3; 
end; 
x2:=x3; 

(*  use  ftewton's  method  to  precisely  isolate  root  *) 

while  abs(x2-xl)>=le-6  do 

begin 

xlisx2* 

x2:=*l-f(xl)/g(xl); 

end; 

qinverse:=x2; 
textcolor(vrtiite) ; 
clrscr; 
end; 

(**) 

(*  inverse  cosine  *) 
function  arccos(x :real ) :real ; 
begin 

if  x=0  then  arccos:=pi/2  else 
arccos:=arctan(sqrt(l-sqr(x))/x) ; 

end; 

(**) 

(*  inverse  sine  *) 

function  aixsin(x  :real )  :real ; 

begin 

if  x=l  then  arcsin:=pi/2  else 
arcsin  :=arctan(x/sqrt(l-sqr(x) ) ) ; 

end; 

(**) 

(*  tangent  *) 

function  tan(x  :real ) :real ; 
begin 

tan:=sin(x)/cos(x); 

end; 

(**) 

(*  plot  gain  pattern  for  horn  *) 
procedure  hom; 

var  flag  (*  plot  first  point  flag  *):boolean; 
fx.fy  (*  scaling  factors  *), 
d  (*  hom  dimension  *), 
f  (*  frequency  *) , 
inc  (*  plotting  resolution  *), 
theta  (*  polar  angle  *), 
x^ain^nax  (*  gain  variables  *):real; 
oldx ,oldy .xplot .yplot  (*  plotting  coordinates  *) -.integer; 
begin 

flag:=false; 

fy:«100;fx:«220; 

textbac kg round (blue) ; 

clrscr; 

repeat 

writeC  ENTER  FREQUENCY  (1  -  100  (Hz)  :'); 


Get  Real; 
f:=realdata; 

if  not((f>=1.0)  and  (f<=100.0))  then  InputError :=true; 
if  Input  Error=true  then  noise; 
until  InputError=false; 
repeat 

write('  ENTER  HORN  DIMENSION  (m) 

Get  Real; 
d:=realdata; 

if  d<=0  then  InputError  :=true ; 
if  Input Error^true  then  noise; 
until  InputEmor=fa1se; 
max :=10*logl0(gainl(d ,f ) ) ; 
graphcolonnode ; 
palette(3); 

graphbackground(black) ; 
theta  :=-pi/2; 
inc:=0.C6; 

(*  sweep  theta  180  degrees  to  generate  polar  plot  *) 

while  theta<pi/2  do 

begin 

x  :=pi*sin(theta) ; 
gain:=$qr(sin(x)/x); 
xplot  :=50mxind(gain*cos(theta)*fx) ; 
yplot  :=100-round(gairr*sin(theta)*fy) ; 

(*  plot  first  point,  then  draw  lines  to  each  successive  point  *) 

if  flag=false  then 

begin 

flag:=true; 

customplot(xplot  ,yplot  ,2) ; 
oldx  :=xplot  ;oldy  :=yplot ; 
end  else 
begin 

customdraw(oldx  ,olc|y,xplot  ,yp1ot  ,2) ; 
oldx  :=xplot  yjldy  :=yplot ; 

end; 

theta  :=theta+inc ; 
end; 

(*  draw  little  horn  to  represent  antenna  *) 
custandraw(  47,94,47,106,3) ; 
customdraw(47,106,40, 103,3) ; 
custandraw(40,103,36,103,3) ; 
custandraw( 36,103,36,97, 3) ; 
custandraw(  36,97,40,97, 3) ; 
customdraw(  40,97, 47, 94,3) ; 
gotoXY(l.l); 

writeln('FR£Q  =  ' ,f :0:1,'  GHz’); 
writeln('K)RN  DIMENSION  =  '^:0:1,’  m'); 
gotoXY(l,20); 

writeln('M*X  GAIN  =  ‘/naxiO:!,*  dB'); 

wr1teln;writeln(’ FIRST  SIDQjOBE  LEVa  «  '  ^nax- 13. 0:0:1,’  dB'); 

end; 


(*  plot  gain  pattern  for  helix  *) 
procedure  helix; 

var  flag  (*  plot  first  point  flag  *), 

check  (*  flag  to  supress  gartage  on  screen  *):boolean; 

fx,fy  (*  scaling  factors  *), 

max  (*  maxinun  gain  *) , 

f  (*  frequency  *) , 

k  (*  wave  number  *) , 

radius  (*  helix  radius  *), 

alpha  (*  pitch  angle  *), 

turns  (*  #  turns  *), 

d  (*  height  of  one  coil  *), 

theta  (*  polar  angle  *), 

inc  (*  plotting  resolution  *), 

AF,PF  (*  array  factor,  pattern  factor  *), 

Hidge(*  plotting  error  adjustment  *):real ; 
oldx.oldy.xplot^yplot  (*  plotting  coordinates  *):integer; 
begin 

textbac kg round (blue) ; 

clrscr; 

fudge  :=0; 

check  :=false; 

flag:=false; 

fy:=500; 

repeat 

writeC  ENTER  FREQUENCY  (.1  -  1.0  GHz)  :'); 

Get  Real; 

f  :=realdata; 

if  not  ((f<=1.0)  and  (f>=0.1))  then  InputError:=true; 

if  Input Error=t rue  then  noise; 
until  InputError=false; 
k:=2*pi*f*le9/3e8; 
repeat 

writeC  ENTER  HELIX  RADIUS  (optimun  =  M/k:0:3,'  m)  :'); 

Get  Real; 

radius  :*real  data; 

if  not  ( ( rad1us>=0. 9/k)  and  (radius<=l.l/k))  then  InputError  :=true; 

if  Input Error=t rue  then  noise; 
until  InputError=false; 
repeat 

writeC  ENTER  PITCH  ANGLE  (12  -  18  deg,  14  =  optimun)  :'); 

Get  Real ; 

alpha  :=real  data; 

if  not  ((alpha>«12)  and  (a1pha<=18))  then  InputError;=true; 

If  Input Error*t rue  then  noise; 
until  InputError*false; 
alpha  :*alpha*dr; 
repeat 

writeC  ENTER  fOBER  OF  TURNS  (N  >  3)  :'); 

Get  Real; 

turns  :»1nt(rea1data); 

if  tums<4,0  then  Input  Error  ;*  true; 

if  InputError*true  then  noise; 
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until  Input Error=false; 
d : =2*pi *  radi us*tan (alpha) ; 
max  :=ga1n2(2*radius  ,tums*d  ,f) ; 
fx  :=fy*2. 2+tums*40. 0; 
hires; 

hirescolor(lightblue) ; 
customdraw( 20,180,620, 180,1) ; 
custandraw(  320, 0,320, 180, 1) ; 

1nc:=l/(6.0*tums);  (*  adjust  resolution  *) 

oldx:=-500; 

theta :=pi +0.1; 

(*  sweep  theta  *) 
while  theta<2*pi  do 
begin 

if  theta< =2*pi - fudge+0. 01  then 

begin  (*  total  pattem=array  factor  X  pattern  factor  *) 

AF:=abs(sin(cos(theta)*tums*pi*d)/(tums*sin(cos(theta)*pi*d))); 
PF:=cos(k*radius*sin(theta)-0.75*pi)*(1.0/(pi*k*radius*sin(  theta))); 
xplot  :=320mxjnd(AF*PF*cos(theta)*fx) ; 
yplot  :=180+round(AF*PF*sin(  theta  )*fy) ; 

(*  suspend  plotting  until  first  change  in  direction  *) 

if  check=false  then  if  xplot>oldx  then  oldx:=xplot  else  if  xplot>300  then 

begin 

check  :=true; 
fudge :=theta-pi ; 

end; 

(*  plot  first  point,  then  draw  lines  to  each  successive  point  *) 
if  check=true  then 
if  flag=true  then 
begin 

customdraw(ol  dx  ,ol  dy^plot  »yplot,l); 
ol  dx :  =xpl  ot  »1  C(y :  =ypl  ot ; 
end  else 
begin 

flag:-true; 

customplot(xplot .yplot ,1) ; 
ol  dx :  =xpl  ot  ;ol  dy  :=ypl  ot ; 

end; 

end; 

theta  :=theta+inc; 
end; 

gotoXY(l.l); 

wr1teln('FREQ  =  * ,f :0:1,'  GHz'); 
writelnC RADIUS  =  '  .radius :0:2,'  m'); 
writelnC  PITCH  =  '  ,alpha:0:l,'  deg’); 
writelnC#  TURNS  =  ',tums:0:0); 
wr1teln(' MW  GAIN  =  'jiax:0:l,'  dB'); 
end; 

(**) 

(*  plot  gain  pattern  of  parabolic  reflector  *) 
procedure  dish; 

var  xl»yltfx,fxl,fy,fyl  (*  scaling  factors  *), 
f  (*  frequency  *) , 
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diameter  (*  dish  diameter  *), 
dl  (*  diameter/wavelength  *), 
x  .gain  jm.  (*  gain  variables  *), 
theta  (*  polar  angle  *), 
inc  (*  plotting  resolution  *):real; 
oldx,ol(ly,xp1ot1yplot  (*  plotting  coordinates  *) integer; 
begin 

fy  :=100. 0;fx  :=130. 0; 

fyl:=20i0;fxl:=10.0; 

textbackground(blue) ; 

clrscr; 

repeat 

writeC  ENTER  FREQUENCY  (1  -  100  GHz)  :'); 

Get  Real ; 
f:=nealdata; 

if  not  ((f>=1.0)  and  (f<=100.0))  then  IhputError :=true; 
if  Input  Error^t  rue  then  noise; 
until  InputError=false; 
repeat 

vriteC  ENTER  DISH  DIAMETER  (m) 

Get  Real ; 

diameter  :=real  data; 
if  diameteK=0  then  InputError  :=true ; 
if  Input  Error=true  then  noise; 
until  InputError=false; 
dl  :=diameter*f*le9/3e8; 
max  :=10*logl0(gai  n3(di  ameter  ,f ) ) ; 
graphcolormode; 
palette(3); 

custandraw(  160, 100, 300, 100, 1 ) ; 

customdraw(  160,40,160,160,1) ; 

oldx:=160;oldy:=100; 

inc:=0.1; 

xl:=0; 

(*  draw  small  parabola  to  represent  antenna  *) 

while  xkl.O  do 

begin 

yl:=sqrt(xl); 

xplot  :=16CHround(xl*fxl) ; 

yplot  :=100-round(yl*fyl) ; 

customdraw(  ol  dx  ,ol  dy  ,xpl  ot  ,ypl  ot ,  1 ) ; 

oldx  :=xplot  ;olc|y  :=yplot ; 

xl:=xl+inc; 

end; 

xl:=0; 

oldx:=160y)ldy:=100; 
while  xKl.O  do 
begin 

yl:=sqrt(xl); 
xplot  :*160mxmd(  xl*fxl ) ; 
yplot  :*100mxind{yl*iyl) ; 
custamdraw(ol  dx  ,ol  dy  .xplot  .yplot ,  1) ; 
oldc  :«xplot  ;ol4y  :«yp1ot ; 
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xl:*xl+1nc; 

end; 

oldx  :*160jolcty  :*100; 
theta ;=-pi/2; 

(*  sweep  theta  over  180  degrees  and  plot  gain  pattern  *) 

while  theta<pi/2  do 

begin 

x :=pi*sin(theta) ; 
gain :=sqr(2*bessel (x)/x ) ; 
xplot:=160m>und(gain*cos(theta)*fx) ; 
yplot  :=100-round(gain*sin(theta)*fy) ; 
custandraw(ol  dx  .oldy  .xplot  .yplot  ,2) ; 
oldx  :=xplot  y)ldy:=yplot ; 
theta  :=theta+inc ; 

end; 

gotoXY(l,l); 

writeln('FREQ  =  ',f:0:l,'  GHz'); 
writeln(' DIAMETER  =  '  .diameter :0:1,'  m'); 
gotoXY(l,20); 

writeln('MAX  GAIN  =  ’.maxiO:!,*  dB'); 

writeln ;writeln( 'FIRST  SIDELOBE  LEVEL  =  '  jnax-17.6:0:l,'  dB'); 

end; 

(**) 

(*  plot  gain  pattern  of  dipole  array  *) 
program  Phased  Array; 
type  lega1=set  of  char; 

var  InputError  (*  illegal  nuneric  input  *):boolean; 
integerdata  (*  integer  input  data  *) : integer; 
realdata  (*  real  input  data  *):real; 
legal  integer  .legal  real  (*  valid  input  sets  *):legal; 
inc  (*  plotting  resolution  *), 
fy.fx  (*  scaling  factors  *), 
m,mag  (*  plot  magnification  variables  *), 
theta  .phi  (*  spherical  coordinate  angles  *), 
leng  {*  length  of  elements  *) , 
k  (*  particular  direction  cosine  *), 

AF.AF1.AF2  (*  array  factors  *) , 

PF  (*  pattern  factor  *) , 

EFtheta.EFphi  (*  element  factors  *), 

Gtheta,Q>hi  (*  gain  for  theta  and  phi  components  *):real; 

n  ,i  j  (*  loop  variables  *) , 

cd  (*  vector  current  direction  *), 

dist  (*  current  distribution  *) , 

x,y  (*  #  subarrays  in  x  and  y  directions  *), 

xplot»yplot,oldx,oldy  (*  plotting  coordinates  *) integer; 

xnun,ynun  (*  #  elements  in  a  subarray  *), 

xspc»yspc  (*  spacing  between  elenents  *), 

xphase^phase  (*  phase  taper  between  elements  *), 

kdx  .kdy.ganmax  .garnnay  (*  array  factor  components  *), 

rum ,derw,r»my .deny  (*  Interim  results  *), 

AFx My  (*  interim  array  factors  *):array[L*10]  of  real; 
title  .subtitle  (*  polar  plot  labels  *):str1ng[20]; 
flag  (*  plot  first  point  flag  *):boolean; 


myflle  (*  dunny  file  *):file; 

(********************) 

§$I  Axustcm.past 
§$I  Ainolse.past 
§SI  A:integer.past 
§$I  A:real.past 

(********************) 

(*  procedure/ function  declaration  area  *) 

(********************} 

(*  define  direction  cosines  *) 
function  kx (theta, phi  .-real) :rea1 ; 
begin 

kx  :=sin(theta)*cos(phi ) ; 

end; 

(**) 

function  ky( theta, phi  :real):real ; 
begin 

Ky:=sin(theta)*sin(phi); 

end; 

(**) 

function  kz( theta, phi  :real ) :real ; 
begin 

kz  :=cos(theta) ; 

end; 

(*  define  both  components  of  elenent  factor  *) 
function  ZEFtheta( theta, phi  :real):neal; 
begin 

ZEFtheta  :=-sin(theta) ; 

end; 

(**) 

function  ZEFph i( theta, phi  :real):real ; 
begin 

ZEFph i  :*0.0; 

end; 

(**) 

function  YEFtheta( theta  .phi  :real )  :real ; 
begin 

YEFtheta  :<os(theta)*sin(ph1 ) ; 

end; 

(**) 

function  YEFphi  (theta  ,phi  :real )  :real ; 
begin 

YEFph1:*»s(phi); 

end; 

(**) 

function  XEFtheta( theta, phi  :real):real ; 
begin 

XEFtheta  :»cos(theta)*cos(phi ) ; 

end; 

(**) 

function  XEFph  1( theta  ,phi  :rea1):real ; 
begin 

XEFphi  :—sin(phi) ; 
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end; 

(*  define  pattern  factors  *) 
function  pa£teml(1ength,k:real):real; 
begin 

pattern  1: =1.0; 

end; 

(**) 

function  pattem2(1eng  ,k  :real ) :real ; 
begin 

if  k=0  then  pattem2:=1.0  else 
pattem2:=leng*sin(pi*leng*k)/(pi*leng*k) ; 

end; 

(**) 

function  pattem3(leng,k:real )  :real ; 
begin 

if  k=0  then  pattem3:=1.0  else 
pattem3:=$qr(sin(pi*1eng/2*k)/(pi*leng/2*k))*leng/2; 

end; 

(**) 

function  pattem4(leng  ,k  :real )  :real ; 
begin 

pattem4:=cos(pi*leng*k)*(2*leng/pi  )/(l-sqr(2*leng*k) ) ; 

end; 

(**) 

function  pattem5(leng,k:real):real ; 
begin 

pattem5:=(sin(pi*leng*k)/(l-$qr(leng*k)))*leng/pi ; 

end; 

(**) 

function  pattem6(leng,k:real):real ; 
begin 

if  k=0  then  pattem6:«1.0  else 

pattem6:=(leng/2)*(sin(pi*leng*k)/(pi*leng*k))/(l-sqr(leng*k)); 

end; 

(**) 

function  pattem7(leng,k:real):real ; 
begin 

pattem7  :=sqr(cos(pi*leng/2*k)  )/(l-sqr(  leng*k)  )*2*leng/pi ; 

end; 

(*  get  nunerical  values  for  direction  cosines  and  element  factors  *) 

procedure  casey; 

begin 

case  cd  of 
l:begin 

k:=kx(  theta, phi); 

EFtheta  :=XEFtheta( theta ,phi ) ; 

EFphi:=X£Fphi  (theta  .phi); 

end; 

2:beg1n 

k  :*Ky(  theta,  phi); 

EFtheta  :*YEFtheta(  theta, phi ) ; 

EFjphi  :*YEFphi  (theta  .phi) ; 

end; 
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3  :begin 


k:=kz(theta,phi); 

EFtheta  :=ZEFtheta( theta, phi ) ; 
EFphi  :=ZEFphi(theta,phi); 


- - » 

case  dist  of  (*  get  nunerical  values  for  pattern  factors  *) 
l:PF:=patteml(leng,k); 

2:PF:=pattem2(1eng,k); 

3:PF:=pattem3(leng,k); 

4:PF:=pattem4(leng,k); 

5:PF:=pattem5(1eng,k); 

6:PF:=pattem6(1eng  ,k) ; 

7:PF:=pattem7(leng,k); 


(*  pattern  computation  procedure  *) 

procedure  compute; 

begin 

AF1:»1.0; 

(*  compute  x  array  factor  *) 

for  i :=1  to  x  do 

begin 

ganmax[i]:=kdx[i]*kx(theta,ph1)4xphase[i]; 
nuw[i]:=sin(xnun[i]#gamnax[i]/2) ; 
denx[  i  ] :  =xnun[  i  ]*si  n  { gamnax[  i  ]/2 ) ; 
if  denx[i]O0  then 
begin 

AFx[  i  ] :  =abs  (  iiihk[  i  ]/denx[  i  ] ) ; 

AFl:=AFl*AFx[i]; 

end; 

end; 

AF2:=1.0; 

(*  compute  y  array  factor  *) 

for  j:=l  to  y  do 

begin 

ganmayij] :  =kdy[  j]*ky  ( theta  ,ph  i )  -typhase[  j] ; 
mmy[  j]:=sin(ynun(  j]*ganmay[  j]/2) ; 
deny[  j]  :=>niin[  j]*sin(gamnayC  j]/2) ; 
if  deny[j]O0  then 
begin 

AFy[j]:=abs(nuny{  j]/deny[j]) ; 

AF2:*AF2*AFy[j]; 

end; 

end; 

AF:=AF1*AF2;  (*  total  array  factor  =  product  of  individual  array  factors  *) 
(*  use  pattern  multiplication  to  get  gain  pattern  *) 

Gtheta  :=AF*PF*abs{  EFtheta) ; 

(*>hi:*AF*PPabs(EFph1); 

end; 

(*  routine  to  draw  6  polar  plots  *) 
procedure  plot  it; 
var  loop,arg  :real ; 
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begin 
m:=1.0; 
while  mOO  do 
begin 

hires; 

hirescolor(lightblue) ; 
custcmdraw(0, 100,639, 100,1) ; 
custandraw( 320, 0,320, 199,1) ; 

(*  choose  appropriate  plane  .component, and  titles  *) 
case  n  of 

l,4:begin 

theta  :=pi/2y>hi  :=0; 

title:='1HETA  =  pi/2  PL AHE'; 

if  n=l  then  subtitles' theta  component'  else 

subtitle:='phi  component1; 

end; 

2,5:begin 

theta  :=0;phi  :=0; 

title:=' PHI  =  0  PLANE' ; 

if  n=2  then  subtitle :=' theta  component'  else 

subtitles' phi  component'; 

end; 

3,6:begin 

theta  :=0;phi  :=pi/2; 

title;*' PHI  =  pi/2  PUVC'; 

if  n=3  then  subtitle:*' theta  component'  else 

subtitle:*' phi  component'; 

end; 

end; 

writeln(title); 

writeln(subtitle); 

inc:=0.C6; 

loop:=0; 

flag:=false; 

(*  sweep  angle  and  perform  plotting  *) 

while  loop<=2*pi  do 

begin 

if  ( (n=l)  or  (n=4))  then  phi:=phi+inc  else  theta  :=theta+inc ; 

casey; 

compute; 

if  n<=3  then  arg:=Gtheta  else  arg:=Q)hi; 
xpl ot :  =320«Tound  ( arg*cos  ( 1  oop-pi  /2  )*fx*m ) ; 
yplot  :=100-round(arg*sin(loop-pi/2)*fy*nv,2) ; 

(*  plot  first  point,  then  draw  lines  between  successive  points  *) 

if  flag* false  then 

begin 

flag:=tnje; 

custcmplot(xplot  .yplot ,  1) ; 
oldb<  :*xp!ot  ;oldy  :=yplot ; 
end  else 

begin 

custandrav^oldx  ,oldy,xplot  .yplot  ,1) ; 
ol  cU :  *xpl  ot  »1  dy :  *yp1  ot ; 
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end; 

loop:=loop+inc; 

end; 

repeat 

*rite{'  ENTER  MAGNIFICATION  FACTOR  :'); 

Get  Real; 
mag: -real  data; 

if  not  (mag>=0)  then  InputError:=true; 
if  Input  Error=true  then  noise; 
until  InputError=false; 
m:=nrtnag; 
end; 
end; 

(**♦*♦******★) 

(*  main  array  computation  section  *) 

(************} 

begin 

legal  integer  :=£'0'..' 9']; 

fy:=200;fx:=440; 

clrscr; 

assign(itiyfile,' a  master.com' ) ; 
repeat 

write('  ENTER  LENGTH  OF  ELEMENTS  (units  of  wavelength)  :'); 

Get  Real ; 

1  eng  :=real  data; 

if  not  ((leng>0)  and  (leng<4)}  then  InputError:=true; 
if  InputError^true  then  noise; 
until  InputError=false; 
repeat 

write('  ENTER  VECTOR  CURRENT  DIRECTION  (l=x,2=y,3=z)  :') 
Get  Integer; 
cd  :=integerdata ; 

if  not  (cd  in  [1..3])  then  InputError:=true; 
if  Input Error=t rue  then  noise; 
until  InputError=false; 
repeat 

writeln  ;writeln ; 
writeln('  1.  UNIFORM'); 
writeln('  2.  RECTANGULAR  PULSE'); 
writeln('  3.  TRIANGULAR  PULSE'); 
writeln('  4.  COSM'); 
writeTn('  5.  Silt'); 
writelnC  6.  COSIfC-SqUARED'); 
writeln('  7.  ABS(SIf€)'); 

writeln ;writeln('  SPECIFY  CURRENT  DISTRIBUTION  :  '); 

Get  Integer; 
dist:=integerdata; 

if  not(dist  in  [1..7])  then  InputError  :=true; 

If  Input Error=true  then  noise; 
until  InputError=false; 
writeln; 
repeat 

wr1te('  ENTER  HMJER  OF  X-OIRECTED  SUBARRAYS  :'); 


Get  Integer; 
x  :=integendata ; 
if  x>10  then  InputError  :=true; 
if  Input  Error=true  then  noise; 
until  InputError=false; 
repeat 

writeC  ENTER  MJM3ER  OF  Y-DIRECTED  SLBARRAYS  :'); 

Get  Integer; 
y:=integerdata; 
ify>10then  InputError  :=true ; 
if  InputError=true  then  noise; 
until  InputError=false; 
for  i :=1  to  x  do 
begin 

repeat 

writeC  ENTER  #  ELEMENTS  IN  X* ,i ,*  :'); 

Get Real ; 

xnun[i]:=nealdata; 

if  (xnun[i]<=0)  then  InputError :=true ; 
if  Input Error=true  then  noise; 
until  InputError=false; 
repeat 

writeC  ENTER  ELEMENT  SPACING  FOR  X\i,'  (units  of  wavelength)  :'); 
Get  Real ; 

xspc[i]:=nealdata; 

if  not  ((xspc[i]>0)  and  (xg>c[i]<l))  then  InputError :=true ; 
if  Input Error=true  then  noise; 
until  InputError=false; 
repeat 

writeC  ENTER  PHASE  TAPER  FOR  X',i,'  (units  of  pi)  :’); 

Get  Real ; 

xphase[i]:=realdata; 

if  not  ((abs(xphase[i])>=0)  and  (abs(xphase[i])<=2))  then 
Input  Error  :=true; 
if  Input Error=true  then  noise; 
until  InputError=false; 
kdx[i]:=2*pi*xspc[i]; 
xphase[i]:=xphase[i]*pi ; 

end; 

writeln; 

for  j :=1  to  y  do 

begin 

repeat 

write('  ENTER  #  ELEMENTS  IN  Y'j,'  :'); 

Get Real ; 

ynun[j]:=realdata; 

if  (ynun[j]<=0)  then  InputError :=true; 
if  Input Error=true  then  noise; 
until  InputError=false; 
repeat 

writeC  ENTER  ELEMENT  SPACING  FOR  Y',j,'  (units  of  wavelength)  :'); 
Get  Real; 

yspc[j]:=realdata; 


if  not  ((yspc[j>0)  and  (yspc[j]<l))  then  InputError:=true; 
if  InputEiror=tn»e  then  noise; 
until  InputError=false; 
repeat 

write('  ENTER  PHASE  TAPER  FOR  Y'J,*  (units  of  pi) 

Get  Real; 

yphase[ j ] : =rea 1 data ; 

if  not  ((abs(>phase[j])>=0)  and  (abs(yphase[j])<=2))  then 
Input  Error  :=true; 
if  InputError=true  then  noise; 
until  Input  Errors  false; 
kdy[j]:=2*pi*yspc[j]; 
yphase[j]:=yphase[j]*pi ; 

end; 

for  n:=l  to  6  do  plotit; 
execute(myfile) ; 
end. 

(**) 

(*  compute  atmospheric  attenuation  *) 
function  air(fneq,elev:neal):real; 
var  x  (*  normalized  frequency  *), 
y  (*  log  of  attenuation  *), 
f  (*  log  of  frequency  *), 

atten  (*  frequency-dependent  part  of  atmospheric  attenuation  *):real 
begin 

f:=loglO(freq); 

if  ((f>=0)  and  (f<1.34))  then 

begin 

x:-f/1.4; 

if  ( ( f>=0)  and  (f<0.48))  then  y:=7.823e-2*f-1.504; 
if  ((f>=0.48)  and  (f<1.0))  then  y:=3.223e-l*f-1.614; 
if  f>=1.0  then  y:=52.3n>wr(x,3)-1.119e2*sqr(x)+80.67*x-20.86; 

end; 

if  ((f>=1.34)  and  (f<1.6))  then 
begin 

x:=(f-1.3)/0.5; 
y : =3. 323*sqr  ( x )  -2. 42*x- 0- 175; 

end; 

if  ((f>=1.6)  and  (f<=1.78))  then 
begin 

x:=(f-1.6)/0.2; 

y:=3.25*sqr(x)-0.1298*x-0.425; 

end; 

if  ((f>=1.78)  and  (f<=2.0))  then 
begin 

x:=(f-I.78)/0.22; 

if  ((f>=1.78)  and  (f<1.81))  then  y:=-3.331*f+7.979; 
if  ((f>=1.81)  and  (f<1.82))  then  y:=-1.031e2*sqr(x)+10.15*x+2.483; 
if  f>=1.82  then  y:=5.114*pwr(x,4)-16.58*pwr(x,3)+20.43*sqr(x)- 
11.19*x+2.373; 

end; 

atten  :=pwr(10»y); 
air:=-atten/sin(e1ev*dr) ; 


(*  compute  azimuth  and  elevation  angles  *) 
procedure  azel ; 

var  alt  (*  satellite  altitude  *), 

latl  (*  latitude  of  first  ground  site  *), 
longl  (*  longitude  of  first  ground  site  *), 
lat2  (*  latitude  of  second  ground  site  *), 
long2  (*  longitude  of  second  ground  site  *) , 
lats  (*  latitude  of  subsatellite  point  *), 
longs  (*  longitude  of  subsatellite  point  *), 
rl,r2  (*  ground  site-subsatell ite  point  great  circle  distances  *), 
tiltl,tilt2  (*  tilt  angles  *), 
p  (*  altitude  normalized  w/resp  to  Earth  radius  *), 
azl,az2  (*  azimuth  angles  *), 
elevl,elev2  (*  elevation  angles  *):real; 
begin 

textcolor(viiite) ; 

textbackground(blue) ; 

clrscr; 

getaltitude; 

alt:=realdata; 

p:=l+alt/te; 

writeln;writeln('  ENTER  C00RDIN(\TES  OF  FIRST  GROUfC  SITE  (deg)1); 
repeat 

writeC  LATITUDE  ; 

Get  Real ; 
latl:=nealdata; 

if  (abs(latl)>90)  then  InputError  :=true; 
if  InputError=true  then  noise; 
until  InputError=false; 
latl:=latl*dr; 
repeat 

writeC  LONGITUDE  : 

Get Real ; 

longl  :=real  data; 

if  (abs(longl)>90)  then  InputError :=true; 
if  Input Error= true  then  noise; 
until  InputError^false; 
longl  :=longl*dr; 

writeln;writeln('  ENTER  COORDIiWTES  OF  SECOND  GROUND  SITE  (deg)  :'); 
repeat 

writeC  LATITUDE  :  '); 

Get  Real ; 
lat2:=nealdata; 

if  (abs(lat2)>90)  then  InputError :=true; 
if  Input Error=  true  then  noise; 
until  InputError=false; 

Iat2:=lat2*dr; 

repeat 

writeC  LONGITUDE  :  '); 

Get  Real ; 

long2:=realdata; 


if  (abs(1cng2)>90)  then  InputError:=true; 
if  Input Error=t rue  then  noise; 
until  Input Error=false; 
long2:=long2*dr; 

writeln;writeln('  ENTER  COORDIMTES  OF  SUB-SATELLITE  POINT  (deg)'); 
repeat 

write('  LATITUDE  :  '); 

Get  Real; 
lats:*realdata; 

if  (abs(lats)>90)  then  IhputError:=true; 
if  InputError=true  then  noise; 
until  InputErnor=false; 
lats:=lats*dr; 
repeat 

write('  LONGITUDE  :  '); 

Get  Real; 

longs  :=real  data; 

if  (abs(longs)>90)  then  InputError  :=true; 
if  Input Error=true  then  noise; 
until  InputError=false; 
longs  :=longs*<fr; 

rl:=arocos(cos(latl-lats)*cos(longl-longs) ) ; 
r2:=arccos(cos(lat2-lats)*cos(long2-longs)); 

(*  take  all  cases  into  account  and  adjust  azimuth  accordingly  *) 
if  latl=lats  then  if  longl-longs<0  then  azl:=9G.O  else  azl:=270.0 
else  azl:=rt*arctan(tan(longl-longs)/sin(latl-lats)); 
if  lat2=lats  then  if  long2-longs<0  then  az2:=90.0  else  az2:=270.0 
else  az2:=rd*arctan(tan(long2-longs)/sin(lat2-lats)); 
if  latl-lats>0  then  azl:=azl+180i0; 
if  lat2-lats>0  then  az2:=az2+180.0; 
if  azKO  then  azl:=azl+36Q.O; 
if  az2<0  then  az2:=az2+360.0; 
tiltl:=arctan(sin(rl)/(p-cos(rl))); 
ti 1 t2:=arctan(sin( r2 )/ (p-oos( r2 ) ) ) ; 
elevl:=arccos(p*sin(tiltl)) ; 
elev2:=arccos(p*sin(tilt2)) ; 

(*  store  elevation  angles  and  ground  site  latitudes  *) 

data.alpha.upl  :=rd*elevl; 

data.alpha.dnl  :=rd*relev2; 

data.1at.up1  :*latl; 

data.lat.cM  :=lat2; 

havealpha  :=true; 

textcol  or(l  ightmagenta ) ; 

wrlteln; 

wr1teln('  FIRST  GROUND  STATION' ); 
wrlteln ;writeln('  azimuth*  ',azl:0:2,'  deg'); 

wr1teln(’  elevation  *  ' ,elevl*rd:0:2,'  deg'); 

wrlteln  ;writeln;writeln('  SECCH)  GROUN)  STATION'); 
wrlteln  ;wr1te1n('  azlnuth  *  ',az2:0:2,'  deg'); 

wr1te1n('  elevation  *  ' ,elev2*nd:0:2,'  deg'); 

end; 

(**) 

(*  compute  steady-state  temperature  *) 


procedure  blackbody; 

const  PHIsun=3.91e26;  (*  average  flux  anitted  from  sun  *) 
var  z  (*  satellite  altitude  *), 

TotArea  (*  total  satellite  surface  area  *), 

Exp  Area  (*  percent  of  total  area  exposed  to  sun  *) , 

Eorbit  (*  irradiance  at  the  satellite's  orbit  *), 

Tsat  (*  steady-state  satellite  temperature  *):real; 
begin 

textbackground(blue) ; 

clrscr; 

getaltitude; 

z:=realdata; 

repeat 

writeln ; 

writeC  ENTER  TOTAL  SATELLITE  SURFACE  AREA  (itP2) 

Get  Real; 

TotArea:=realdata; 

if  TotArea<=0  then  InputError  :*true ; 
if  Input Error=true  then  noise; 
until  InputError=false; 
repeat 

^iteln  * 

writeC  ’ENTER  PERCENT  OF  TOTAL  SURFACE  AREA  EXPOSED  TO  SUN 
Get  Real; 

ExpArea  :=realdata; 

if  not  ((ExpArea>0)  and  (ExpArea<*100))  then  InputError  :=true; 
if  Input  Error-true  then  noise; 
until  InputError=false; 

ExpArea  : -ExpArea/ lOOMbtArea ; 

Eorbit  :-PHIsun/(4*p1*sqr(1.5ell-z)) ; 

Tsat : -sqrt  ( sqrt  ( Elorb  i  t*ExpArea/  ( Tot  Area*5. 67e-8 ) ) ) ; 
textcolor(l  ightmagenta) ; 

writeln ;writelnjwriteln('Steady-State  Temperature  =  '  ,Tsat:0:2,'  deg  K') 
end; 

(**) 

(*  compute  Shannon  rate  and  fyquist  rate  *) 
procedure  capacity(bandwidth,l¥,Teq,rate^:real); 
var  shannon  (*  Shannon  rate  *) , 
nyquist  (*  Hyquist  rate  *)  :real ; 
begin 

sharwonflag:=false; 
nyquistflag: -false; 

shannon  :=bandwidth*le6*log2(l-wlog(Pr)/(bandwidth*Ie6*lc*Teq)  )/le6; 
nyquist  :-2*bandw1dth*le6*log2(M)/  le6; 

(*  use  Shannon  rate  if  It  is  limiting  factor  *) 

if  shamorKnyquist  then 

begin 

limit : -shannon; 
shannonflag:*true; 
end  else 

(*  use  f^qulst  rate  If  it  is  limiting  factor  *) 
begin 

limit  :*«\yqu1st; 


nyqu1stflag:*true; 

end; 

(*  otherwise  use  previously  computed  data  rate  *) 

if  1  tmit>rate  then 

begin 

limit  :*rate; 
s  ha  nncn  fl  ag:=*  false; 
nyquistflag  :*false; 

end; 

end; 

(**) 

(*  compute  circular  orbit  parameters  *) 

procedure  circular; 

var  h  (*  satellite  altitude  *), 

alpha  (*  mlninun  elevation  angle  *), 
theta  (*  subtended  angle  *), 
velocity  (*  satellite  velocity  *), 
period  (*  orbital  period  *) , 
coverage  (*  area  of  coverage  *) , 
le  (*  Earth's  rotation  period  *), 
visibility  (*  duration  of  satellite  visibility  *);real; 
begin 

textbackground(blue) ; 
clrscr; 
getaltitude; 
h:=realdata; 

velocity  :=sqrt(mj/ (Re+h) ) ; 

period  :=2*pi/sqrt(mu  )*pwr(Re4h,  1. 5)/3600; 

repeat 

writeln; 

writeC  ENTER  HINIHLH  TEFMINAL  ELEVATION  ANGLE  (deg)  :'); 

Get  Real; 

alpha  :*realdata; 

if  ((alpha<=0)  or  (alpha> 360.0))  then  LnputError:=true; 
if  Input Error=true  then  noise; 
until  InputError=false; 
repeat 

writeln; 

vrite('  [R]etrograde  or  [Pjrograde  orbit  ?  '); 
readln( letter) ; 

if  not  (letter  In  [' p'  ,'P'  ,'r'  ,'R'])  then  noise; 
until  letter  in  ['p'  ,'P'  ,'r'  ,'R'j; 

If  letter  in  ['p','P']  then  le:=-24.0  else  Te:«24.0; 
theta  :^rccos(R»*cos(alpha*dr)/(Re+h)  )*rd-alpha ; 
coverage  :*2*p1*sqr(Re)*(l-cos(theta*dr)  )/le6; 
v1sib111ty:«2*theta*per1od/(360*(l+period/Te)); 
textcolor(  1  Ightmagenta ) ; 
writeln  ;*r1teln; 

wr1teln( 'Velocity  «  ' ,veloc1ty:0:2,'  ny'sec'); 
wr1teln( 'Orbital  Rjriod  «  ' ,per1od;0:2,'  hr'); 
writeln  ('Coverage  Area  «  '  .coverage  :8,'  sq  ton'); 
write(' Duration  of  Visibility  ■  '); 

If  iw3.5784e7  then  »r1teln('CGNTIkJ0US')  else  writeln(vis1b111ty:0:2,’  hr'); 
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(**) 

(*  compute  elliptical  orbit  parameters  *) 
procedure  elliptical; 
var  apogee  (*  satellite  apogee  *), 
perigee  (*  satellite  perigee  *), 
aj>  (*  major  and  minor  axes  of  orbit  *), 
r,phi  {*  polar  coordinates  of  orbiting  satellite  *), 
alpha  (*  elevation  angle  *), 
theta  (*  subtended  angle  *) , 

VA,VP  (*  velocities  at  apogee  and  perigee  *), 

Vavg  (*  average  velocity  *), 

velocity  (*  velocity  at  a  specified  reference  angle  phi  *), 
eccentricity  (*  orbit  eccentricity  *), 
period  (*  orbital  period  *), 
coverage  (*  coverage  area  *)  :real ; 
begin 

textbackground(blue) ; 

clrscr; 

repeat 

^iteln  * 

writeC  ’ENTER  ALTITUDE  @  APOGEE  (m) 

Get  Real ; 

apogee  :-rea1data ; 

if  not  (apogee>0)  then  InputError  :=true ; 
if  Input Error=true  then  noise; 
until  InputError=false; 
apogee :=apogee+fte; 
repeat 

write ( '  ENTER  ALTITUDE  @  PERIGEE  (m)  :'); 

Get  Real; 

perigee  :«real  data; 

if  not  (perigee>0)  then  InputError  :=true ; 
if  InputError=true  then  noise; 
until  InputError=false; 
perigee :*perigee+R»; 
a  :*(apogee4perigee)/2; 
b:=sqrt(apogee*perigee) ; 
eccentricity  :=sqrt( l-(sqr(b)/sqr(a) ) ) ; 
va  :*sqrt(mu/a*{  ( 1-eccentric  ity)/  ( 1+eccentricity) ) ) ; 
vp  :*sqrt(nu/a*(  ( l-*eccentricity)/  { 1-eccentri  city) ) ) ; 
repeat 

writeln; 

writeC  ENTER  REFERENCE  ANGLE  (deg)  :'); 

Get  Real; 
phi  :=realctrta; 

if  not  (phi>«0)  and  (phK-360)  then  InputError  :»true; 
if  In  put  Error- true  then  noise; 
until  IhputError*false; 
vavg:«sqrt(vp^va); 

r  :«a*(  l-sqr(eccentr1c  ity) )/  ( l-*eccentr1c1ty*cos(ph1*dr) ) ; 

velocity  :«sqrt((2/r-l/a)*mu) ; 

period  :«2*pi*pwr(a ,  1. 5)/sqrt(n»j)/3600; 


m 


repeat 

wrlteln  * 

write('  ’ENTER  TEW  INAL  ELEVATION  ANGLE  (deg)  :*); 

Get  Real; 

alpha:*realdata; 

if  not  ((a1pha>-0)  and  (alpha<=360))  then  iiputError:=true; 
if  Input  Error* true  then  noise; 
until  InputError=false; 

theta  :=arocos(Re*cos(alpha*dr)/(Re+r) )*rd-alpha ; 
coverage:=2*pi*sqr((te)*  ( l-cos(  theta*dr ) )/  le6; 
textcolor(lightmagenta) ; 
writeln  ;wr1teln ; 

writelnC Velocity  @  Apogee  =  ',VA:0:2,'  nVsec'); 

writeln(' Velocity  @  Rerigee  =  ',VP:0:2,'  nVsec'); 

writeln( 'Average  Velocity*  ',Vavg:0:2,'  nVsec'j; 

if  phi OO  then  writeln( 'Velocity  @  phi  =  ',phi:0:0,'  deg  =' ,velocity:0:2, 

'  m/sec'); 

wri  tel  n  ( '  Eiccentric  ity  =  ’  eccentricity  :0:2); 
vriteln( 'Orbital  Period*  '  .period :0:2,'  hr'); 

*riteln( 'Coverage  Area  =  '.coverage#,'  sq  Ion'); 
end; 

(**) 

(*  confute  escape  velocity  at  altitude  *) 

procedure  EscapeVelocity; 

var  2  (*  launch  site  altitude  *), 

Vesc  (*  escape  velocity  *):real; 
begin 

textbackground(blue) ; 

clrscr; 

repeat 

an  tel  n  * 

write(' ’ENTER  LAUNCH  SITE  ALTITUDE  (m)  :'); 

Get  Real; 
z:*realdata; 

If  z<0  then  InputError  :*true ; 
if  Input  Err  or=true  then  noise; 
until  InputError=false; 
vesc  :=sqrt(2*tru/  ( Re+z ) ) ; 
textcolor(l  ightmagenta) ; 

writelnpvriteln^ritelnC Escape  Velocity*  \Vesc:0:2,'  nVsec'); 
end; 

(**) 

(*  confute  gain  for  horn  *) 
function  ga1nl(d,f:rea1):real; 
begin 

gainl  :«4*p1*sqr(d*f*Le9/3e8) ; 

end; 

(**) 

(*  confute  beamwldth  for  horn  *) 
function  bwl(d,f:rea1):real; 
begin 

bv»l:*0.8B*3e8/(f,le9*d)*rd; 

end; 
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(*  compute  gain  for  helix  *) 
function  ga1n2(d,l  ,f  :rea1 )  :real ; 
begin 

gain2:*15.0*pi*d*l*sqr(f*le9/3e8); 

end; 

(**) 

(*  compute  beamwidth  for  helix  *) 
function  bw2(d,l  ,f:  real):  real; 
begin 

bw2:=52.0»3e8/(f*le9*p1*d*sqrt{1 )) ; 

end; 

(**) 

(*  compute  gain  for  parabolic  dish  *) 
function  gain3(d,f  :real):real; 
begin 

gai  n3  :=sqr(p1*d*f*le9/3e8) ; 

end; 

(**) 

(*  compute  beamwidth  for  parabolic  dish  *) 

function  bw3(d,f:real):real; 

begin 

bw3:=l.Q2*3e8/(f*le9*d)*rd; 

end; 

(**) 

(*  compute  gain  for  phased  dipole  array  *) 
function  gai n4(n,d,f: real ): rea  1 ; 
begin 

gain4:=pi*n*d*f*le9/(1.4*3e8) ; 

end; 

(**) 

(*  compute  beamwidth  for  phased  dipole  array  *) 

function  bw4(n,d,f:real):real; 

begin 

bw4:=50.0*3e8/(n*d*f*le9); 

end; 

(**) 

(*  compute  required  beamwidth  for  global  coverage  *) 

function  bw5(alt  ,elev  :real )  :real ; 

begin 

bw5:*(2*arcsin(Re#cos(elev*dr)/(Re+alt) ) )*rd; 

end; 

(**) 

(*  compute  required  gain  for  global  coverage  *) 
function  gain5(bw:real):real; 
var  c:real; 
begin 

c:«1.27;  (*  typical  value  *) 
ga1n5:*sqr(p1*c/(b*dr)); 

end; 

(**) 

(*  compute  required  beamwidth  for  spot  coverage  *) 
function  bw6(alt  .nadir  :real):real; 


begin 

bw6:»(2*arctan(nadir/(2*alt)))*rd; 

end; 

(**) 

(*  compute  required  gain  for  spot  coverage  *) 
function  gain6(alt,nadir:real):real; 
var  c:nea1; 
begin 

c:*1.27;  (*  topical  value  *) 
gain6:=sqr(pi*c*alt/ nadir) ; 

end; 

(**) 

(*  compute  EIRP  *) 

function  isotropic(Rxit  Jx&1n,TxCircuitLoss :rea1 ) :real ; 
begin 

isotropic  :=10*loglO(Rout)+TxG&in+TxCircuitloss ; 

end; 

(**) 

(*  confute  equivalent  temperature  *) 

function  eqtap(l&MTemp,NoiseFigure,ffcCincuitLoss  :real )  ‘.real ; 
begin 

eqtenp  :=ftxAitTanp+(MDiseFigure/alog(RxCi  rcuit  Loss)-l.  0)*29G. 0; 

end; 

(**) 

(*  compute  (VT  figure  of  merit  *) 
function  fignerit(lfc&rin,Teq:real):rea1; 
begin 

fignerit  :=RxG&in-10*logl0(Teq) ; 

end; 

(**) 

(*  compute  received  useful  power  *) 

function  power(Rout  JxC1rcuitLoss,Tx&1n,TxEffic1ency,Rath|jDSS,AirAtten, 
RainAtten , Rx Ef f ici ency ,  fa  Gain ,  faCi  rcu  it  Loss  :real )  :real ; 

begin 

power  :=10*logl0(Pout)+TxCircuitloss+TxGbin+10*logl0(TxEfficiency)+ 
PathLoss+AirAtten+fiainAtten+10Mogl0(faEfficiency)+fa&in+ 
faCircuitloss; 

end; 

(**) 

(*  compute  CM  *) 

function  cnratio(Pr,Teq:real):real; 
begin 

cnratio  :*Pr-10*logl0(k*Teq) ; 

end; 

(**) 

(*  compute  available  Eb/No  *) 
function  EbftAvail(Pr,Teq,datarate:rea1):real; 
var  to  (*  noise  power  density  *):real; 
begin 

fb:«k*Teq; 

Ebtofcall  :«Pr-10*logl0(N)*datarate*le6); 

end; 


(*  compute  bit  error  rate  *) 

function  bi  terror  (modtype  integer  jn,Arai  1  EbH)  :real )  :real ; 
begin 

/Wall  Ebfb  :*alog(/Va1  1  Ebfto) ; 

If  (((m>l)  and  (/VailEMb*log2(m)<=86.0))  or  (/VailEbftK=86.0))  then 
case  modtype  of 

1 ,2,3,4:biterror  :*q(sqrt(2*Avai  1  Ebfto) ) ; 

5:biterror:=q(sqrt{2*  Avail  £bNo*1og2(m))*sin(pi/m))/log2(m) ; 
6,9:biterror  :=q(sqrt(Avai  1  Ebto) ) ; 

7:biterror  :=Q.5*exp(-0.5*Avail  EbNo) ; 

8:biterror  :=«V4.0*exp(-Avai  1  Ebfb*log2(m) ) ; 

10:biterror  :=((m-l)/nf  q(sqrt(6Mog2(m)*  Avail  EbNo/(sqr(m)-1.0) ) ))/ 
log2(m); 

11  :biterror  :=0. 5*exp(-Avai  1  Eblto) ; 
end  else  biterror :=0; 
nosound; 
end; 

(**) 

{*  confute  required  Eb/No  *) 

function  EbN3Req(modtype:integer;m,BER:real):neal; 

var  Ebhb  :real ; 

begin 

case  modtype  of 

l,2,3,4:Eblt>:=sqr(q1nverse(BER))/2.0; 
5:Ebfb:=sqr(qinverse(log2(m)*BER)/sin(pi/!i)))/(2.0*log2(m)); 
6,9:EbNa  :=sqr(qinverse{BER)) ; 

7  :Ebfb  :=-2.0*ln(2.0*BER) ; 

8:EbNo  :*-ln(4.(Vhi*BER)/log2(m) ; 

10:EbNo:=(sqr(m)-1.0)/(6.0*log2(m))*sqr(qinverse(m7(m-l)* 

log2(m)*BER)) ; 

ll:EbNo:=-ln(2*BER); 

6nci  * 

EbftoReq:=10*logl0(Ebto); 

nosound; 

end; 

(**) 

(*  compute  max  datarate  *) 
function  maxrate(Pr,Teq,ReqEbft:real):real; 
var  It)  (*  noise  power  density  *)  :real ; 
begin 

fto:*k*Teq; 

maxrate  :=(alog(Pr)/  {No*alog(ReqEblt)))  )/le6; 

end; 

(**) 

(*  compute  rainfall  attenuation  *) 
function  ra1n(rainrate,f,elev,alt,lat :real )  :real ; 
var  x  (*  normalized  frequency  *), 
y  {*  dum\y  variable  *), 
ajb  (*  coefficients  a(f)  and  b(f)  *), 
gamna  (*  empirical  constant  *), 
fteltlayer  (*  melting  l^yer  height  *), 

StormHelght  (*  effective  storm  height  *) , 
tothlength  (*  path  length  through  storm  *) , 


EffDist  (*  total  effective  distance  *):real; 
begin 

If  ra1nrate>0  then 
begin 

(**************) 

(*  calculate  a(f)  and  b(f)  *) 

(**************) 
if  rainrate<=25.0  then 
begin 

if  f<=10.0  then 
begin 

x:=(f-1.0)/9. 0; 

a  :=3. 475e-3*pwr(x ,6)-l. 802e-2*pwr(x,5)+2. 76Be-2*pwr(x ,4)- 
1. 036e-2*pwr(x  ,3) +7.878e-3*sqr(x)+9. 815e-4*x+€.  50Be-5 ; 
end  else 
begin 

x:=f/100.0; 

a  :=~5.231*pwr(x  ,6)+20. 12*pwr(x  ,5)-28. 12*pwr(x  ,4)+15. 28*pwr(x  ,3)- 
1. 067*sqr(x)  +0i  1094*x-l.  707e-3 ; 
end; 

if  f<8.0  then 
begin 

x:=(f-1.0)/7.0; 

b:=1.238*pwr(x,5)-2.749*pwr(x,4)+1.495*pwr{x,3)+7.341e-2*sqr(x)+ 
Q.2389*x+Ol890B; 
end  else 
begin 

x:*(f-8.0)/92.0; 

b  :=3. 734*pinr(x  ,8)-5. 39*pwr(x  ,7)-l.  103*pwr(x  ,6)+0.  2915*pwr (x  ,5)+ 
7.045*pwr(x,4)-5.547*pwr(x,3)+1.116*sqr(x)-0^919*x+1.188; 

end; 
end  else 
begin 

if  f<=10.0  then 
begin 

x:*(f-1.0)/9.0; 

a  :*1.463e-2*pwr{x ,3)-5. 037e-3*sqr(x)  +1. 9e-3*x+4. 672e-5; 
end  else 
begin 

x:«f/ 100.0; 

a:=28.64*pwr(x,10)-49.0B*F*r(x,9)-24.49*pvr(x,8)+57.94*pwr(x,7)+ 
27. 15*pwr(x  ,6)  -37.87*pwr(x  ,5)  -27. 01*pwr(x  ,4)  +29. 23*pwr(x  ,3)- 
3. 764*sqr(x )  +0. 2253*x-2. 531e-3 ; 
end  * 

if  (*(f>«1.0)  and  (f<=6.0))  then 
begin 

x:*(f-1.0)/5.0; 

b  :«1. 719*p»r(x  ,6)-4. 176*pwr(x  ,5)+3. 0B3*pwr(x  ,4)-0. 9744*pwr(x  ,3)+ 
a  4254*sqr(x ) +a  2611*x+a  9468 ; 
end  else  if  ( ( f>6. 0)  and  (f<«30.0))  then 
begin 

x:«(f-6.0)/24.0; 

b  :«4. 418*pwr(x  ,5)-10. 51*p<r  (x  ,4)-*8. 006*pwr(x  ,3)-l.  7B3*sqr(x )- 


0.4557*x+1.288; 
end  else 
begin 

x:-(f-30.0)/7Q.0; 

b:*a5882*pwr(x,4)-1.709*pwr(x,3)+1.868*sqr(x)-0.9e7*x+0.964; 

end; 


(*****************) 

(*  calculate  effective  distance  *) 

(*****************} 

if  abs(lat)<=30.0  then  Me1tLa,yer:=4.8  else  Meltla)«r:=7.8-0,l*abs(lat) ; 
if  rainrate<=10i0  then  StormHeight:=Melt  Layer  else 
Stormlteight  :=MeltLajer+logl0(rainrate/10.0) ; 

Path  length  := ( StormHei  ght-al  t/  le3 )/ s  i  n  (el  ev*dr ) ; 
garnna:=l.(V22.0; 

if  rainrate<=lQ.O  then  EffDist:=Path length  else 
begin 

y:=garana*b*ln(rainrate/10.0)*cos(elev*dr) ; 

EffDist  :=(  l-exp(  -y*Path  Length)  )/y ; 

end; 

(*  compute  total  rainfall  attenuation  *) 
rain  :=-a*pwr(ra1nrate  ,b)*EffDist ; 
end  else  rain:=0; 
end; 

(**) 

(*  ground  station  antenna  analysis  *) 

procedure  antennal; 

var  gain  (*  antenna  gain  *), 

BW  (*  antenna  beamwidth  *), 
rho  (*  antema  efficiency  *), 
f  (*  frequency  *) , 

d,l  (*  antenna  dimensions  *)  :real ; 

Ant  Type  (*  antenna  type  *)  integer; 
begin 

textbackground(blue) ; 
clrscr; 

writeln;writeln('A  NTENNA  ANALYSIS':®); 

goto  XY(  1,5); 

writeln(title[site]:50) ; 

gotoXY(l.lO); 

vriteln('  1.  Him'); 

wr1teln('  2.  Helix'); 

writeln('  3.  fferabolic  Dish'); 

writeln('  4.  biased  Array  of  Dipoles'); 


vriteln('  1.  Him'); 

wr1teln('  2.  Helix'); 

writeln('  3.  fferabolic  Dish'); 

writeln('  4.  biased  Array  of  Dipoles'); 

repeat 

write  in  ;vriteln;vr1te(*  ENTER  ANTEJKA  TYPE  :'); 

Get  Integer; 

Ant  Type  :*1ntegerdata ; 

if  not  (Ant Type  in  [1..4])  then  InputError  :=true; 

"  InputError*true  then  noise; 
irti1  I*putError*false; 


rep(  it 


write!'  ENTER  ANTEWiA  EFFICIENCY  :'); 


IV*'. 

& 

$ 


if 


K,»V 


Get  Real; 
rho:=rea1data; 

if  not((rho>=0)  and  (rho<=1.0))  then  InputError:=true; 
if  Input  Error=t  rue  then  noise; 
until  Input&ror=  false; 

if  uplink=true  then  f:=data.freq.up1  else  f:=data.fneq.dnl ; 
case  Ait  Type  of 

lrbegin  (*  horn  *) 
repeat 

write('  ENTER  DIMENSION  OF  HORN  (m)  :'); 
Get Real ; 
d  :=realdata ; 

if  d<=0  then  InputError:=true; 
if  Input  Error=t  rue  then  noise; 
until  InputError=false; 
gain:=gainl(d,f); 

BW:=bwl(d,f); 

end; 

2:begin  (*  helix  *) 
repeat 

write('  ENTER  DIAMETER  OF  HELIX  ;’); 

Get Real ; 
d:=nealdata; 

if  d<=0  then  InputError:=true; 
if  Input Error=t rue  then  noise; 
until  InputErnor=false; 
repeat 

writeC  ENTER  LENGTH  OF  HELIX  :'); 

Get Real ; 

1  :=realdata ; 

if  1<=0  then  InputError:=true; 
if  Input Error= true  then  noise; 
until  Input Error=  false; 
gain  :=gain2(d  ,1  ,f); 

BW:=bw2(d,l  ,f); 

end; 

3:begin  (*  parabolic  dish  *) 
repeat 

writeC  ENTER  DIAMETER  OF  DISH  (m)  :'); 

Get  Real ; 
d :=realdata ; 

if  d<=0  then  InputError:=true; 
if  Input Error*t rue  then  noise; 
until  Input Error*  false; 
ga1n:=gain3(d,f); 

BW:=bw3(d,f) ; 

end; 

4:begin  (*  dipole  array  *) 
repeat 

writeC  ENTER  #  LINEAR  ELEMENTS  :'); 

Get Real ; 

n:«1nt(realdata); 

if  n<2.0  then  Input  Error  :*true; 
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if  Input Error=t rue  then  noise; 
until  Input Error=  false; 
repeat 

writeC  ENTER  ELEMENT  SPACING  (tn) 

Get  Real ; 
d:=realdata; 

if  not({d>0)  and  (d<1.0))  then  Input  Error  :=true; 
if  Input Error=t rue  then  noise; 
until  Input Error=false; 
gain:=gain4(n,d,f); 

BW:=bw4(n,d,f); 

end; 

end; 

textcolor(lightmagenta) ; 

writeln;writeln; 

gain:=10*logl0(gain); 

writeln('MAX  GAIN  =  ',gain:0:2,'  dB'); 

writeln('3-<iB  BEAMWD7H  =  ',BW:0:3,'  deg'); 

(*  store  antenna  data  *) 

if  uplink=true  then  with  data  do 

begin 

Antenna Type.upl  :=int(AntType); 
argl.upl  :=d; 
arg2.upl  :=1 ; 
arg3.upl  :=n; 

TxGain.upl  :=gain; 
beamwidth.upl  :=BW; 

TxEfficiency.upl  :=rho; 
end  else  with  data  do 
begin 

Antenna Type.dnl  :=int(AntType) ; 
argl.dnl  :=d; 
arg2.dnl  :=1 ; 
arg3.dnl  :=n; 

IfcGain.dnl  :=gain; 
ftoEfficiency.dnl  :=rho; 

end; 

textcolor(vrf)ite); 

end; 

(**) 

(*  satellite  antenna  analysis  *) 
procedure  antenna2; 
var  beam  (*  beam  type  *):integer; 
gain  (*  actual  antenna  gain  *) , 
reqgain  (*  required  antenna  gain  *), 

BW  (*  actual  antenra  beamwidth  *) , 
reqBW  (*  required  antenna  beamwidth  *), 
f  (*  frequency  *) , 
elev  {*  elevation  angle  *), 
rho  (*  antenna  efficiency  *) , 
d  (*  dish  diameter  *), 
radir  (*  desired  nadir  diameter  *)  :real ; 
begin 
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if  up1ink=true  then  f:=data.freq.upl  else  f :=data.freq.dnl ; 
if  uplink=true  then  elev:=data.a1pha.upl  else  elev:=data.alpha.dnl ; 
textbackgnound(blue) ; 
clrscr; 

writelnjwritelnCA  NTENNA  ANALYSI  S':55); 

gotoXY(l,5); 

writeln(title[site] :50) ; 

gotoXY(l.lO); 

writelnC  1.  Earth  Coverage'); 

writeln('  2.  Spot  Coverage'); 

repeat 

writelnjwritelnjwriteC  ENTER  DESIRED  BEAN  TYPE  :'); 

Get  Integer; 
beam:=integerdata; 

if  not  (beam  in  [1.  .23)  then  Input  Error  :=true; 
if  Input Error=true  then  noise; 
until  InputError=false; 
repeat 

writeC  ENTER  ANTEMiA  EFFICIENCY  :'); 

Get Real ; 
rho:=realdata; 

if  not((rho>0)  and  (rho<=1.0))  then  Input  Error  :=true; 
if  Input Error=true  then  noise; 
until  InputError=false; 
repeat 

writeC  ENTER  ANTEWA  DIAMETER  (m)  :'); 

Get Real ; 
d:=realdata; 

if  d<=0  then  Input  Error  :=true; 
if  Input Error=true  then  noise; 
until  InputError=false; 
case  beam  of 

l:begin  (*  global  coverage  *) 

reqBW  :=bw5  ( data  .al  t  i  tude  ,el  ev ) ; 
reqgain  :=gain5(reqEW) ; 

BW:=bw3(d,f) ; 
gain:=gain3(d,f) ; 

end; 

2:begin  (*  spot  coverage  *) 
repeat 

writeC  ENTER  NADIR  DIAMETER  (m)  :'); 

Get  Real ; 

nadir  ;=realdata ; 

if  nadir<=0  then  InputError:=true; 
if  Input Error=true  then  noise; 
until  InputError=  false; 
reqBW:=tw6(data.altitude, nadir) ; 
reqgain :=gain6(data. altitude .nadir) ; 

BW:=bw3(d,f); 

gain:=gain3(d,f); 

end; 

end; 

textcolor(  1  ighbnagenta ) ; 


writeln;write1n; 
gain:=10*logl0(gain) ; 
reqgain:=10*1ogl0(neqgain); 
writeln('MAX  GUN  =  '  .gain :0:2,'  dB’ ) ; 
writeln('3-dB  BE/WIDTH  =  ',BW:0:3,'  deg'); 

(*  check  if  required  parameters  are  attainable  *) 
if  (gairKneqgain)  and  (uplink=false)  then 
begin 

textcolor(yellow) ; 

writeln  ;writeln( 'NOTE:  Required  Gain  =  1  ,reqgain:0:2, 

'  dB  for  25-dB  sidelobes' ) ; 

end; 

if  (BWCreqBW)  and  (uplink=false)  then 
begin 

textcolor(yellow) ; 

writeln  ;wri  tel  n( 'NOTE:  Required  Beamwidth  =  ',reqBW:0:2, 

'  deg  for  25-dB  sidelobes'); 

end; 

(*  store  antenna  data  *) 

if  uplink=true  then  with  data  do 

begin 

beamtype.upl  :=int(beam) ; 
arg4.upl  :=radir; 
arg5.upl  :=d; 

IfcGain.upl  :=gain; 

RxEfficiency.upl  :=rho; 
end  else  with  data  do 
begin 

beamtype.dnl  :=int(beam) ; 
arg4.dnl  :=nadir; 
arg5.dnl  :=d; 

TxGain.dnl  :=gain; 

TxEfficiency.dnl  :=rho; 
beamwidth.dnl  :=BW; 

end; 

textcolor(v*iite) ; 
end; 

(**) 

(*  receiver  analysis  *) 
procedure  receiver; 
var  NF  (*  noise  figure  *), 

loss  (*  receiver  circuit  losses  *), 
temp  (*  receiving  antenna  noise  temperature  *):real; 
begin 

textbackground(blue) ; 
clrscr; 

writeln  jwriteln('R  ECEIVER  ANALYSIS'  :56) ; 
gotoXY(l,5); 

if  uplink=true  then  *riteln(title[5]:47)  else  writeln(title[6]:44); 

gotoXY(l.lO); 

repeat 

writeln; 

writeC  ENTER  RECEIVER  AMPLIFIER  NOISE  FIGURE  :’); 


Get Real ; 

NF:=rea1data; 

if  not  (NF>=1.0)  then  InputError:=true; 
if  Input  Error=t  rue  then  noise; 
until  InputError=false; 
repeat 

wpitcln * 

writeC  ’ENTER  RECEIVER  CIRCUIT  LOSSES  (-dB) 

Get  Real ; 
loss  :=realdata ; 

if  not  (loss<=0)  then  InputError:=true; 
if  Input Error=true  then  noise; 
until  Input Error=false; 
repeat 

wf*i  tel  n  * 

writeC  ’ENTER  RECEIVING  ANTOWA  TEMPERATURE  (deg  K)  :'); 

Get  Real ; 
temp  ;=real  data; 

if  not  (temp>0)  then  In  put  Error  :=t  rue; 
if  Input Error= true  then  noise; 
until  Input Error=false; 

(*  store  input  data  and  perform  computations  *) 

with  data  do 

begin 

if  uplink=true  then 
begin 

NoiseFigure.upl  :=NF; 
i&CircuitLoss.upl  :=loss; 
ffc/VitTemp.upl  :=tanp; 

Teq.upl  :*eqtemp(RxAntTemp.upl  .foiseFigune.upl  .ffcCircuitLoss.upl); 
GT.upl  :=fignerit(ffcGain.upl  .Teq.upl ) ; 
textcolor(lightmagenta) ; 

wrltelnjwritelnC  SYSTEM  EQUIVALENT  TEMPERATURE  =  \ 

Teq.upl  :0:2,'  deg  K' ) ; 

writelnC  G/T  FIGURE  OF  MERIT  =  '  .GT.upl  :0:2,'  dB'); 
end  else 
begin 

toiseFigure.dnl  :=NF; 

FfcCircultLoss.dnl  :=loss; 
ftcAntTanp.dnl  :®temp; 

Teq.dnl  :=eqtemp(ffcAntTemp.dnl  .fbiseFigune.dnl  .IfcCircuitLoss.dnl); 
GT.dnl  :=figner1t(ffcGain.dnl  .Teq.ciil ) ; 
textcolor(lightmagenta) ; 

writeln;writeln('  SYSTEM  EQUIVALENT  TEMPERATURE  = 

Teq.dnl  :0;2,'  deg  K'); 

writelnC  G/T  FIGURE  OF  MERIT  =  ’  .GT.dnl  :0: 2,’  dB'); 

end; 

end; 

wait; 

textco1or( white) ; 
end; 

(**) 

(*  transmitter  analysis  *) 
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procedure  transmitter; 

var  power  (*  transmitter  power  *) , 

loss  (*  transmitter  circuit  losses  *):real ; 
begin 

textbackground(blue) ; 
clrscr; 

wr1telnjwrite1n('T  RAN  SMITTER  ANALYSI  S':60); 
gotoXY(l,5); 

if  uplink=true  then  writeln(title[5]:47)  else  writeln(title[6]:44); 

goto XY( 1,10); 

repeat 

writeln ; 

writeC  ENTER  OUTPUT  /AMPLIFIER  POWER  LEVEL  (W)  :*); 

Get Real ; 

power  :=realdata; 

if  not  (power>0)  then  InputError:=true; 
if  Input Error=t rue  then  noise; 
until  Input Error=  false; 
repeat 

writeln  * 

writeC  ’ENTER  TRANSMITTER  CIRCUIT  LOSSES  (-dB)  :'); 

Get  Real ; 
loss  :=nealdata ; 

if  not  (loss<=0)  then  Input  Error  :=true; 

If  Input Error=t rue  then  noise; 
until  Input Error=false; 

(*  store  input  data  and  perform  computations  *) 

with  data  do 

begin 

if  uplink=true  then 
begin 

tout. upl  :=power; 

Tx Circuit  Loss. upl  :=loss; 

EIRP.upl  :*isotropic(Pout.upl  .Tx&in.upl  .TxCircuit Loss. upl); 
textcolor(l ightmagenta) ; 
writeln  *r1teln('  EIRP  =  '  .EIRP.upl  :0:2,*  dB'); 
end  else 
begin 

tout.dnl  :=power; 

TxCIrcultloss.dnl  :=loss; 

EIRP.dnl  :=1sotropic(Pout.dnl  .TxGiin.dnl  .TxCircuit Loss .dnl); 

textcol or(  1 i ghtmagenta ) ; 

writeln jwr1teln('  EIRP  *  '  .EIRP.dnl :0:2,'  dB'); 

end; 

end; 

textcolor( white) ; 

wait; 

end; 

(**) 

(*  link  receiving  performance  analysis  *) 
procedure  Rec  tower; 
var  1  (*  counter  *)  integer; 
rate  (*  rain  rate  *)  :real ; 


sv 


<5? 


begin 

textbac  kg  round  (blue) ; 

for  1 :»1  to  2  do 

begin 

clrscr; 

gotoXY(l,5); 

writeln(t1t1e[6+i]:46) ; 

gotoXY(l.lO); 

vritelnC  1.  Negligible  (  <  .25  mq/hr  )'); 

writeln('  2.  Drizzle  (  .25  -  1.25  rmv'hr  )*); 

wrltelnC  3.  Light  tain  (  1.25  -  12.5  rnn/hr  )'); 

writeln('  4.  Medium  Rain  (  12.5  -  25  imVhr  )')*» 

writeln('  5.  Heavy  ffein  (  25  -  100  mn/hr  )'); 

wrltelnC  6.  Tropical  Dtwipour  (  >  100  mVhr  )'); 

wrltelnC  7.  SPECIFY  EXACT  RAIffiATE  (0-150  mn/hr  )'); 
repeat 

writeln*nriteln;writeC  ENTER  RAIM1A7E 

Get  Integer; 

choice  :=integerdata ; 

if  not  choice  in  [1..7]  then  IhputError:=true; 
if  Input Error^true  then  noise; 
until  InputEmor-false; 
case  choice  of 
l:rate:=0. 125; 

2:rate:=0.75; 

3:rate:=6.88; 

4:rate:*18.75; 

5:rate:=62.5; 

6:rate:=125.0; 

7:beg1n 

repeat 

wrltelnjwrlteC  ENTER  RAINRA1E  {m/hr)  :'); 

Get  Real ; 
rate: -real  data; 

if  not((rate>*0)  and  (rate<=150))  then  InputError:*true; 

If  Input  Error-true  then  noise; 
until  Input  Error- false; 

if  rate-0  then  if  1«1  then  data.RalnAtten.upl  :=0  else 
data.RainAtten.dnl  :=0; 

end; 

end; 

if  rate>0  then  with  data  do 

if  1»1  then  RainAtten.upl  :«nain(nate,freq.upl  .alpha.upl  .Ho.upl  Jat.upl ) 
else  RalnAtten.dnl  :«ra1n(rate,freq.dnl  .alpha.dnl  ,rt).dnl  .lat.dnl); 
if  1«1  then  data.ralnrate.upl  :*rate  else  data.rainrate.dnl  :«rate; 
end; 

data.A1rAtten.upl  :«a1r(data.freq.upl  /lata.alpha.upl) ; 
data.A1rAtten.dnl  :*a1r(data.fneq.dnl  .data.alpha.dnl); 

(*  confute  slant  range  .path  loss.recvd  power  .and  CNR  *) 

with  data  do 

begin 

Slanttonge.upl  :«(sqrt(sqr(Re)#sqr(s1n(alpha.up1*dr))+2*lte*data.altitude+ 
sqr(data.alt1tude))-Re*sin(a1pha.upl*dr))/le3; 


Slanttenge.dnl  :«(sqrt(sqr(Re)*sqr(sin(alpha.dnl*dr))+2*Re*data.altitude+ 
sqr(data.alt1tude))-Re*s1n(a1pha.dnl*dr)  )/le3; 

Path  Loss,  up! :=-36.6-20*logl0(S1anttenge.upl*le3*6.214e-4*freq.upl*le3) ; 
tethloss.dnl  :*-36.6-20*1ogl0(Slantfenge.dnl*le3*6.214e-4*freq.dnl*le3); 
Pr.upl  :*power(tout.upl  .TxCircuitLoss.upl  .TxGain.upl  .TxEfficiency.upl , 
Path  Loss,  upl  .AirAtten.upl  ,Ra1n/>tten.upl  ,ftc  Efficiency. upl  ,te&in.upl , 
RcCIrcultLoss.upl ) ; 

fY.dnl  :=power(Bout.dnl  .TxCincuitLoss.dnl  .Txteln.dnl  .TxEfficiency.dnl , 
PathLoss.dnl  .AlrAtten.dnl  .RalnAtten.dnl  .ffcEfficiency.dnl  .RcGain.dnl , 
foCircuitloss.dnl); 

Cl«t.upl  :=cnratio(Pr.upl  .Teq.upl ) ; 

CNR.dn1  :=cnratio(Pr.dnl  Jeq.dnl) ; 

end; 

end; 

(**) 

(*  link  transmission  performance  analysis  *) 
procedure  bltrate; 
var  rate  (*  cbtarate  *) , 
err  (*  bit  error  rate  *) , 
nun  (*  runber  of  signalling  levels  *);real; 
modulation  (*  modulation  type  *):integer; 
begin 
clrscr; 
gotoXY(l,5); 

if  uplink=true  then  writeln(t1tle[l]:50)  else  write1n(title[3]:50); 

gotoXY(l.lO); 

repeat 

write(‘  ENTER  DESIRED  DATA  RATE  (fb/sec)  :'); 

Get Real ; 
rate  :=realdata ; 

if  rate<le-3  then  IhputError:*true; 
if  Input Error*true  then  noise; 
until  InputError=fal$e; 
repeat 

vriteln  >r1teln ; 
wrltelnC  1.  BPSK* ) ; 

vr1teln(‘  2.  qPSK'); 

writeln('  3.  OqPSK'); 

writeln('  4.  MSK'); 

wrltelnC  5.  M-PSK'); 

wrltelnC  6.  Coherent  FSK'); 

wrltelnC  7.  Non -Coherent  FSK'); 

wrltelnC  8.  M-FSK' ) ; 

wrltelnC  9.  ASK'); 

wrltelnC  10.  M-ASK'); 

wrltelnC  11.  DPSK'); 

uriteln  • 

writeC  ’ENTER  MODULATION  TYPE  :'); 

Get  Integer; 

modulation  :»1ntegerdata ; 

If  not  (modulation  In  [1..11])  then  lnputError:«true; 

If  Input Error»t rue  then  noise; 
until  InputError»false; 
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nun  :=2.0; 

case  modulation  of 
5,8,10:begin 
repeat 

writeC  ENTER  M 
Get Real ; 
nun:=realdata; 

if  ((nun<2.0)  or  (frac(nun)OO))  then  Input  Error: -true; 
if  Input Error=t rue  then  noise; 
until  In  put  Error=  false; 
end; 

end; 

repeat 

writeC  ENTER  DESIRED  BER  :'); 

Get  Real; 
err:=realdata; 

if  not((err>0)  and  (erKl.O))  then  Input  Error: -true; 
if  Input Error=true  then  noise; 
until  Input  Error-false; 

(*  store  input  data  and  perform  computations  *) 

with  data  do  if  uplink-true  then 

begin 

datarate.upl  :=rate ; 
modtype.upl :=int(modulat1on) ; 

M.upl  :-nun; 

BER.upl  :=err; 

/Wail  B>fo.upl  :*Ebto/Wail  (Pr.upl  .Teq.upl  .datarate.upl ) ; 

ReqEbfo.upl  :=EbNjReq(trunc(modt^pe.upl )  ,M.upl  ,BER.upl ) ; 
margin.upl  :=AvailEbfo.upl-Req£bK>.upl ; 
ftaxDatarate.upl  :^naxrate(fV.up1  .Teq.upl  .ReqEbf'b.upl ) ; 
capac  ity  ( bandwi  dth  .upl  .Pr.upl  .Teq.upl  /texC&tarate.upl  ,M.upl ) ; 
textcolor(l  ightmagenta) ; 

writelnjwritelnC  Available  B>/No  3  '  ,/VailEbfo.upl  :0:2,'  dB'); 
writeln('  Required  Eb/lt>  »  '  .ReqEbffc.upl  :0:2,'  dB'); 
writeln('  tex  Data  tote  =  limit :0:2,'  H)/sec'); 
if  shannon  flag-true  then 
begin 

noise; 

writeln('***  SHWICN  RATE  USED  ***'); 
ffex  Datarate.upl  :*1  imit ; 

end; 

if  nyquistflag-true  then 
begin 

noise; 

writeln('***  NYgUIST  RATE  USED  ***'); 

Ftox  Datarate.upl  :-l  imit ; 

end; 
end  else 
begin 

datarate.dnl  :«rate ; 
modtype.dnl ;«int(modulat1on); 

M.dnl  :*run; 

BER.dnl  ;-err; 
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/Wall  a>fo.dil  :=a>fo/Vail  (Pr.dnl  .Teq.dnl , datarate.dll ) ; 
ReqB>to.dn1  :*Q>Nt>RBq(trunc(modt.ype.4il)  ,M.dnl  ,BER.dnl ) ; 
margln.dnl  :=/Wa11  Ebfo.dnl-ReqEbito.Ail ; 
texCbtarate.dnl  :=raxrate(Pr.dnl  .Teq.dnl  .ReqEbto.dnl); 
capacity(  bandwidth. dnl  .Rr.dnl  .Teq.dnl  ^bxUatarate.dnl  «M.dn1 ) ; 
textcolor(lightmagenta); 

writelnjwriteln('  Available  Eb/No  *  '  •A/all Bbhb.dnl  :0:2,*  dB'); 
writeln(*  Required  Ba/No  *  '  .ReqEbfto.dnl  :0:2,‘  dB* ) ; 
writeln('  Muc  feta  Rate  =  ' .limit :0:2,'  Ha/sec' ); 
if  shannonflag=true  then 
begin 

noise; 

writeln('***  SHAWCN  RATE  USED  ***'); 
ffex  Data rate. dnl :=1 imit ; 

end; 

if  qyquistflag=true  then 
begin 

noise; 

writelnC***  NYQUIST  RATE  USED  ***' ); 
ftex  Data  rate. dnl  :*1  imit ; 

end; 

end; 

textcolor( white) ; 

if  upl1nk=true  then  wait; 

end; 
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The  product  of  this  research  effort  is  a  comprehensive 
and  integrated  software  package  for  the  analysis  of  digital 
satellite  links.  The  system  is  designed  to  run  on  the  IBM 
PC,  and  should  also  run  on  most  compatibles. 

The  system  performs  three  basic  classes  of  functions: 
satellite  orbital  analysis,  antenna  gain  pattern  plotting, 
and  link  analysis.  The  first  class  of  functions  includes 
the  computation  of  such  quantities  as  velocity,  orbital 
period,  and  coverage  area  for  satellites  in  circular  and 
elliptical  orbits.  The  second  class  of  functions  is 
concerned  with  plotting  the  gain  patterns  for  horns, 
helixes,  parabolic  reflectors,  and  phased  arrays  of 
dipoles.  The  last  class  of  functions  represents  the  major 
thrust  of  the  system,  and  entails  computing  such  items  as 
the  G/T  figure  of  merit,  received  useful  power,  carrier-to- 
noise  ratio,  bit  error  rate,  maximum  data  rate,  and  power 
margin.  Inherent  within  this  class  are  mathematical  models 
for  computing  the  attenuation  due  to  rainfall  and 
atmospheric  absorption . 

The  link  budget  itself  appears  as  a  color-coded  display 
with  two  columns:  one  for  the  uplink  path,  and  one  for  the 
downlink  path.  The  user  also  has  the  capability  to  change 
certain  key  inputs,  and  then  have  the  system  automatically 
recompute  the  entire  link  budget  with  the  modified  data. 
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